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Foreword 


T his volume is one of a series which summarize the 
progress made during the period 1958 through 1964 
in discipline areas covered by the Space Science and 
Applications Program of the United States. In this 
way, the contribution made by the National Aero- 
nautics and Space Administration is highlighted against 
the background of overall progress in each discipline. 
Succeeding issues will document the results from later 
years. 

The initial issue of this series appears in 10 volumes 
(NASA Special Publications 91 to 100) which describe 
the achievements in the following areas; Astronomy, 
Bioscience, Communications and Navigation, Geodesy, 
Ionospheres and Radio Physics, Meteorology, Particles 
and Fields, Planetary Atmospheres, Planetology, and 
Solar Physics. 

Although we do not here attempt to name those who 
have contributed to our program during these first 6 
years, both in the experimental and theoretical research 
and in the analysis, compilation, and reporting of 
results, nevertheless we wish to acknowledge all the 
contributions to a very fruitful program in which this 
country may take justifiable pride. 

Homer E. New'ell 
Associate Administrator for 
Space Science and Applications, NASA 
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Preface 


G eodesy is the study of the size and shape of the Earth 
and of locations on the Earth. In its crudest forms, it 
is, like astronomy, almost as old as history. The beginning of 
the era of space geodesy may be set in 1958 with the announce- 
ment, based on the analysis of observations of Vanguard I 
(1958/3), that the flattening of the Earth’s poles is significantly 
smaller than had been derived from teiTestrial geodesy. This 
result implies the possible existence of considerable stress dif- 
ferences in the Earth’s interior which may be supported by 
internal mechanical strength or by convection currents. 

The first definite evidence that the Earth’s gravitational 
field was irregular was derived from observations of several 
satellites early in 1959. These observations and analyses 
showed that the Northern Hemisphere of the Earth contains 
slightly more material than the Southern Hemisphere. There- 
fore the equipotential surface (i.e., the fictitious surface on 
which the pull of gravity is the same everywhere) is farther 
from the Equator at the North Pole than it is at the South Pole. 
Since water follows such an equipxjtential surface, the oceans 
define a pear-shaped Earth. 

This result was quickly followed by further analyses of the 
Earth’s gravitational field as it affects the orbits of satellites. 
Mathematicians find it convenient to describe this field in 
terms of components which vary with latitude and components 
which vary with longitude. The latter are somewhat harder 
to determine, since they affect satellite orbits in much the 
same way as atmosplieric drag, light pressure, and other 
disturbing factors. . 

The increasingly accurate description of the Earth’s gravi- 
tational field will eventually lead to a more precise under- 
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standing of the internal structure of the Earth. If, as has beerf* 
proposed, the irregularities in the terrestrial gravity field are 
correlated with the rate of heat conduction from the interior 
of the Earth, we may have an important clue to physical con- 
ditions in various parts of the Earth’s mantle. Moreover, a 
more detailed analysis of these results should enhance our 
understanding of seismic and volcanic activity. On a more 
direct and practical level, most oceanic navigation is based on 
a knowledge of local gravity and of the gravitational shape of 
the Earth. Thus, an improved knowledge of the terrestrial 
gravity field can lead to improved navigation accuracy in 
areas not reached by modern electronic navigation beacons. 

Not only can a detailed knowledge of the Earth’s gravita- 
tional field be derived from the analysis of orbits of satellites 
but the process can and must be reversed to provide us more 
accurate predictions for satellite- and space-probe trajectories. 
Such knowledge is vital, for example, for proper direction of 
manned satellites which will rendezvous with other satellites 
launched earlier. 

Four satellites have been of outstanding importance for 
geodetic studies. The two Echo balloon satellites provided a 
readily visible target which could be photographed against a 
star background from widely separated sites on the surface of 
the Earth. By photographing the satellite simultaneously 
from several sites, the lengths and relative orientations of the 
baselines between the sites can be derived by standard triangu- 
lation techniques. The first satellite designed and launched 
specifically for geodetic purposes was the Anna satellite 
launched by DOD in 1962. This satellite carried a high- 
intensity xenon lamp which could be commanded to provide a 
sequence of five closely spaced flashes. These flashes pro- 
vided starlike images which were observed against a back- 
ground of stars. The images could be measured with pre- 
cision, and the active nature of the satellite’s flashes insured 
that observations were indeed simultaneous. The Anna 
satellite also carried a highly stable transmitter for Doppler 
studies of satellite motion. The SECOR electronic ranging 
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system did not work satisfactorily on this satellite, but has since 
'been successfully tested on other satellites. Anna provided a 
useful and encouraging test of both the flashing-light and 
Doppler-beacon systems, and many of the lessons learned from 
Anna have been used in the design of the Geos satellite, whose 
launch was planned for 1965. 

Another satellite of special geodetic usefulness, Syncom, was, 
like the Echos, launched as a communication satellite. Be- 
cause of its special orbit which keeps it essentially fixed with 
respect to a given longitude, it is particularly sensitive to cer- 
tain terms in the Earth’s gravitationail field which are difficult 
to derive from an observation of satellites in lower orbits. It 
is encouraging to find that the analysis of the orbital motions 
of Syncom substantially confirmed analyses based on other 
satellites. Much progress has been made in satellite geodesy 
during the years 1958 to 1964, although it is likely that sut> 
stantial improvements will follow the launch of the Geos and 
Pageos satellites. 

The tracking of satellites, and particularly of space probes, 
has resulted in improvements in both terrestrial and other 
planetary constants sufficientiy great that the International 
Astronomical Union (lAU) undertook a review of the con- 
stants used in the various national ephemerides as well as in 
computing the orbits of satellites and space probes. On the 
basis of the detailed analysis of all the evidence to date, the 
lAU adopted, at its meeting in Hamburg in August 1964, a 
new set of astronomical constants. This was the first major 
change in the astronomical constants in this century, and will 
significandy improve astronomical predictions in the future. 

The use of Geos and further analysis of Syncom data should 
substantially improve our knowledge of the Earth’s gravita- 
tional field, but the most outstanding accomplishments ex- 
pected in satellite geodesy within the relatively near future 
oTC 1X1 geometrical geodesy (i.e., the accurate location of various 
sites on the same Earth-centered reference frame), a field 
which has barely been touched by satellite techniques. Ob- 
servations of both Geos and Pageos will permit, for the first 
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time, accurate mapping throughout the world on a single 
coordinate system, and a determination of the relative loca-* 
tions and orientations of the various geodetic datums which in 
the past have been derived separately for each country or 
region with little or no possibility of interconnection. This 
unification will, of course, lead to appreciably more accurate 
maps and better defined distances between pgints on different 
continents, and to accurate location of isolated islands. These 
improved maps will also permit more accurate tracking and 
guidance of interplanetary probes, since tracking stations will 
be better located with respect to the center of the Earth. 

Although NASA has provided the majority of the satellites 
used for geodetic studies, and much of the tracking data on 
which the analyses have been based, there has been substantial 
participation by many groups throughout the world. Prob- 
ably the work done directly by NASA or with NASA support 
represents only about one-quarter of the activity in the field. 

So far, almost all the geodetic studies have referred to the 
Earth. However, the lunar-orbiter program will extend the 
techniques of satellite geodesy to the Moon. Mariner II 
already has provided an improved value for the mass of Venus. 
Improved techniques for observing the Martian satellites, 
stimulated by NASA interest in that planet, will improve our 
knowledge of the Martian gravity field. Geodetic studies of 
planets other than the Earth by measurement of and from 
planetary orbiters will be very important in the next 
decade. 
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Introduction 


F or the purpose of the present review, satellite geodesy 
is defined as those areas of dynamical and geometrical 
geodesy to which satellite techniques have contributed exten- 
sively. This review indicates the contributions to our under- 
standing of the Earth’s gravitational field which have come 
from analyses of satellite motions, and also the contributions 
to the lesser but potentially more important use of satellites 
for geodetic triangulation. There are many related areas of 
geodesy and geophysics which have important interfaces with 
satellite geodesy; these have not been included in this review. 

Many authors have reviewed various aspects of the progress 
in the field of satellite geodesy. Rather than rewrite their 
material, representative publications have been collected and 
combined to form this report. Except where otherwise noted, 
the editing has been limited to the occasional selection of 
appropriate p>ortions from longer articles. It should be em- 
phasized that while an attempt has been made to include 
appropriate articles, in many cases, other equally valid choices 
might have been included instead. Some articles are included 
not because they give significant results in themselves, but 
because they describe an important technique. Together, 
it is believed that these articles give a reasonably complete 
review of the progress and activity in the field of satellite 
geodesy during the period 1958 to 1964. 

More specifically, the following material is covered in the 
report: 

(1) The state of knowledge in geodesy at the time of NASA’s 
creation. 

(2) Early results in satellite geodesy. 

(3) Dynamical geodesy based on optical observations. This 
is the first technique employed in satellite geodesy and has 
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produced the most extensive results to date. Three modern 
summaries include recent determinations of both zonal and * 
tesseral harmonics. 

(4) Dynamical geodesy based on nonoptical tracking. 
Doppler tracking was applied to determining the orbit of the 
first sputnik. It has recently proven an important geodetic 
tool. The SECOR ranging system and the various range-and- 
range-rate systems are beginning to provide geodetic data. A 
particularly important contribution has resulted from the use 
of a range-and-range-rate system to track Syncom. 

(5) Astronomical constants: This important area is included 
because of the role of satellite geodesy in determining the ter- 
restrial constants and because of the closely related techniques 
used for determining other constants. 

(6) Geometrical geodesy using optical techniques. 

(7) Geometrical geodesy using nonoptical techniques. 

(8) Special satellites (Anna and Echo) for geodesy and tests 
of techniques for using these satellites. 

No attempt has been made to include the theory used in 
deriving the terrestrial gravity field from observations of arti- 
ficial satellite orbits. This has been extensively reviewed by 
Cook,* Kaula,^ Mueller,^ and others. The interested reader 
is referred to these authors both for the mathematical bases of 
the geodetic analyses and for more detailed background on the 
applications of the various mathematical techniques. 


Late in 1958, a conference on “Contemporary Geodesy" was 
held in Cambridge, Mass., to review the rapid impact of elec- 
tronic and satellite tracking techniques on this ancient field. To 
set the background for further discussion. Dr. Alwyn R. Robbins 
of Oxford University presented a review of the history and current 
status of geodesy. The following was originally published in 


^Cook, A. H.; Space Science Reviews, vol. 2, Sept. 1953, pp. 355-437. 

^Kaula, W. M.: Celestial Geodesy. In: Advances in Geophysics. Vol. 9, Academic 
Press, 1962, pp. 192-293. 

^Mueller, Ivan I.: Introduction to Satellite Geodesy. Frederick Ungar Pub. Co. 
(New York), 1964. 
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Contemporary Geodesy, Geophysical Monograph No. 4, American 
4 Geophysical Union of the National Academy of Sciences- 
National Research Council Publication No. 708, 1959. 

Evolution of the Geodetic Concept 

Atwyn R. Robbins 
Oxford University 
Oxford, England 

Geodesy is a very old profession. If you look in the Bible, I think 
the Book of Numbers or Deuteronomy where there is a list of curses, you 
will find one which in effect says; “Cursed be he who moves his neighbor’s 
boundary stone.” So it goes back some way. 

Next, the ancient Greeks thought the Earth was a plane supported by 
four elephants on the back of a turtle. Aristotle went a step further and 
said it was a sphere. Later Eratosthenes noticed that the Sun shone directly 
down a well at noon at the summer solstice; he observed the sun somewhere 
else at the same time, made a traverse by camel caravan, and computed 
the radius of the Earth. 

Then things stood still for a few centuries. With the coming of the 
telescope and the use of logarithms, triangulation was originated. Finally 
the size and shape of the Earth was measured by triangulating along 
meridians to compute the semiaxis and the flattening of the ellipsoid. 
One measurement appeared to show that the Earth was a prolate spheroid; 
this was disputed by Newton and others and then we had the famous 
French arcs in 1735 and 1736 which proved it was an oblate spheroid. 

In the nineteenth century, with the realization of the need for maps, 
many countries observed the national framework of triangulation and 
some of them determined their own spheroid, or figure of the Earth, that 
happened to fit their country best. That was all right but when com- 
munications improve, national barriers become meaningless and geodetic 
networks must be international. The spheroid that fits one country 
doesn’t necessarily fit others. 

During this century, the-International Union of Geodesy and Geophysics 
has given geodesy international recognition. During the last thirty years, 
especially since the last world war, these national triangulations have been 
linked more and more and many datums have been tied together by 
triangulation. 

When triangiilation over large areas is computed on an ellipsoid or 
spheroid whose size and shape are known, position of the sphci oid in rela- 
tion to the geoid must also be determined. A datum contains seven 
constants: two, the semiaxis and the flattening of the spheroid, and two 
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additional constants to make the minor axis parallel to the axis of rotation 
of the Earth. These last two you do not use per se but you use them, 
without noticing as it were, when computing geodetic azimuth from 
astronomical. As far as these four constants are concerned, you can 
compute on any spheroid you choose but you still will not necessarily be 
on the same datum. Finally, you define the latitude and longitude and 
geoid-spheroid separation at the origin. Now the datum is completely 
defined. Nothing else can be defined; everything else must be computed. 
So if you have two disconnected triangulation systems on the same 
spheroid, they are still on different datums in that they have different origins. 

The definition of latitude and longitude and geoid-spheroid separation 
at the origin is completely arbitrary. You can assume that the separation 
is zero and that the geodetic latitude and longitude are the same as the 
astronomical. If you do that and if you happen to be in an unlucky spot 
where the geoid rises or falls slightly, then as the network extends hundreds 
of miles, this tilt will become more pronounced and the separation of the 
two surfaces will increase. Generally, you will reduce your bases to mean 
sea level but you should reduce them to the spheroid. If you have enough 
deviations of the vertical you can compute along section lines and calculate 
the separation and its effect on scale. But one seldom has enough 
information. 

There is, however, one way of going about it: You can compute devia- 
tions of the vertical on one world datum if you have enough information 
on the intensity of gravity over the world. However, there is not enough 
gravity information so there will probably be some residual errors left 
in computation because of insufficient data. Perhaps some will disagree 
with that statement. 

Be that as it may, it is important to recognize that deviations of the 
vertical from Stokes’ theorem are on one datum, and any others computed 
on other datums are different. The two cannot agree except by chance. 
So we have a multiplicity of datums. The task of the geodesist is to 
reduce these and combine them into one world geodetic datum. You can 
do it by having more observations of gravity and so on, or you can make 
intercontinental ties between triangulation systems. Then you have the 
problem of computing the separation of the geoid and spheroid across the 
sea gaps. One way is to use gravity and interpolate. Ways of doing it 
are now being studied; some research is being done on that at the moment. 
So it does not really matter what datum you have, as long as you have one 
which fits reasonably well and as long as you have enough observations. 
It is the lack of sufficient observational information that is holding things 
up to some extent at the moment. The objective is to have a world 
datum and to portray the geoid on it. 

We have come a long way since the introduction of the telescope. I do 
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not recall the date of the early triangulation in Great Britain but I remem- 
* her their geodetic theodolite had a sixty-inch circle and they had to put it 
on top of St. Paul’s Cathedral on a scaffolding. Nowadays one can get 
better results with a five inch. Of course we also have shoran, and more 
recently still the satellite. 

On the gravity side, for the pendulum we have come up with more 
accurate timing devices. The gravimeter is being improved and new 
types are being developed which can be used aboard surface ships as 
well as under water. 

Finally we come to the Earth satellites which are, perhaps, controversial. 
How much can we get out of them? I would personally like to see it the 
other way around. How much information can we geodesists give to the 
physicists? If we know the gravity on the Earth and then tell the physicist 
what gravity is doing to the satellite, the physicist can determine what 
the effect of atmospheric drag is and so on. That, again, is the reverse 
of what many people are thinking. The other way around is to try to 
find out from the physicist what the drag is doing, whence to determine 
gravity all over the Earth. It all depends on the relative sizes of the 
effects we get from one source or the other. 
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Part 1 

Early Results in Satellite Geodesy 


The effect of the Earth’s gravitational field on the 
motion of near-Earth satellites is so marked that even 
comparatively crude tracking techniques sufficed to pro- 
duce the first results in satellite geodesy. An improved 
value for the flattening of the Earth and a lack of sym- 
metry between Northern and Southern Hemispheres 
were reported in 1958 and early in 1959, respectively. 
The following are early papers in this field. 
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The following •paper •was originally published in Nature, 
vol. 18S, 1958, pp. 198-199. 

N66 37347 

Motion oj the Nodal Line o] the Second Rus- 
sian Earth Satellite (1957^) and Flatten- 
ing of the Earth 

E. Buchar 

Observatory of the Technical University 
Karlovo ndm. 13, Prague II 


D espite the fact that the ohservation of an artificial satellite can 
be carried out only on a small region of the Earth’s surface, we can 
use these measurements for determining the geocentrical co-ordinates of 
the satellite, if we know the approximate value of its distance from the 
Earth’s centre. The most suitable observations are those in the neigh- 
bourhood of the zenith. If we know the inclination of the orbit, we can 
then determine also the position of the orbital node. 

In order to determine the motion of the orbital node, we used 33 
visual observations, made between December 7, 1957, and March 21, 
1958, at various places in Czechoslovakia. By adjustment of the posi- 
tions of the node grouped about five normal positions, the following 
expression for the right ascension of the ascending node was derived : 

a = 256.59° ± 0.20° - (2.9007° ± 0.0046°) (<- <o) 

- (0. 00234° ± 0.00012°) (< - <o) == 

Epoch to = January 22.0, 1958, U.T. Time t is expressed in days. The 
mean errors were computed from the deviations of the normal places. 
The daily motion of the node is then given by: 

^ - 2.9007 ± 0.0046° - (0.00468° + 0.00024°) {t - to) 


It is well known that the theoretical value of this motion is given by 
the approximate equation: 


dQ 

dt 


3fc(C-A) 

2ma^ 


(1 +2e^) cos i 
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where C — A is the difference between the moments of inertia of the ter- 
restrial spheroid, m the mass of the Earth and k the gravitational con- 
stant. The symbol a stands for the major semi-axis of the orbit, expressed 
in units of the equatorial radius of the Earth, e and i being the eccentricity 
and inclination of the orbit, respectively. 

It is obvious that, by using this equation, we can derive from the 
observed motion of the node the quantity K={C — A)/m. If, for the 
time t = to, we assume a= 1.1127 + 0.0003 and for e and i the values 
published by D. G. King-Hele,^ namely, e = 0.0731 +0.0005, f = 65.29° + 
0.03°, we then arrive at: 

K = 0.0010856 + 0.0000024 

The mean error was computed as the total effect of the errors of the 
basic data. The oblateness of the Earth, a, can be computed from the 
expression: 

a = ^K+^+l(SK^+Kh-h^) 

where h is the ratio of the centrifugal force to the value of gravity at 
the equator. 

The reciprocal value of the Earth’s flattening is then: 

- = 297.7±0.3 

a 

If we use four more values for a and e determined^ between February 
10 and March 3, 1958, we arrive finally at the result: 

E: = = 0.0010883 + 0.0000014 

m 

i = 297.4 + 0.2 

a 

In these numerical results approximate account was taken of the 
second term of the disturbing function. 

These preliminary results show that satellite observations of greater 
precision will give us the possibility of determining correctly the oblate- 
ness of the Earth. 


* King-Hele, D. G., Nature, 181, 738 (1958). 

’ Smithsonian Astrophys. Obs. Circulars (Cambridge, U.S.A.). 
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The following paper was originally published in Science, 
vol. 129, No. SSJfS, Feb. 27, 1959, p. 565. 

N66 37348 

• — m 

Vanguard Measurements Give Pear-Shaped 
Component oj Earth* s Figure 

J. A. O’Keefe, Ann Eckees, and R. K. Squires 
Theoretical Division 

National Aeronautics and Space Administration 
Washington, D.C. 

T he determination of the orbit of the Vanguard satellite, 1958|32, 
has revealed the existence of periodic variations in the eccentricity 
of that satellite (i). Our calculations indicate that the periodic changes 
in eccentricity can be explained by the presence of a third zonal harmonic 
in the earth’s gravitational field. The third zonal harmonic modifies 
the geoid toward the shape of a pear. In the present case, the stem of 
the pear is up — that is, at the North Pole. According to our analysis, 
the amplitude of the third zonal harmonic is 0.0047 cm/sec^ in the surface 
acceleration of gravity, or 15 meters of undulation in the geoid. 

Figure 1 shows the observed variation in eccentricity. The p>eriod of 
the variation in eccentricity is 80 days, approximately equal to the 
period of revolution of the lines of apsides. The eccentricity is a maxi- 
mum w’hen the perigee is in the Northern Hemisphere. The amplitude 

.19100 
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^ .19000 

u 

I .18950 
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.18850 
.18800 
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Days since launch 

Figure 1 . — Eccentricity of satellite 1958/3j (Vanguard). 
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of the variation is 0.00042 + 0.00003. Similar perturbations may exist 
in the angle of inclination of the orbit, although the data for them are 
much less accurate. No perturbations of this magnitude appear to 
exist in the semimajor axis. 

In principle, the perturbation might be caused by both odd and even 
harmonics. However, the even harmonics can be excluded because the 
observed effect has opposite signs in the Northern and Southern hemi- 
spheres. Furthermore, we can also exclude tesseral harmonics (those 
which depend on longitude as well as latitude) because these also are the 
same in the Northern and Southern hemispheres, apart from a shift in 
longitude. We are left with the zonal harmonics (those which depend 
only on latitude) of odd degree. 

Of the odd zonal harmonics, the first degree is forbidden; and those of 
higher degree are unlikely to have a large effect because they die out 
inversely as the (n-|-l) power of the distance. The effect is therefore 
due mostly to the third zonal harmonic, with a possible contribution 
from the fifth. 

Accordingly, a calculation was made of the effect of the third zonal 
harmonic on the orbit elements of 1958/32, by methods developed by 
O’Keefe and Batchlor (2). In the resultant expression for the eccen- 
tricity, the dominant terms were those whose argument was the mean 
motion of perigee. These were larger than the others by a factor of 10®. 
Keeping only the large terms, we find 


,3, 

c — eo+oi4.3.o — 

2 na® 


sin i 


'■ (l - I sin® 


Sin CO 


( 1 ) 


where Aa.o represents the coefficient of the third zonal harmonic in the 
notation of Jeffreys (S), n is the orbital mean motion and n' is the mean 
motion of the perigee, e is the .eccentricity and Co the mean eccentricity, 
i is the angle of inclination, oj is the argument of perigee, and a is the 
semimajor axis. 

Setting in the constants of the orbit and the observed amplitude of e, 
we find 


A3,o=( 2.5±0.2)X10®» (2a) 

in meter-second units. Utilizing the relation given by Jeffreys, 

A _ 


(where A„,, is the coefficient of the disturbing potential, g„,, is the 
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acceleration of gravity at the surface of the earth, and c is the earth’s 
equatorial radius), we find that the third zonal harmonic of gravity at 
the earth’s surface, in milligals, is 


Sr3.o = 4.7±0.4 (2b) 

Equation 2 is relevant to what Vening Meinesz (4) and Heiskanen call 
the “basic hypothesis of geodesy.” These authors assume that the 
earth’s gravitational field is very nearly that of a fluid in equilibrium. 
They consider that the deviations from such an ellipsoid, in any given 
area, do not exceed about 30 milligal-megameter units — ^that is, they 
assume that one will not find deviations of more than 30 milligals over 
an area of 1000 kilometers on a side, or deviations of more than 3 milligals 
in an area 3000 kilometers on a side. 

Our determination of the third-degree zonal harmonic shows that the 
h 3 q)othesis of Vening Meinesz and Heiskanen is not justified ; for example, 
each of the polar areas has a value of about 120 milligal-megameters, 
and each of the equatorial belts a value more than twice as great. 

The presence of a third harmonic of the amplitude (2) indicates a very 
substantial load on the surface of the earth. Following the arguments 
of Jeffreys, w’e may calculate the values of this load and the minimum 
stress required in the interior to support it. We find a crustal load of 
2 X 10’ dy/cm’. W’e can choose between assuming that stresses of approxi- 
mately this order of magnitude exist down to the core of the earth, or 
that stresses of about 4 times that amount exist in the uppermost 700 
kilometers only ( 3 , p. 199). These stresses must be supported either by 
a mechanical strength larger than that usually assumed for the interior 
of the earth or by large-scale convection currents in the mantle (5). 
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Part 2 

Derivation of the Earth’s Gravity Field From Optical 
Photographs of Satellites 


The network of Baker-Nunn wide-angle cameras 
which was established during the International Geo- 
physical Year has been tracking American and Russian 
satellites by photographing them at dawn and at dusk 
against a star background. Drive motors permit the 
camera to compensate for the motion of the satellite, thus 
allowing the photography of very faint satellites. In 
this mode, stellar images are short trails. So that they 
can be measured at a given instant of time, these trails 
are chopped by an accurately timed shutter. Although 
various techniques can be and have been used to track 
satellites and geodetic information has been derived from 
most of these techniques, the Baker-Nunn cameras have 
undoubtedly produced the most important body of data 
on which such analyses have been based. The follow- 
ing are the results of three recent analyses of these optical 
data. The section by Kaula is a coalition of two of his 
papers which appeared in the Journal of Geophysical 
Research in 1963. 
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papers originally pvhlished in Journal of Geophysical 
Research, vol. 68, 1963, pp. 473-484 and 5183-5190. 


Current Knowledge oj the Earth* s Gravita- 
tional Field From Optical Observations 

W. M. Kaula 

Goddard Space Flight Center, NASA 

T hts paper describes the analysis of Baker-Nunn camera observa- 
tions listed by Veis et al. [1962] and Haramundanis [1962] by the 
methods described by Kaula [1961a, b; 1962a, b]. An IBM 7090 com- 
puter was used. Solutions were made for all geodetic and gravitational 
parameters estimated to have effects of more than ±20 meters on 
satellite orbits. The intent of the analysis was to apply all devices 
short of allowing for covariance of observations at different times. 
This intent resulted in programs so complicated that most of the time 
spent on the work was consumed by purely computational difficulties. 

OBSERVATIONS 

The Baker-Nunn camera system, its accuracy, and experience in its 
operation by the Smithsonian Institution Astrophysical Observatory are 
described by Henize [1960], Lassovszky [1961], Weston [1960], and Veis 
and Whipple [1961]. That the random error of the plate measurements 
is of the order of ± 2 " has been confirmed in this analysis by the accuracy 
with which a line can be fitted to plotted residuals with respect to an 
orbit of observations close together in the same pass. Since the signifi- 
cant timing error is virtually constant throughout a pass, no such test of 
timing error is possible because of the dominant effect of drag error in 
the orbit. 

The precisely reduced Baker-Nunn camera observations of 1959a, 
1959 j;, and 1960i2 from launch until the end of 1961, of 19615, from 
launch until the middle of 1961, and of 1961a5, in the spring of 1962 
were analyzed. The observations through mid-1961 have been pub- 
lished in the catalogs compiled by Veis et al. [1961-1962] and are referred 
to in the 1950 mean positions of the stellar catalog. For this analysis, 
the epoch of the right ascension and declination was updated to the 
epoch of the orbital arc fitted to the observations, taking into account 
precession plus nutational terms of more than 0.25' amplitude — i.e., the 


17 



SATELLITE GEODESY 


18.6 year and semiannual terms. All times are given for the observa- ■. 
tions, which are treated as equivalent to ephemeris time. A small 
correction was applied in calculating Greenwich sidereal times to allow 
for the precession and nutation between the epoch of the orbital arc and 
the instant of observation. 

The above-mentioned ± 2" accuracy of fitting of a line to residuals is 
appreciably smaller than the residuals themselves, which indicates that 
extra observations within a pass did not add extra weight to the analysis 
of the orbit. Hence, to conserve computer time and to avoid over- 
weighting certain passes, observations were omitted which were neither 
terminal observations of a pass nor observations interior to a pass at 
intervals of 2 minutes or more. 

The final rejection criterion applied was to pmit observations on 
days of appreciable atmospheric disturbance, as measured by the geo- 
magnetic index A^. For the 1960t2 analysis observations were omitted 
on days for which Ap exceeded 50; for 1959ai and 1959 tj, when Ap ex- 
ceeded 70. In some cases additional observations on adjacent days 
were omitted in order that an orbital arc would not bridge across days 
of high Ap index. 

The principal defect in the observations is, of course, their poor dis- 
tribution — due to the dependence on reflected sunlight, and to the 
limited number of tracking stations (twelve). The number of observa- 
tions of each satellite used is given in Table 1. 


Table 1 . — Satellite Orbit Specifications 



Satellite 

1959a 1 

1959>) 

1960i2 

19615i 

1961a«i 

Epoch _ 

1959 Feb. 

1959 Sept. 

1960 Sept. 

1961 Feb. 

1962 Mar. 


28.5 

28.5 

22.0 

20.0 

8.5 

Semimajor axis 

1.304585 

1.334500 

1.250057 

1.252779 

1.568136 

Eccentricity 

0.16582 

0.19008 

0.01146 

0.12135 

0.01197 

Inclination 

0.57381 

0.58212 

0.82434 

0.67835 

1.67316 

Argument of perigee 

3.36062 

3.20403 

2.26377 

2.02733 

4.28853 

Longitude of node _ 

2.52442 

3.48304 

2.28139 

2.76786 

5.7i;i36 

Mean anomaly 

6.00463 

3.82408 

2.72868 

5.96587 

1.51124 

Perigee motion /day 

-t-0.09181 

-1-0.08501 

-1-0.05186 

-hO.08315 

-0.01733 

Node motion/day _ _ 

-0.06108 

-0.05712 

-0.05413 

-0.06347 

-t-0.00367 

Max. A/m, cmVg-- 

0.21 

0.27 

0.27 

15.9 

0.08 

Min. A/m, cm^/g._ 

0.21 

0.04 

0.08 

15.9 

0.02 

Perigee height, km. 

560 

510 

1500 

640 

3500 

Number of days 

1032 

792 

480 

150 

54 

Number of observa- 






tions . _ 

3513 

3034 

2502 

1395 

552 
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Geometry . — The observation equation used was expressed in terms of 
the meridian and prime vertical components of the plate measurement, 
assuming that the satellite was on the camera axis [Kaula, 1961a, 1962c]. 
It consists of the first two rows of the matrix equation 



+Rbx{Rxqq 


+Ciede+Ci,wn d<— Ra(— fl)duo)/r (1) 


In (1), a is the declination, a the right ascension, r the camera-satellite 
range. The first two row’s of (b/r)*^ are zero if the observed 5, a are 
used in Rbx, the matrix w’hich rotates the inertial coordinate system to 
a rectangular system with the 3-axis coinciding with the camera-satellite 
line and the 1-axis coinciding with the meridian, q is the satellite 
position in orbit-referred coordinates, with the 1-axis toward osculating 
perigee and the 3-axis normal to the osculating orbit ; R,, is the rotation 
from orbit-referred to inertial coordinates; is a 3X6 matrix of partial 
derivatives of the inertial rectangular coordinates with respect to the 
osculating Keplerian elements, corrections to which are symbolized by 
de; is the row of C^e corresponding to the mean anomaly; n is the 
mean motion; dt is a correction to the time of observation; Rs ( — 0) is 
the geodetic to inertial rotation matrix, with the Greenwich sidereal 
time, 6, as argument ; and dUo is a vector of corrections to station position. 
(Derivations of all these v-ariables are given in equations (46), (47), and 
(52) to (60) of Kaula [1961a], or equations (3.1) to (3.8) and (3.11) to 
(3.15) of Kaula [1962c].) 

The partial derivatives in (1), 


c.={ 


( dS/dt ) 

(1,0,0) 

\a(a cos S)/dt) 

(0, 1,0/ 


xCxMfi/r 


( 2 ) 


were not actually used to determine timing corrections, but rather in 
three other ways: first, to apply a correction rC,/c for the time of travel 
of the signal, where c is the velocity of the light; second, to give lower 
weight to the along-track component than to the across-track component 
of the observation, by giving each observation a 2X2 covariance matrix: 

= (3) 


where a/ is the variance of the direction measurement, <r,^ is the variance 
of the timing, and the superscript T denotes transpose; and, third, to 
compute residuals in along-track and across-track components by applying 
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to (3, a cos S) residuals the rotation 

c,/Vc7^\c,/Vc7w 

i-c,/Vc7+^\c^/^/c?+c? 

where Ci, C 2 are the two elements of C,. 

Table 2 . — Tracking Station Data 


[In length units of 6.378165 meters] 


Station 

Lat. and 
Long., 
deg 

Datum 

Starting 

coordinates 

Pre- 

assigned 

Prelimi- 

nary 

solution 

Final 

solution 

Organ Pass 

32.4 

Am. 

-240778.9 

±3.0 

-14.8 

-02.8 ±0.9 


253.4 


-810109.7 

+ 3.9 

-05.6 

-03.8 ±1.8 




533234.2 

+ 3.1 

+03.3 

-00.2 ±1.0 

Arequipa 

-16.5 


304591.7 





288.5 


-909989.8 







-281725.5 




Curacao . . 

12.1 


353050.9 





291.2 


-912004.8 







+208082.5 




Jupiter 

27.0 


+ 153068.1 





279.8 


-878214.3 







+451581.1 




OUfantsfontein _ _ 

-26.0 

EAST 

792726.2 

+ 3.4 

+ 16.3 

+05.4 ±1.3 


28.3 


425915.7 

+ 2.9 

-17.8 

-08.0 ±1.4 




-435196.0 

+ 2.9 

+00.1 

+02.6 ± 2.5 

San Fernando _ _ 

36.5 


800487.5 





353.8 


- 87038.6 







591033.8 




Naini Tal. 

29.4 


159630.5 





79.5 


857808.2 







487532.6 




Shiraz 

29.6 


529444.8 





52.5 


690490.0 







491723.2 




Woomera 

-31.1 

Au. 

-624562.7 

+ 11.3 

+06.4 

-15.0 ±4.0 


136.8 


586884.9 

+ 14.5 

+ 14.5 

+04.4 ±8.3 




-513573.3 

+ 13.2 

+02.8 

+08.7 ±6.5 

Tokyo -- 

35.7 

JKM 

-618774.5 

±5.2 

+00.0 

-08.5 ±2.7 


139.5 


527787.3 

±6.9 

+ 15.2 

+06.8 ±3.4 




579917.5 

±5.2 

+01.9 

+00.9 ±2.4 

Villa Dolores 

-31.9 

Ar. 

357509.6 

±28.4 

+2.7.1 

+36.9 ±3.9 


294.9 


-770550.4 

±22.8 

-02.2 

+03.4 ±5.0 




-526083.5 

±26.2 

+01.2 

+03.3 ±6.6 

Maui 

20.7 

H 

-857008.8 

±21.7 

+01.8 

+01.5±6.1 


203.7 


-376954.1 

±35.8 

+28.5 

+ 14.1 ±7.9 




351587.3 

±37.1 

-53.3 

-50.1 ±4.7 







' 
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Since effects were investigated that were expected to be as small as 
• ± 20 meters, all stations were assumed to have position error, but those 
stations connected to the same geodetic system were assumed to shift 
together. Hence the twelve cameras were referred to six datums: four 
to the Americas (Am.) system; four to the Europe-Africa-Siberia-India 
(EASI) system; and one each to the Austraha (Au), Japan-Korea- 
Manchuria (JKM), Argentina (Ar) and Hawaii (H) systems. The 
initial station positions used are the solutions given in Table 2. Cor- 
rections for errors in the computed positions of three stations relative to 
the principal datums were provided by I. G. Izsak of the Smithsonian 
Institution Astrophysical Observatory. Corrections to coordinates u,, 
U 2 , and M3 in earth radii are listed in sequence (all values times 10~®): 

San Fernando.. -|-5.6 —5.0 —8.5 

Xaini Tal +2.7 —5.0 +4.9 

Curacao —1.0 0.0 +2.8 

For the Am., EASI, and JKM systems, the starting station positions 
were those obtained in the solution for a world geodetic system by 
Kaula [1961c]. For the Au, Ar, and H systems the positions calculated 
by Veis [1961b] were taken and shifted by placing tangent to the datum 
origin an ellipsoid of flattening 1/298.3 and an equatorial radius of 
6,378,165+A'^o meters, where No is the geoid height in the vicinity of the 
datum origin as given by Kaula [1961c]. The initial station positions 
are given in Table 1 in length units of 6.378165 meters referred to the n 
coordinate system, with axes toward 0°, 0°; 0°, 90° E; and 90° N, 
respectively. 

The datura shifts listed in Table 3 apply to the starting coordinates 
in column 4 of Table 2. 

Dynamics . — Variables in the obser\'ation equation (1) that are depend- 
ent on the djmamics of the satellite orbit are 

Riq = R3(— fl)Rl( — f)R3( — Mj) (5) 

! a(cos E — e)'l 

a\/l— sin E? (6) 

where E, a, e, i, w, and Q are the osculating eccentric anomaly, semi- 
major axis, eccentricity, inclination, argument of perigee, and longitude 
of the ascending node, respectively; and 


de — J deo'-hCe,dp,-|-Ceidpi-|-Ceddpd+Cepdpp (7) 

where Cg' denotes the elements of an intermediate orbit at epoch ; Po are 
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Table 3 . — Datum Shifts 

[In length units of 6.378165 meters] 


Datum 

Coor- 

dinate 

1959«i 

1959^ 

196012 

196l5i 

1961q:5i 

Weighted 

mean 

Americas . . 

Au \ 

-02.5 

-02.6 

-03.7 

-03.8 

-06.4 

-03.8 + 1.0 


AU 2 

-04.7 

-05.2 

-09.6 

-11.3 

-04.2 

-05.1+0.8 


Am 3 

-00.9 

-00.5 

+00.6 

-01.9 

-00.3 

-00.4+0.2 

Europe-Africa- 

Aut 

+06.5 

+07.3 

+06.6 

+ 11.6 

+04.6 

+05.8 + 0.7 

Siberia-India 

AU 2 

-07.8 

-07.8 

-09.3 

-04.1 

-10.2 

-08.9+0.5 


Am 3 

+02.0 

+01.3 

+02.2 

-01.4 

+02.1 

+01.9 + 0.2 

Australia. . 

AMi 

-16.3 

-19.6 

-19.6 

-11.2 

-26.6 

-17.3 + 1.5 


AM2 

+09.6 

+06.0 

+06.7 

+07.5 

+03.0 

+05.2 + 1.7 


Am 3 

+ 10.8 

+ 14.6 

+ 10.0 

+ 14.7 

+09.4 

+ 10.5+0.4 

Japan-Korea- 

Ami 

-08.9 

-11,5 

-08.3 

-08.5 

-06.5 

-08.9+0.5 

Manchuria 

Am 2 

+04.1 

+05.2 

+ 13.0 

+08.7 

+09.3 

+09,4+0.7 


Am 3 

+01.4 

+00.1 

+01.8 

-00.1 

+04.4 

+01.5+0.8 

Argentina 

Ami 

+35.6 

+37.9 

+39.9 

+50.7 

+34.4 

+38.3 + 1.6 


AM2 

-03.7 

+03.8 

-02.4 

-02.1 

+00.0 

-02.3+0.6 


Aua 

+ 10.0 

+07.5 

+09.9 

+04.2 

-06.3 

+05.7+3.5 

Hawaii . . . 

Am I 

+03.3 

+01.3 

-05.2 

+00.4 

-00.3 

-04.0 + 1.6 


AM2 

+06.1 

+04.5 

+ 16,0 

+01.2 

+ 15.2 

+09.2 + 2.8 


AM3 

-45.4 

-48.6 

-47.7 

-67.9 

-25.0 

-45.5 + 3.6 


parameters expressing variations in the earth’s gravitational field (i.e., 
spherical harmonic coefficients) ; p, are arbitrary polynomials of the 
Keplerian elements ; are param3ters of an atmospheric model and the 
interaction of the satellite therewith; and p^ are parameters expressing 
radiation pressure effects. 

The procedure used to compute the osculating elements M, a, e, i, u, 
and n, and the partial derivatives matrices J, Cej,, Ce(, Ced, Cep, was 
as follows. 

Preliminary orbits were determined by iterated differential correction 
fitted to the observations based on the following parameters: (a) the 
constants of integration of the orbital theory of Brouwer [1959]; (b) 
the gravitational field parameters kM and the zonal harmonics J2, Jz, 
J i] and (c) arbitrary polynomials in terms of the Keplerian elements. 
The principal purpose of this preliminary orbit determination was to 
obtain osculating elements at the instant of each observation close 
enough to the true values so that the corrections could be considered 
linear. 

The intermediate orbit elements defining the preliminary orbit were 
used to generate Fourier series to express the effects of the several 
perturbations and the partial derivatives of the osculating elements 
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with respect to the parameters of the perturbations described below. 

For Ce„ the effect of spherical harmonics of the earth’s gravitational 
field, the disturbing function developed by Kaula [1961a] was used: 




X 


( c, 
1--S, 


^ (n— »»)even 

cos {(n-2p)u + {n-2p+q)M 

nmj («— m)odd 


(fs ^ (n— m)even 

+m(0-e) sin {(n—2p)o> 

(n— m)odd 

+ (n-2p+5)7lf+m(Q-0)}l (8) 


where ko=1; k„ = 2, m^O. This disturbing function was used in the 
Lagrangian equations of motion [Brouwer and Clemence, 1961, p. 289] 
and integrated under the assumption that a, e, and i remained constant 
and that M, oj, and 12 changed secularly. The program automatically 
determined for each spherical harmonic all terms above a specified 
minimum in absolute magnitude and stored the results as subscripted 
numerical arrays to be multiphed by the sines and cosines evaluated at 
the instant of each observation. An example of one of the 210 such 
partial derivatives formed for satellite 1960t2 follows : 

56/3(731 = 1.850 cos (lo+H— 5) 

— 0.001 cos (co+Af+Q — 5) 

+5.058 cos ( — <0+12—5) 

+0.002 cos ( — <0 — A/+12 — 5) 

— 0.609 cos (— <0 — 2il/+12— 5) (9) 

Using a rejection criterion of O.ln’-^ applied to partial derivatives of 
the elements Af+<o+ 12 cos i, e^(<o+12 cos i), 12 sin i, e, i, and a between 
one and six significant periodicities were found for each term. 

The harmonics listed in Table 4 were selected as having an rms antici- 
pated effect on the satellite orbit of +20 meters or more, using the 
degree variances given by Kaula [1959], 

As expected, the partial derivatives indicated poor separation of even- 
degree harmonics of the same order m, causing principally along-track 
{perturbations of frequency m{tl — 6). However, the effect of odd-degree 
harmonics, especially third, was unexpectedly distinct, differing not only 
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in frequency but also in perturbing mainly the eccentricity (or perigee, 
height) of a nearly circular orbit. 


Table 4 . — Gravitational Coefficient Data 


[Multiply all numbers by a scaling factor of 10 


Coef- 

ficient 

Starting 

value 

Preassigned 

<r 

Preliminary 

solution 

Final 

solution 

ACoo 

0.00 

+ 10.00 

4.52 

1.23 ±3.29 

AC20 

-0.054 

+ 0.07 

0.01 

-0.06 ±0.01 

C21 

0.00 

±2.00 

-0.04 

Fixed 

*S2V 

0.00 

+ 2.00 

-0.07 

Fixed 

C22 

0.00 

+ 2.00 

2.96 

1.84 ±0.19 

<822 

0.00 

+ 2.00 

-1.71 

-1.71 ±0.28 

C30 

0.970 

+ 0.02 

1.07 

0.98 ±0.01 

C3I 

0.00 

+ 1.26 

1.89 

1.77 ±0.21 

S3I 

0.00 

+ 1.26 

0.28 

-0.11 ±0.20 

C32 

0.00 

+ 1.26 

-0.37 

0.34 ±0.26 

*S*32 

0.00 

+ 1.-26 

0.30 

0.08 ±0.35 

C33 

0.00 

+ 1.26 

-0.32 

-0.31 ±0.46 

^33 

0.00 

+ 1.26 

0.61 

0.74 ±0.46 

C.o 

0.613 

+ 0.12 

0.79 

0.55 ±0.10 

C 4 I 

0.00 

±0.63 

-0.10 

-0.21 ±0.16 

*^41 

0.00 

±0.63 

0.60 

0.46 ±0.15 

C42 

0.00 

±0.63 

0.66 

-0.03 ±0.19 

iS'42 

0.00 

±0.63 

0.25 

0.32 ±0.19 

C43 

0.00 

±0.63 

1.23 

0.50 ±0.21 

<^48 

0.00 

±0.63 

0.10 

0.16±0.19 

Cu 

0.00 

±0.63 

-0.48 

-0.24 ±0.27 

Su 

0.00 

±0.63 

1.12 

0.55 ±0.29 

C50 

0.019 

±0.02 

-0.30 

0.03+0.01 

Csi 

0.00 

±0.39 

0.35 

0.08 ±0.14 

S,1 

0.00 

±0.39 

0.70 

0.26 ±0.15 

Ceo 

-0.110 

±0.10 

0.00 

-0.10±0.02 

Cei 

0.00 

±0.28 

0.06 

0.02 ±0.08 

>^61 

0.00 

±0.28 

-0.05 

-0.18 ±0.07 

Cfi2 

0.00 

±0.28 

0.21 

0.00 ±0.06 

St2 

0.00 

+ 0.28 

0.07 

0.06 ±0.07 

Ces 

0.00 

±0.28 

1.11 

0.13 ±0.08 

*Se3 

0.00 

±0.28 

0.31 

0.21 ±0.16 

Cm 

0.00 

±0.28 

0.00 

0.13 ±0.11 

*S64 

0.00 

±0.28 

-0.28 

-0.24 ±0.10 

C7O 

0.121 

±0.02 

-0.13 

0.10 ±0.01 


Notes. — C„m, (,S„„ are coefficients of spherical harmonic terms A:iV//r(a/r)"/f„m 
such that f/f/nm’’ (iff=4ir for integration over the sphere [Kaula, 1959, equations 
(16) to (18)]. 

ACoo, ACjo are corrections to 0.3986032X10^' (1.0— O.OOIO8236P2) cgs. 
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; For tesseral harmonic coefficients, initial values of zero were assumed; 
for zonal harmonic coefficients, the values of Kozai [1962a] were used. 

For the gravitational effects of the sun and moon, a similar disturbing 
function was used [Kaula, 1962d]. All secular terms w'ere retained, plus 
periodic terms of more than 2X10-® amplitude, of which twm to nine 
were found for each orbit. 

For the radiation pressure effect of the sun, the disturbing function 
given by Kaula [1962d], was used. Because of the irregular effect of 
the earth's shadow, the perturbations were not integrated analytically, 
and a numerical harmonic analysis was applied instead. A harmonic 
analysis interval of 15 days (or a minimum period of 30 days) was found 
sufficient to reflect all variations in amplitude of more than 2X10“^ 
Partial derivatives were formed only for one parameter; the mean 
reflectivity times the cross-sectional area. 

For drag, the effect of an empirical atmospheric model was applied, 
with density in the form [Jacchia, I960]; 

'’"'’‘'(l^) cos"|| (10) 

In (10), S is the solar flux of 10.7- or 20-cm wavelength, h is the 
height above the earth’s surface, and xf/ is the angle from the center of 
the diurnal bulge, determined by 

cos { 1, 0, 0}R3(X*)Ri(«)R3(x)Rx,q/r (11) 

where X* is the sun's longitude, « is the inclination of the ecliptic, and 
X is the lag of the atmospheric bulge behind the sun. 

The atmosphere was assumed to rotate with the solid earth and to 
have the oblateness of a fluid. The customary assumption of the drag 
force being proportional to the square of the velocity was made. The 
force components that are radial, transverse, and normal to the satellite 
and its orbital plane were used in the Gaussian equations of motion 
[Brouwer and Clemence, 1961, p. 301], and numerical Fourier series were 
developed for the effects on the Keplerian elements. In generating 
these series, second-order effects on the angular elements that are de- 
pendent on the secular motions due to the oblateness were included. 
\V ith an analysis iulerv'al of 3 days, variations in amplitude as small as 
3X 10^® were obtained in M. 

For satellites 1959ai and 1959 tj the values of the parameters in (9) 
that w’ere determined by Jacchia [1960] were used. For satellite 1960t2, 
c, u, and k were set equal to zero, and po, m. H, b, n, and x vvere deter- 
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mined so as to fit the atmospheric models of Harris and Priester [1962].; 
For 1959a 1 and 1959?;, the Jacchia model absorbed most of the long-period 
drag variations but did not fit variations characterized by periods of 
less than 10 days. For 1960t2, the hlarris and Priester model did not 
reduce residuals significantly and had a negligible effect on the values 
determined for the geodetic parameters, and hence it was omitted. 

In computing the effects of arbitrary polynomials or the partial 
derivatives with respect thereto, the second-order effects of the accelera- 
tion based on the assumption of constant perigee height [equations (5) 
to (14), O’Keefe et ah, 1959; equation 2.100, Kaula, 1962c] were applied. 

In the partial derivatives J with respect to the intermediate orbit 
elements at epoch (equation 7), the effects of secular motions due to 
oblateness were included [Kaula, 1961a, equation 49]. To ensure that 
the + 20-meter specification was met, the extension of Brouwer’s theory 
to periodic terms of order J-^ by Kozai [1962b] was examined but was 
found not to be needed. 

In the final analysis of the orbit all the perturbations were added to 
the osculating elements as determined from the preliminary orbit at 
each observation. The long-period and secular perturbations that are 
due to luni-solar attraction, radiation pressure, and drag by a specified 
atmospheric model were omitted. For the orbital arc lengths of 10 to 
20 days, it was found that these effects were adequately absorbed by an 
arbitrary acceleration in the mean anomaly. Their inclusion made little 
difference in the solutions obtained for tesscral harmonics or station 
shifts — if anything, they may have distorted the results by shifting 
computed satellite directions farther from those observed. 

Data analysis . — As discussed in earlier papers [Kaula, 1961b, 1962a, 
b], difficulties are created by (1) the nonuniform distribution of observa- 
tions; (2) the similarity of effects on the observations of different gravi- 
tational coefficients and station-position errors; (3). the inadequacy of 
the atmospheric model; and (4) the prohibitive amount of computing 
time which would be required by a solution that takes into account 
serial correlation between different times. Three methods were used to 
overcome these difficulties : 

(1) Preassigning variance and covariance V for the starting values of 
parameters to which corrections z are being determined so that the 
solution becomes [Kaula, 1961a] 


g=(MnV-'M4-F“')“'MnVo-T (12) 

where W is the covariance matrix of the observations, M is the matrix 
of partial derivatives in the observation equations, and f is the vector 
of residuals. 
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: (2) Assigning higher weight to the across-track than to the along-track 

component of an observation, as described by equation 3. 

(3) Using arbitrary polynomials. 

It was found that inclusion or omission of effects which were secular 
or of periods of more than a few days had very little influence on the 
values determined for the station shifts or tesseral harmonics. The most 
troublesome inadequacy was the inability of the empirical atmospheric 
models to.explain orbital variations in the 1.0 to 0.1 cycle per day part 
of the spectrum. The principal improvement which might be made 
would be to utilize the correlation of corpuscularly caused density vari- 
ations wdth the A-p index [Jacchia, 1962]. 

It was found necessary to apply the device, specifying variance and 
covariance for the starting values of the parameters, in order to avoid 
absurdly distorted results due to the ill-conditioning caused by the non- 
uniform distribution of observ’ations and the inadequate accounting 
for drag effects. For the stations on the Am., EASI, and JKM geodetic 
systems, the 9X9 covariance matrix generated in the solution of Kaula 
[1961c] was used. For the three isolated datums, the assigned covariance 
matrices were based on assumed error ellipsoids of + 35-meter vertical 
semiaxes in all three cases, and horizontal semiaxes of ± 100 meters for 
Au., ±200 meters for Ar., and ±250 meters for H. The smaller un- 
certainty for the Australian system is based on the improvement of its 
position by adjusting deflections of the vertical [Veis, 1961]. For the 
zonal spherical harmonic coefficients of the gravitational field, the pre- 
assigned variances were based on the uncertainties given by Kozai 
[1962a], multiplied by 4. For the tesseral harmonics n, m = 2,l and 2,2, 
the preassigned variance of (2.0 X lO”®)* was based on the order of magni- 
tude of earlier determinatioi^ of J-^hy Kozai [1961], Kaula [1961b], and 
Newton [1962], except that C 21 and S 21 were held fixed. For the tesseral 
harmonic coefficients of the third and higher degrees, the preassigned cr’s 
in Table 4 were computed from the degree variances o-/ [A^] [Kaula 
1959] : 


<r'^\Cnm or Sp 


g-n [Ag] 

' in-l)Y{2n+l) 


(13) 


The observational variance employed was (0.026 sec)- time and (9.2 
sec)2 direction for 12-day arcs of 1959ai and 1959ij, (0.047 sec)^ time and 
(13.4 sec)^ direction for 20-day arcs of 1960i2; (0.146 sec)^ time and 
(43.8 sec)2 direction for 10-day arcs of 1961 5i; and (0.047 sec)^ time and 
(13.4 sec)2 direction for 18-day arcs of 1961a5i. The principal criterion 
used in determining the observational variances was the’x^ test; i.e., the 
quantity 
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s=(frw-if-z^n¥-'f/n-p) (14): 

should average 1 for several orbital arcs, where f is the vector of observa- 
tion equation residuals; W is the covariance matrix of observations; 
z is the vector of corrections to parameters; M is the matrix parameter 
coefficients in the observation equations; n is the number of observa- 
tions ; and p is the number of free parameters. In forming the covariance 
matrix W, observations in the same pass were treated as having the same 
timing error. 

The arc lengths used were chosen after some experimentation as giving 
a reasonable compromise between magnitude of residuals and number 
of observations. 

The use of arbitrary polynomials was held to a minimum; i.e., the 
only one used was a f variation in the mean anomaly. 

In determining the estimated mean value and its standard deviation 
from several orbital arcs of the same satellite, the weighting of a partic- 
ular arc was considered to be proportionate to its degrees of freedom. 
The computer program limited to fifteen the number of arcs that could 
be combined at a time. In combining the results of several sets of 
fifteen (or fewer) arcs, the weight ascribed to the mean of each set was 
considered to be the inverse of its variance, or standard deviation squared. 

In order that the final mean and standard deviation reflect as much as 
possible any systematic differences which were functions of orbital 
specifications, all sets were combined, with inverse-variance weighting, 
into four groups: 1959ai and 1959?;, twelve sets; 1960t2, two sets; 1961 5i, 
one set; and 1961a5i, one set. The final means and standard deviations 
given in Table 5 are the result of an inverse-variance weighted com- 
bination of these four group solutions. However, for most of the 
variables, the standard deviations from combining the four groups were 
smaller than the standard deviations combining all sixteen sets at once, 
primarily because the differences between the 1960t2 mean and the 
195901 and 1959?) mean were smaller than the scatter of 1959oi and 
1959?j solutions about their own mean. 

To avoid the tendency to prejudge the order of magnitude of the 
solution, which is the main defect of the preassigned-variance technique, 
some computer experimentation was tried in determining the amplitudes 
of specified periodic variations, in place of harmonic coefficients, in 
holding the reference orbit fixed, and in analyzing residuals. Applying 
these methods to one satellite at a time did not give as good results 
as the preassigned-variance method, to judge by the scatter of solutions. 
To apply them to data from more than one satellite simultaneously 
required considerable program revision which did not seem worthwhile 
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Table 5 . — Gravitational Coefficient Solutions 


[Multiply all numbers by a scaling factor of 10“*] 


Coeffi- 

cient* 

1959 ai 

1959 r; 

1960 i 2 

1961 «i 

1961 a«i 

Weighted mean 

ACoo 

4.96 

- 8.88 

- 0.75 

- 18.50 

- 9.85 

- 2.46 + 2.36 

AC20 

- 0.06 

- 0.06 

- 0.05 

- 0.29 

0.00 

- 0 . 0310.02 

C2S 

1.30 

1.36 

1.99 

1.80 

2.52 

1 . 88 + 0.29 

-Sjs 

- 1.74 

- 0.76 

- 1.63 

- 0.32 

- 0.89 

- 1 . 38 + 0.17 

Cso 

0.97 

0.96 

0.98 

1.01 

0.97 

0 . 97 + 0.01 

C31 

1.30 

1.62 

1.53 

- 0.96 

1.18 

1 . 52 + 0.03 


0.29 

0.99 

- 0.10 

- 0.34 

0.46 

0 . 14 + 0.16 

C32 

- 0.14 

- 0.13 

0.29 

2.35 

- 0.84 

- 0 . 02 + 0.26 

Si 2 

0.49 

0.29 

0.38 

- 0.16 

0.98 

0 . 4210.06 

C33 

0.36 

1.11 

0.42 

2.36 

1.70 

0 . 7010.26 

tS 33 

0.83 

1.11 

0.89 

0.43 

- 1.33 

0 . 7610.29 

C40 

0.68 

0.67 

0.61 

- 0.35 

0.62 

0 . 6710.02 

C,1 

- 0.38 

- 0.38 

- 0.33 

- 0.48 

- 1.00 

- 0 . 3310.01 

S.I 

0.43 

0.53 

0.45 

0.39 

0.45 

0 . 3710.15 

C42 

- 0.10 

- 0.10 

0.02 

0.03 

0.47 

0 . 0110.02 

S42 

0.52 

0.68 

0.36 

- 0.43 

0.06 

0 . 3510.15 

Ci , 

0.18 

0.35 

0.50 

0.44 

0.17 

0 . 1710.02 

Saz 

0.29 

0.11 

- 0.00 

0.16 

0.42 

0 . 4110.03 

Cu 

0.12 

0.01 

- 0.20 

0.20 

- 0.24 

- 0 . 0110.08 

N44 

0.11 

0.22 

0.36 

0.29 

0.32 

0 . 1810.05 

Cso 

0.02 

0.03 

0.03 

0.01 

0.02 

0.02 + 0.01 

Csi 

- 0.14 

- 0.02 

- 0.01 

- 0.63 

(•■) 

- 0 . 1310.02 

*Ssi 

- 0.06 

- 0.03 

- 0.01 

0.23 

(•>) 

- 0 . 01 10.01 

C60 

- 0.09 

- 0.08 

- 0.04 

1.10 

- 0.10 

- 0 . 0910.02 

Cel 

- 0.01 

- 0.03 

0.00 

- 0.26 

- 0.09 

- 0 . 0510.03 

<Sei 

- 0.09 

- 0.02 

- 0.07 

- 0.49 

- 0.06 

- 0 . 0610.01 

Ces 

— 0.04 

0.05 

- 0.01 

- 0.07 

0.05 

0 . 0110.01 


- 0.09 

- 0.18 

- 0.01 

- 0.07 

0.01 

- 0 . 0210.03 

Cea 

- 0.02 

- 0.10 

0.15 

- 0.02 

(•>) 

0 . 1510.01 

*Se 3 

- 0.12 

- 0.01 

- 0.08 

- 0.06 

(-) 

- 0 . 0810.01 

Ce 4 

- 0.00 

0.06 

0.06 

- 0.19 

- 0.01 

- 0 . 01 10.01 

Su 

- 0.06 

- 0.09 

- 0.42 

- 0.30 

0.03 

- 0 . 0310.07 

Cto 

0.12 

0.12 

0.07 

0.09 

0.12 

0 . 1210.01 


' Cnm and Stim are coefficients of spherical harmonic terms kM /r/a/r)" Hnm such 
that da = Air for integration over the sphere. ACoo and AC» are corrections 

to 0.3986032X10*1 (l.O-O.WlOSmP,) cgs. 

No determinations of Csi, Sn, Cst, Sei were made from l%laSi because the 
partial derivatives of the orbit with respect to these coefficients were all smaller than 
the criterion O.ln* * [Kaula, 1963a]. 


because this method has been applied extensively by Izsak [1963]. Other 
changes tried and dropped as unnecessary were deleting orbital segments 
for which observations are scanty and holding fixed the station shifts 
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obtained from the previous analysis of 1960t2 observations. The device 
of weighting observations inversely as their density with respect to the* 
phase angle (node-GST) was tried and dropped. 

Results . — The analysis described above took much time to apply to 
the large quantity of 1959o:i and 19597? data. The attempt to combine 
solutions from different sets of arcs was not made until this analysis 
had been completed. Consequently, the good agreement shown by 
Table 5 between the results from 1960i2 on the one hand and from 
1959ai and 1959 t/ on the other came as a pleasant surprise. The com- 
bination of results is not as good, of course, as is suggested by the formal 
standard deviations given in the table; in particular, the errors in dif- 
ference of position between stations in North and South America — or 
between stations in Europe, Africa, and India — which were held fixed 
with respect to each other, are probably several times as great as some 
of the stated uncertainties. The good agreement is even more marked 
for the spatial representations given in Figures 1 and 2; e.g., for the seven 



Figure 1 . — Vanguard geoid. Geoid heights, in meters, referred to an ellipsoid of 
flattening 1/298.24, determined from observations of satellites 1959ai and 


1959t;. 


most extreme maximums and minimums in the Vanguard geoid of Figure 
1, there are maximums and minimums in the Echo rocket geoid of Figure 
2 agreeing within 10° in location and within 11 meters in magnitude. 
The degree of independence in these solutions is fairly satisfying. The 
orbits differ by 0.23 in inclination, and 0.16 in eccentricity; the arc 
lengths used differed in a ratio of 5 to 3, and the observational weighting 
differed in a ratio of 3 to 2. It would be very desirable, however, to 
obtain comparable series of observations of a satellite of much higher 
inclination. 
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Figure 2 . — Echo rocket geoid. Geoid heights, in meters, referred to an ellipsoid 
of flattening 1/298.24, determined from observations of satellite 1960ij. 


The principal sources of systematic error likely to be common to 
satellites 1959ai, 1959?;, and 1960t2 seem to be (1) that the magnitudes 
of the results will be influenced by the preassigned variances and (2) 
that the relative positions of tracking stations on the same geodetic 
datum may be appreciably in error. 

For a parameter whose effects are fairly distinct in periodicities, etc., 
from those of other parameters, it is implausible that its preassigned 
variance could cause a correction that is too lai^ or of wrong sign, but 
it might cause a correction that is too small. However, the variance 
actually used in the analyses is not the estimated squared magnitude of 
the correction <ri(c), but rather Na^{c), where N is the number of orbital 
arcs in a set. Since N was always between 10 and 15, this seems to be 
no more than a mild restraint preventing occasional ill-conditioned arcs 
from obtaining absurdly large corrections beyond the range of hnearity. 

Distortion caused by the preassigned variances seems most likely to 
occur in separating gravitational coefficients whose principal effects are 
of the same period; i.e., coefficients Jnm and J^i, such that m=l and 
n—k is even. The most prominent set of such coefficients is J 22 , Ja, 
and Je 2 , all of which cause semidaily variations of argument 2(12 — fl). 
A way of removing some (but not all) of the influence of the preassigned 
covariances would be to assume that what we have determined is not 
the coefficients themselves but the amplitudes of semidaily variations 
in the orbital elements: e.g., for the cos 2(Q— fl) term in the variation 
of the inclination 


Ai= 


C22+ 


dC 




di X 

062 


dC, 


(15) 
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The semimajor axis and the eccentricity have no semidaily variation. 
If we omit the 19615i and the 1961a3i results and assume that the simi-'. 
lar 1959ai and 1959»; orbits should be combined, we have two sets of 
eight equations for three unknowns. Using values Cm = 1.315X10“®, 
Sm=- 1.473X10“', C42 = 0.101X10“', ^42 = 0.567X10“', Ce2=-0.009X 
10“', /S 62 = ~0.104X10“' for the combined 1959ai and 19597/ solution 
(corresponding to Figure 1) and using values from Table 5 for 1960i2 we 
get the computed amplitudes of periodic perturbations in columns 6 and 


Table 6 . — Semidaily Perturbations oj Satellite Orbits 


Satellite 

Ele- 

ment 

(El) 

3E1 

dC22 

3E1 

3 C 42 

3E1 

dCii 

10«X 

Comp 

AEle 

dEl 

dS‘22 

3E1 

3>S42 

3E1 

dSti 

10‘X 

Comp 

aei. 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

1959cri and 

M 

-2.92 

0.36 

5.74 

- 3.93 

2.92 

- 0.36 

- 5.74 

- 3.91 

1959>; 

i 

3.62 

- 5.89 

3.83 

5.33 

3.62 

- 5.89 

3.83 

- 9.12 

combined. 

03 

1.68 

12.21 

-24.58 

1.20 

-1.68 

-12.21 

24.58 

- 7.09 


Q 

-5.54 

4.22 

4.88 

- 7.76 

5.54 

- 4.22 

- 4.88 

-11.10 

1960.2 

M 

-6.04 

0.00 

- 0.86 

- 12.05 

6.04 

0.00 

0.86 

- 9.86 


i 

5.48 

- 5.08 

- 0.34 

10.81 

5.48 

- 5.08 

- 0.34 

- 10.76 


0) 

-2.59 

20.86 

- 2.61 

- 4.76 

2.59 

-20.86 

2.61 

-11.70 


Q 

-5.07 

- 3.26 

7.66 

-10.12 

5.07 

3.26 

- 7.66 

- 7.07 


10 of Table 6. Using these amplitudes as the observation equation 
constants and solving by the rule of minimizing 2(dAEl)^ yields 

C22 = 1.85X10“' ,§22= -1.75X10“' 

C42=0.05X10“' 542 = 0.34X10“' 

C62 = 0.10X10“' S62= -0.22X10“' 

All the coefficients are increased over the mean in Table 5 except 
S 42 , which hints of ill conditioning. However, it looks as though only 
Ce 2 and §62 might have been significantly reduced by the preassigned- 
variance method. 

The assumption that the relative positions of tracking stations on the 
same datum should be known through triangulation networks with a 
rms error of +20 meters or less was based on standard methods of esti- 
mating triangulation accuracy [Bomford, 1962, pp. 143-159], as is con- 
firmed by the misclosures of large loops of triangulation: (1) 15 meters 
in the 4000-km loop around the western Mediterranean [Whitten, 1952] ; 
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(2) less than 25 meters in the 10,000-km loop around the Caribbean 
[Fischer, 1959]; and (3) within 15 meters for the 10,000-km loop around 
the Black Sea and Caspian Sea, through Turkestan, and connecting 
in northwest India [Fischer, 1961]. 

The connections between the three northern hemisphere stations of 
the EASI system are closely associated with loops 1 and 3, and the con- 
nections between the three northern hemisphere stations of the Am. 
system are closely associated with loop 2. More in doubt are the posi- 
tions of the stations in the southern hemisphere in Peru and South Africa, 
which depend on long single arcs of triangulation. A test run was there- 
fore made on all the 1960i2 data, in which the stations Arequipa (in 
Peru) and Olifantsfontein (in South Africa) were assumed to be on sepa- 
rate datums. The results of this test corroborated the assumption as to 
triangulation accuracy; the station in Peru moved 24 meters with respect 
to those in North America, while the station in South Africa moved 14 
meters with respect to those in Eurasia. The changes in the gravitational 
coefficients were insignificant — C 22 , from 1.99 to 2.11 X10~®; *§ 22 , from 
-1.63 to -1.60X10-'; Cn, from 1.53 to 1.49X1Q-'; Cn, from -0.33 
to —0.28X10-'; etc. — and the maximum effect on any geoid height in 
Figure 2 was 4 meters. 

There still exists the possibility of errors in the local connection of 
tracking stations to the triangulation systems, a matter in which better 
standardization of procedures is needed [Kaula, 1963b]. To check this 
type of error for stations on the major datums we calculate the geometric 
geoid heights corresponding to the final positions in rectangular coor- 



Figure 5. — Combined geoid. Geoid heights, in meters, referred to an ellipsoid 
of flattening 1/298.24, determined &om observation of satellites 1959ai, 1959i|, 
196O12, 1961«i, and I 96 laSi. 
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dinates and then compare these heights with the gravitational geoid 
heights in Figure 3. To estimate the size of discrepancies to be expected, 
we have the geoid height variance of 1076 m^ from autocovariance 
analysis of gravimetry [Kaula, 1959] and a mean square height of the 
satellite geoid of 466 m^, obtained from the sum of the squares of the 
coefficients in Table 5. If the station positions and the equatorial radius 
were correct, the rms expected discrepancy between the geometric and 
gravitational geoid heights due to the inability of the satellite orbits to 
pick up the shorter-wave variations would be (1076 — 466)''^= ± 25 meters. 

The results of the comparison are shown in Table 7. Applying the 
mean correction of -|-31 meters yields a mean equatorial radius of 
6,378,196+11 meters and a rms discrepancy of +38 meters, which 
implies a rms radial position error of (382 — 252) *' 2 = +29 meters. Of 
the stations on the major datums, the 60-meter discrepancy for San 


Table 7 . — Comparison oj Geometric and Gravitational Geoid Heights^ 


Station 

Datum 

Geometric 
geoid 
height, m 

Gravita- 
tional 
geoid 
height, m 

Discrepancy for 
6,378,196 m 
radius, m 

Discrepancy for 
6,378,196 m 
radius, m 

Organ Pass 

Am. 

-4 

-10 

+ 6 

-25 

Arequipa.. 


-16 

-8 

-8 

-39 

Curagao _ _ . 


-25 

-10 

-15 

-46 

Jupiter. 


-18 

-9 

-9 

-40 

(Uifantsfontein. 

EASI 

+22 

+5 

+ 17 

-14 

San Fernando. _ 


+ 117 

+26 

+91 

-1-60 

Naini Tal . . 


-17 

-41 

-f24 

-7 

Shiraz 


+ 14 

-18 

+32 

+ 1 

Woomera . . . 

Au. 

+47 

+24 

+23 

-8 

Tokyo .. _ 

JKM 

+54 

+4 

+50 

+ 19 

Villa Dolores 

Ar. 

+ 104 

+4 

+ 100 

4-69 

Maui.. . . . 

H 

+54 

-12 

+66 

-t-35 


Note. — Geoid heights referred to ellipsoid of equatorial radius 6.378,165 meters, 
flattening 1 /298.24. 


Fernando causes suspicion of local connection error; however, there is 
also a 69-meter discrepancy for Villa Dolores, which was free to move to 
its correct position. 

The solution in Figure 3 agrees with astrogeodetic [Fischer, 1961] 
and gravimetric [Uotila, 1962] solutions, particularly in showing a more 
pronounced negative in the western Atlantic. The discrepancies which 
exist may in part be ascribed to the method of analysis of the terrestrial 
data, since the agreement is appreciably better with the combination of 
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astrogeodetic, gravimetric, and satellite zonal harmonic data of Kaula 
_ [1961c], especially for western Europe. 

’ The reference flattening of 1/298.24 is used in Figures 1, 2, and 3 
to facilitate comparison with the results of Kaula [1961c, 1963a]. The 
flattening equivalent to the solution obtained for AC 20 is 1/298.28. The 
Ji equivalent is 1082.48X10“®. 

In conclusion it can be said that better explanations are needed for 
the systematic discrepancies indicated by Tables 5 and 7. However, 
considering that the observations used herein depended on reflected 
sunlight; that they were all made more than 3 years before the minimum 
of solar activity; and that the orbital specifications are far from ideal, 
the prospects are bright for extracting more information on the gravita^ 
tional field from more recent and anticipated satellites. It will be of 
particular interest to push the analysis to a good determination of some 
sixth- or eighth-degree harmonics to see whether or not they corroborate 
other indicators of a weak upper mantle. 
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INTRODUCTION 

I N THE PAST THREE TEARS Humerous attempts have been made to 
determine the small, large-scale, longitude-dependent irregularities in 
the gravitational held of the earth from their effects on the motion of 
artificial satellites. As often happens with methods still under develop- 
ment, the first investigations on this subject were handicapped by an 
oversimplified analysis of insufficient observational material; they led to 
contradictory results and are already obsolete. It was only recently that 
determinations based on optical data [Kaula, 1963a, 6; Izsak, 1963] and 
Doppler data [Anderle and Oesterwinter, 1963; Guier, 1963] have shown 
encouraging agreement in the main features of the geopotential as inde- 
pendently obtained by these authors. At the same time their results com- 
pare favorably with those derived from an analysis of surface-gravity data 
[Kaula, 19615; Uotila, 1962]. Despite the reported progress, much work 
must be done before we can reach definite conclusions about the longi- 
tudinal dependence of the earth’s gravitational field. The purpose of 
this paper is to present new evidence concerning the subject. No at- 
tempt will be made to evaluate critically the methods and results of the 
authors mentioned above, for this would seem premature at this time. 

The basic philosophy adopted in the present investigation is that a 
statistical problem of such complexity should be considered from an 
experimental rather than from a theoretical point of view. Instead of 
trying to anticipate what might happ>en in the process of a very intricate 
least-squares fitting, I think it preferable to carry out many different 
solutions under conditions as different as possible and then to choose the 
best solution by comparing residuals, standard deviations, correlation 
coefficients, etc. Questions regarding the convergence of the geopoten- 
tial’s development into an infinite series of spherical harmonics seem 
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irrelevant to me. A least-squares fitting of a finite number of terms can 
merely give an approximation in the mean, which, in practical applica- 
tions, is what we really are concerned with. 

Our main interest is in the external gravitational field of the earth. 
However, since the interpretation of satellite observations is affected by 
the errors in the assumed coordinates of the Baker-Nunn camera stations, 
it is imperative to improve the latter. Experimental evidence shows that 
the correlations among harmonic coefficients and station coordinates are 
generally weak, and those among corrections to the coordinates of 
different stations are entirely negligible. Since the original connection 
of the stations to major geodetic systems proved to be unreliable in 
several instances, I decided to treat the 12 tracking stations individually. 

COMPUTATIONAL PROCESS 

The analysis of the observations begins with a painstaking computa- 
tion of satellite orbits. 

DOI 3, our differential orbit improvement program [Gaposchkin, 1964], 
accounts analytically for the first-order short-period perturbations result- 
ing from the oblateness and for the long-period perturbations caused by 
the equatorial unsymmetry of the earth, as well as for lunar perturba- 
tions with biweekly periods. Secular perturbations of gravitational 
origin, the effect of atmospheric drag and solar radiation pressure, how- 
ever, are as a rule determined empirically in the process of least-squares 
fitting. The length of an orbital arc represented by a single set of mean 
orbital elements varies between one and four weeks, depending on the 
particular object ; in such a time interval there are 60 to 500 observations. 
On completing the computations, DOI 3 produces for each observation 
a binary card that contains among other useful information the sidereal 
time at Greenwich di, the along-track and across-track residuals dU t 
and dll i of the observation, as well as the instantaneous values AT,, a^, Cj, 
/ (, coj, and n, of the mean orbital elements. The binary cards pertain- 
ing to several satellites in several time intervals are collected, and their 
content is put on a magnetic tape, which will serve as input in the 
subsequent analysis. 

The residuals df q and dW’i comprise, besides components not suffi- 
ciently accountable, the positional and timing errors of the Baker-Nunn 
observations, the effect of the tesseral harmonics in the geopotential, 
and the effect of small errors in the adopted station coordinates. By 
analyzing many thousand observations one can arrive at reasonable 
least-squares estimates of the tesseral harmonics’ numerical coefficients 
and of the corrections to the station coordinates. This is the purpose 
of the tesseral harmonics and station coordinates program, a brief 
description of which follows: 
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Let the geocentric position vector of a satellite as computed by DOI 3 
.be denoted by r^, and a station vector as used in the computations be 
Rc : the corresponding topocentric position vector of the satellite is Pc = 
ft ~ Rc- If the correction to be applied to the station vector is dR = R — R^, 
and the perturbations caused by the hitherto neglected tesseral harmonics 
are given by dr = r— r^, the corresponding change in the topocentric posi- 
tion vector will be dp = dr — dR. Having only angular residuals at our 
disposal, we disregard the line-of-sight component of this vector by 
projecting it onto the directions of the along-track and across-track unit 


vectors Cv and Cir. For the angular corrections in these directions 
have then 

we 

du = du'‘—du‘ 

where 

dw = diD’‘—dw‘ 

(1) 

du‘ = p~^evdR 

and 

dwj*=p-ieir-dR 

(2) 

du’‘ = p~^ev‘dT 

dw"=p~^ew’dr 

(3) 


The positions of the Baker-Nunn stations [Veis, 1961, 1963a] being 
given in an earth-fixed rectangular coordinate system X,Y, Z, we put 


dR — dX-\-er dY-\-ez dZ (4) 

Here the unit vector ex lies in the intersection of the equator with the 
meridian of Greenwich; Cz points toward the north pole; and er is 
perpendicular to both of them. 

The explicit representation of the vector dr is much more involved, 
of course. It requires the elaboration of a perturbation theory for the 
tesseral harmonics in the geopotential, a topic that will be touched upon 
below. At present it suffices to say that any perturbations of the orbital 
elements result in a vectorial displacement, 

and that the perturbations in question are linear in the numerical co- 
efficients of the tesseral harmonics. Full details concerning equations 1 
throiffeh 5 were given in an earlier paper [Izsak, 1962]. 

In this way we obtain the observation equations 

duj = dUi"— du‘ = dUi 

and 

dwi — diVi‘~diL'i=dWi ( 6 ) 

To be minimized is 

[{dUi — dut)^+{dWi — dWiy] 
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where the summation is extended over the available observations of a 
number of satellites at all participating stations. An earlier version of_ ■ 
the program had the facility of applying different weights to the along- 
track and across-track residuals, but the advantages of such an approach 
seemed rather dubious in practice. 

Concerning the perturbation theory for tesseral (and sectorial) har- 
monics I confine myself to a few remarks. The disturbing function, 
which determines Lagrange’s equations for the variation of orbital 
elements and thus the perturbations of the latter, is written in the form 

= • %i^nm COS mX+*S„„ sin mX)P„m(sin vp)| 

the adopted definition of the Legendre associated functions being 
Pnm{x) = d”Pn{x)/dx^ 

To utilize Lagrange’s equation, we must develop the disturbing function 
in terms of the orbital elements. This is conveniently done in two steps. 

We note first that any rotation of the X, Y, Z coordinate system 
induces a linear transformation of the nth-degree spherical harmonics 
[Wigner, 1959], while the elements of the transformation matrix are 
easily calculated as functions of the Eulerian angles of rotation by means 
of the corre.sponding Cayley-Klein parameters. Let the X' and Z' axis 
of a new coordinate system point toward the perigee and pole of the 
satellite’s orbit, so that the Eulerian angles coincide with Q-d, I, and u. 
In this system the polar coordinates of the satellite are X' = r and v' = 0; 
therefore the spherical surface harmonics degenerate into trigonometric 
functions. Implementing these considerations we arrive at the expansion 


P,„.(sin ip) exp = exp {i[{n-2j){v+w)+m{a- 6)]} 

J-O 

where 

"2"i!(n-i)! ^ 


'2n-2jl 

1 

--I 

1 

Vn — m — k^ 


2n— »> V /n\ 

y <T \J) 


7 = cos(//2), <r = sin(//2), v = m-n-\-2j-\-2k, and the range of the k 
summation is 

max {0, n-m-2j| </c<min {2n—2j, n-m] 

The polynomials (7) have the symmetry 

K„A(7lcr) = (-l)"-’"K„„."-V|7) 
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In the second step, the functions 
* exp [f(n— 2j)i’] 

are to be developed into a. Laurent series of exp {iM), a process that 
invokes the well-known Hansen coefficients from celestial mechanics. 
This and related subjects were dealt with recently [Izsak et al., 1964]. 

For the actual computation of the effect of tesseral harmonics on the 
orbital elements, we used a program that William M. Kaula of the 
Goddard Space Flight Center kindly made available to us; its mathe- 
matical formulation has already been described by Kaula [1961a]. Once 
given the orbital elements a, e, and I, this program yields for the several 
harmonics apart from the at-this-stage-unknown numerical coefficients 
Cnm and the perturbations dM , • • •, dQ in the form of trigonometric 
poljmomials. Three terms are used if n is even, and four if n is odd. 
Thus, for instance, the effect of the second-, third-, and fourth-degree 
harmonics upon each orbital element is represented by a total of 54 
trigonometric terms. The complete expressions for the perturbations are 


n.m 

: ( 8 ) 
dn = X(C„„Pu'‘”+Snmqa’‘”) 

n.m 


the symbols p.v"'", •••, pa"", qo”” standing for the trigonometric 

polynomials given by Kaula’s program. As examples, pertaining to 
the mean orbital elements a = 8306.8 km, e = 0.16446, and 7 = 32.883° of 
the satellite 1959ai (Vanguard 2), we choose 

Pa/“ = 4.008 sin [2w+2M-f-2(n—0)] 

-2.979 sin [2(12-0)] 

-0.025 sin [-2w-2M+2(n-e)] 

and 


9a®*=— 2.255 cos [— u-l-(0— 0)] 

-f5.481 cos [w-f (12— 0)] 

— 1.113 cos [— w— il7-f (12— 0)] 

0.135 cos [w-[-A7-[- (12 — 0)1 

As a matter of fact, Kaula’s program was adapted to fully normalized 
spherical harmonics, the integral of whose square over the unit sphere is 
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4t. The relation between the coefficients Cnm, S„m of the ordinary, and 
C„m, S„m of the fully normalized, spherical harmonics is 


where 


Cnm — ^ nmC nm ^nm ^nm^n 


Nno={2n+l) 
nn 


1/2 


2(2n + l)(n — m)! 
(n+m)! 


1/2 




Let us now return to the observation equations 6. They are set up 
as follows. Comparison of (3) and (4) gives for the station that secured 
the observation 


du’ = dX-\-v7 dY +u^ dZ 

dw’ = w^ dX-\-w''^ dY-Yw^ dZ 


where 

M-*’ = p-i(ec7-ex), • • •, = 

As to the harmonics, from (2) and (5) we have 

du'‘ = u^ da+u‘ de+u^ dl+u“ do>+u° dQ 

dw" = 0.dM+w‘‘ da+w‘ de+w^ dl+w‘^ dw+w° dO 

with the abbreviations 

= p~^{eu‘Sr/dM),’ • •, vP = p~^{ew‘dr/dQ) 

Then, using (8) and collecting terms that belong to the same harmonic, 
we get 

du'‘— '^{uc'™Cnm-\-Us’‘"‘S„m) 

n,m 

dw''=XiWc’'’"Cnm + Ws’‘’^Snm) 

n,m 

where 

= = h wV"” 

= 

Once the observation equations are set up, they are solved by the 
computer using standard least-squares techniques. In its present version 
our program can simultaneously handle up to 150 orbital arcs of 10 
different satellites; this corresponds to about 20,000 individual observa- 
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tions. The selection of the unknowns to be solved for is arbitrary. 
Their number, however, is limited to 44. Of these unknowns not more 
Tthan 38 may be tesseral-harmonics coelEcients, and not more than 36 
may be corrections to station coordinates. The printout of the program 
has been devised to give detailed information concerning the intricate 
least-squares process under discussion. Its main features are : 

(1) A list of the original (DOI 3) ; the stations — ^improved (if any) ; 
the harmonics — ^improved (if any) ; and the improved residuals 


and 


dUt dWt 

SUt’‘=dU,-dUi'‘ 8Wi’‘=dWt-dWi'‘ 


bUi = dUi — dUi 


&Wi = dWt—dwt 


(2) The mean values 

! N -i 1/2 

|(2A)-';g[(db^.)=-b(dTF,)==]| 

and 

|(2A)-‘^[(6C7d^+(5IF<)2]|''‘ 

of the original and of the improved residuals per orbital arc, as well as 
for all the observations. We call their ratios the respective improve- 
ment factors of a solution. 

(3) The least-squares estimates of the fully normalized tesseral and 
sectorial harmonics coefficients and of the corrections to station coor- 
dinates together with their standard deviations. 

(4) The matrix of the correlation coefficients, that is, the normalized 
inverse matrix of the normal equations. 


DATA ANALYSIS 

As indicated by the foregoing description of the computational proc- 
ess, the mathematics involved is rather straightforward; the difficulty 
of the problem lies in the analysis of real data. First of all, the informa- 
tion contained in the original residuals is contaminated by unaccountable 
effects with a total that is significantly larger than the well-established 
accuracy of the observations : 2 seconds in the across-track and 3 seconds 
in the along-track direction. Second, while the perturbations caused 
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by the presence of the various spherical harmonics are continuous func- 
tions of time, optical observations necessarily result in discrete data 
of a distribution that, for a single orbital arc, is often anything but 
uniform. This situation, though especially troublesome in case of per- 
turbations with high frequencies, is comparatively harmless so far as 
station coordinates are concerned. Third, the effect of the different 
harmonics is not very much different in size and, for some terms, not 
even in shape. Therefore the perturbations caused by the various 
harmonics are hard to separate, a fact borne out by eventually consider- 
able correlations among the computed values of the unknowns. About 
the only thing one can do to minimize the influence of these disadvanta- 
geous circumstances is to use a huge number of observations of several 
satellites with a variety of orbital elements, particularly with different 
inclinations. The self-explanatory Table 1 constitutes a breakdown of 
the 15,191 observations on which this analysis is based. 

Table 1 . — Observational Material 


Satellite 

7, deg 

e 

a, km 

Arcs 

Number of 
observations 

1959ai, Vanguard 2 - - 

32.9 

0.165 

8301 

15 

2091 

1959a2, Rocket of 

32.9 

0.183 

8492 

2 

229 

1959t;, Vanguard 3 . 

33.4 

0.189 

8496 

12 

1824 

1961 5i, Explorer 9_ 

38.9 

0.115 

7985 

16 

2446 

1962aei, Telstar 1... 

44.8 

0.242 

9672 

3 

269 

1960i 2, Echo 1, rocket 

47.2 

0.012 

7972 

34 

4324 

1962/3^1, Anna IB .. 

50.1 

0.007 

7508 

3 

296 

1961oi, Transit 4A 

66.8 

0.008 

7318 

13 

971 

1961o 2, Injun 1 — Greb 3 

66.8 

0.008 

7319 

7 

494 

1961a5i, Midas 4 

95.9 

0.011 

10005 

6 

2247 


In the present computational approach, we need the coordinates of 
the Baker-Nunn camera stations in some uniform rectangular system. 
The method of reducing the fundamental geodetic coordinates of the 
stations, based on field surveys and referred to five different geodetic 
datums to such a system, has been given by Veis [1961]. The actual 
values used in our computations, however, come from later, yet unpub- 
lished, work of his. They were determined on the basis of improved 
geodetic information available June 15, 1962, and (instead of the inter- 
national ellipsoid) reduced to the ellipsoid (recommended by Kaula) with 
the characteristics a£ = 6378.165 km and /= 1/298.3. For the sake of 
easy reference, we list in Table 2 the initial coordinates of our tracking 
stations as derived by Veis. 
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Table 2 . — Observing Stations 



X, km 

y, km 

Z, km 

Datum 

Number of 
observations 

1. Organ Pass. . . 

-1535.702 

-5167.026 

3401.108 

Am. 

1826 

2. Olifantefontein. - 

5056.123 

2716.523 

-2775.799 

Eu. 

1742 

3. Woomera 

-3983.602 

3743.226 

-3275.656 

Aus. 

2023 

4. San Fernando 

5105.623 

-555.194 

3769.670 

Eu. 

1315 

5. Tokyo 

-3946.522 

3366.453 

3698.855 

Jap. 

965 

6. Xaini Tal 

1018.135 

5471.207 

3109.519 

Eu. 

161 

7. Arequipa _ 

1942.762 

-5804.082 

-1796.838 

Am. 

931 

8. Shiraz 

3376.872 

4404.022 

3136.250 

Eu. 

1123 

9. Curasao.. . . 

2251.841 

-5816.928 

1327.236 

Am. 

831 

10. Jupiter 

976.319 

-5601.410 

2880.311 

Am. 

1458 

11. Villa Dolores 

2280.645 

-4914.512 

-3355.441 

Arg. 

1276 

12. Maui . 

-5466.118 

-2404.068 

2242.437 

Am. 

1540 


The satellite observations set forth in Table 1 have been used recently 
in a great variety of computer runs. In some cases the required machine 
time was well over an hour. There is little point in giving a full account 
of the numerical results obtained. Instead, we select samples of them 
in the form of tables that we think are pertinent to the problems under 
discussion. All (fully normalized) harmonics coefficients will be expressed 
in units of 10“®, and the station corrections will be given in meters. 
Throughout the computations we used the value = = 3.986032X10“ 

cm’ sec”’ recommended by Kaula. 

Table 3 shows how futile it is to use a single satellite even for the 
determination of one pair of coefficients. The results per satellite are 
incompatible, probably because the effect of the several harmonics 


Table 3. — Second-Degree Solutions per Satellite 

(Scaled by 10') 


Satellite 

(?22 

^22 

Residuals, sec 

Original 

Improved 

1959ai 

0.97 + .08 

-2.17 ±.08 

14.27" 

12.95' 

1959 i 7 

1.22 ±.10 

-2.32 ±.10 

13.53 

12.49 

19618, 

1.50 + .06 

-1.64 ±.06 

15.56 

13.63 

1960,2 

1.51 ±.03 

-1.37 ±.03 

16.02 

12.67 

1961o, 

0.81 ±.04 

-0.52 + .04 

23.61 

20.43 

196102 

0.85 ± .06 

-0.26 ±.06 

24.36 

21.50 

1961aS, 

2.29 ±.05 

-0.52 ±.05 

9.84 

7.99 

All ten 

1.08 ±.02 

-0.89 ±.02 

15.67' 

13.86' 
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averages out differently according to the inclination of the objects. The 
considerable differences in the mean angular residuals result primarily' 
from the variety in the mean topocentric distances of the objects. 

Table 4 illustrates the interaction of harmonic coefficients and station 
corrections in a least-squares adjustment. In the first run only the 
harmonics up to the third degree were considered unknown; in the second 
only the stations were varied. In the third and definitive run these 
two problems were solved simultaneously. Clearly, the presence of tes- 
seral harmonics influences the computed station coordinates much more 
than small errors in the latter affect the computed harmonic coefficients. 
The most conspicuous correlation coefficients among harmonics and 


Table 4. — Third-Degree Solutions for Harmonics (^Scaled by 10^) and 
Corrections to Station Coordinates {in meters') 


C 22 ^22 

1.03 - 

1.14 




1.02 + .02 -1.13 + .02 


Csi Szi 

1.02 - 

0.32 




1.35 + .02 -0.18 + .02 


Cs2 *Ss2 

0.30 - 

0.45 




0.35 + .03 -0.73 + .03 


C 33 >^33 

-0.31 

0.99 




-0.37 + .04 1.01 +.04 


dX' dYi dZi 



-10 

19 

51 

-40+ 5 

34+5 

-20 + 

5 

dXi dYi dZi 



31 

64 

-10 

11+ 5 

-26+ 5 

-32 + 

6 

dX) dY> dZi 



2 

-20 

13 

-67+ 4 

-16+ 4 

10 + 

4 

dXi dY, dZt 



42 

-72 

10 

-24 + 11 

-34+ 8 

68 + 

8 

dX t dY s dZ ^ 



-47 

-89 

-41 

-101+ 7 

-80+ 7 

12 + 

8 

dX, dY, dZ, 



49 

-64 

84 

44 + 16 

-76 + 26 

75 + 20 

dX^ dYj dZ^ 



-59 

-44 

-65 

42+7 

-1+ 8 

-103 + 

8 

dX, dY, dZ, 



77 

-39 

115 

33+7 

-40+ 7 

59 + 

7 

dX, dY, dZ, 



-102 

-55 

-39 

-10+ 7 

6+ 9 

-64 + 

8 

dXiQ dY 10 dZio 



-21 

5 

0 

-41+ 5 

5+ 6 

-36 + 

6 

dXii dY ii dZ\\ 



-23 

-43 

9 

-19+ 6 

-53+ 7 

-21 + 

7 

dX\2 dY 12 dZ\2 



53 

-152 

-138 

37+6 

-152+ 5 

-163 + 

6 

Improved 










residuals 

13.03" 



14.85" 



12.11" 




stations are (S32, dF 4 )=— 0.25 and (<832, dXio) = 0.25 and those among 
the coordinates of different stations turn out to be (dXi, dY ^ = —0.07 
and (dXio, dFi 2 ) = 0.07. These data refer, of course, to the third run, 
whose results concerning the station coordinates we regard as the best 
the present version of our computer program can yield in conjunction 
with the observational material used in this work. A generalization of 
the program now under development will be indicated below. 

At this point it is interesting to compare the results of the third run 
with those obtained by Kaula [19636] and Veis [1963c], who were both 
using dynamical methods. In addition, we consider the purely geo- 
metrical method of simultaneous observations, initiated and pursued at 
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the Smithsonian Astrophysical Observatory by Veis [19636], This tech- 
nique permits the determination of two comjxments and of relative 
displacement of the participating stations, perpendicular to the direction 
of one station as seen from the other. Simultaneous observations have 
already been successful in the instance of four pairs of neighboring sta- 
tions, thus providing a valuable check on results arrived at by dynamical 



Direction 7 to 11 
( 1 826 km ) 


Direction 9 to 10 
(2021 km) 


Scale 



Figure 1 . — Comparison of satellite results for the relative direction of Baker-Nunn 
camera stations 1° [Veis, 1963*]; 2° [Kaula, 1963c]; 3° [Veis, 1963c]; 4“, Izsak, 
Table 4. 
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methods. Figure 1 compares these results. The origin of the individual 
diagrams corresponds to the assumed direction of the respective stations! 
Small circles represent relative displacements according to 1° [Veis, 
1963h], where error ellipses are also given; 2° [Kaula, 19636]; 3° [Veis, 
1963c]; and 4°, our Table 4. 

The average agreement of the geometric results, 1°, seems to be best 
with the dynamical results, 4°. In case of the direction of Arequipa- 
Villa Dolores (7-11), the discrepancy between the determinations 2° 
and 1° is so striking as to call for a careful investigation. 

During the experimental phase of this work the harmonic coefficients 
were computed in a variety of combinations. The purpose of such 
numerical experiments was merely to establish the capabilities of our pro- 
gram. We soon concluded that although the individual coefficients C„m 
and Snm may vary significantly from run to run, the combined contribu- 
tion of the several harmonics to the gravitational field of the earth is 


Table 5. — Fourth-, Fijth-, and Sixth-Degree Solutions for Harmonics 

(Scaled by 10') 


Cil 

1.06 - 

0.73 

1.08 -0.70 1.17 ±.02 

-0.95 ±.03 

Cz\ >?3l 

0.88 - 

0.26 

0.84 -0.21 0.81 ±.02 

-0.25 ±.02 

C32 jSsj 

0.22 - 

0.36 

0.30 -0.31 0.24 ±.03 

-0.25 ±.03 

C33 Sii 

-0.30 

0.79 

-0.53 0.87 -0.50 ±.04 

0.93 ±.04 

C4I *Si 4 l 

-0.18 - 

0.24 

-0.19 -0.20 -0.18 ±.01 

-0.25 ±.01 

C42 *^42 

-0.16 

0.52 

-0.19 0.46 -0.11 ±.02 

0.23 ±.02 

C43 'S43 

0.34 - 

0.08 

0.32 -0.08 0.28 ±.02 

-0.08 ±.02 

C44 *S *44 

0.06 

0.59 

0.20 0.57 -0.08 ±.05 

0.29 ±.06 

Csi Sh\ 



-0.12 0.20 -0.09 ±.01 

0.19 ±.01 

Cbi Sb 2 



0.30 -0.49 0.31 ±.03 

-0.50 ±.03 

Cbt Sbz 



-0.91 0.17 -0.72 ±.05 

0.11 ±.05 

Cbi Sbi 



-0.21 0.53 -0.18 ±.06 

0.51 ±.06 

Cbb Sbb 



0.17 -0.35 0.18±.10 

-0.42 ±.10 

Cbi *Sei 




-0.01 ±.01 

0.13 ±.01 

Cb 2 




0.16 ±.02 

-0.37 ±.02 

CbS Sb 3 




0.14 ±.02 

-0.17 ±.02 

Co 4 iS «4 




-0.20 ±.04 

-0.41 ±.05 

Cbb Sbb 




-0.40 ±.04 

-0.28 ±.04 

Cbb Sbb 




-0.53 ±.08 

-0.41 ±.08 

Improved residuals 

12.50" 


12.27’ 

12.10’ 


rather well determined by sufficiently rich observational material. Table 
5 is a compilation of recently obtained fourth-, fifth-, and sixth-degree 
solutions for harmonics, the last of which we regard as the definitive 
one. There seems little justification at present for going beyond sixth- 
degree terms. The standard deviations of the computed coefficients, 
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• whatever their significance, increase sj'stematically with the index m, so 
that they are largest for the sectorial terms. We enumerate here only 
the most notable correlation coefficients; 


{Ci,, C42) = 0.72 
(S22, S42> = 0.78 
<C«, Co3) = 0.3G 
(C4t, C6 i)=-0.38 

<C42, C62) = 0.61 
(C 44 , C64)~0.65 


(C22, C62) — 0.55 
(S22, S62> = 0.56 

{533, 5=3) =0.36 
{S41, 5ei)= — 0.36 
< 542 , 562 > = 0.63 
< 544 , 564 > = 0.67 


To facilitate the visualization of numerical results, level curves of 
geoid heights pertaining to our sixth-degree solution are plotted in Figure 
2. The reference ellipsoid of oblateness 1/298.28 and the purely latitu- 
dinal deviations from it are defined bj" the zonal harmonics coefficients of 
Kozai [1962]. The most interesting deviations of this figure from that 
in Izsak [1963] and Kaula [19635] occur over the Atlantic and Xorth 
Africa; their reality, however, cannot be ascertained without further 
evidence. 

Extensions of the present computer program are feasible in a number 
of ways. First of all, one could increase the number of unknowns to 52, 
so that a simultaneous solution for the 12 stations of our network and 
second-, third-, and fourth-degree harmonics would become possible. 
Or, wishing to bring about a better separation of the geometric and 
dynamical part of the problem, one might be inclined to iterate previously 
obtained solutions. For instance, one could start with the results of the 
third run in Table 4, take its improved residuals and use these in a solu- 
tion for higher-degree harmonics. The effect of the harmonics thus 
derived, in turn, vould be subtracted from the original residuals when 
starting a new solution for the stations. 


Appendix . — Normalization Coefficients 


nm 

nm 

nm 

A nm 

22 

0.64550 

51 

0.85635 



52 

0.16183 

31 

1.0801 

53 

0.033034 

32 

0.34157 

54 

0.0077863 

33 

0.13944 

55 

0.0024622 

41 

0.94868 

61 

0.78680 

42 

0.22361 

62 

0.12440 

43 

0.059761 

63 

0.020734 

44 

0.021129 

64 

0.0037855 



65 

0.00080707 



66 

0.00023298 
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Figure Z— Level curves of geoid heights at 10-m intervals. The small triangles indicate the positions of the 12 Baker-Nunn camera 
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Abstract . — From Baker-Nunn observations of nine satellites, whose inclinations 
cover a region betw’een 28° and 95°, the following values were derived for the zonal 
harmonics coefficients of the earth’s gravitational field: 


Jt =1082.645X10-', 
±6 

Jt = -1.649X10-', 
±16 

J, = 0.646X10-', 

+ 30 

J»= -0.270X10-', 
±50 

J,„= -0.054X10-', 
±50 

Ja= -0.357X10-', 
±47 

7u= 0.179X10-', 

±63 


J, =-2.546X10-', 
±20 

Ji =-0.210X10-', 
±25 

■J, =-0.333X10-', 
±39 

J, =-0.053X10-', 
±60 

Jn= 0.302X10-', 
±35 

Ji3 = -0.114X10-', 
±84 


1. INTRODUCTION 

I N A PREVIOUS PAPER (Kozai, 1963) I derived a set of values for the 
coefficients of zonal spherical harmonics in the earth’s grativational 
potential from the available observations of artificial satellites. How- 
ever, at that time I did not give much weight to observations of high- 
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inclination satellites simply because accurate observations for such 
satellites were not available. 

We now have precisely reduced Baker-Nunn observations for some of 
the high-inclination satellites, and I have found that secular motions of 
ascending nodes of these satellites cannot be accurately expressed by my 
previous values of zonal harmonics. Therefore, I had to improve my 
previous values by adding observations of the high-inclination satellites 
and higher-order harmonics to the expression of the earth’s potential. 

In this paper I have tried to eliminate any accidental errors in observa- 
tional data, by using many more observations of a given satellite than in 
my previous paper. I have used fourteen sets of observations for 1959 al 
and ten sets for 1959??, in contrast to the single set of data used for each 
satellite previously. Consequently, I believe that the data reported here 
are more reliable than those in the previous paper even for low-inclination 
satellites. Although we still lack sufficient observations for satellites 
with inclinations of between fifty and eighty degrees, this gap in the data 
will probably be filled in the near future. 

2. METHOD OF REDUCTION 

The observations used in this determination were made by Baker-Nunn 
cameras, and the first steps in the reductions were made by Phyllis Stern 
by the Differential Orbit Improvement program, in which first-order 
short-periodic perturbations due to the oblateness of the earth are taken 
out. The mean orbital elements of each satellite for every two days or 
four days were obtained from observations covering four or eight days. 
Luni-solar periodic itnd solar radiation perturbations in the orbital 
elements were then computed and subtracted from the mean orbital 
elements. 

To derive secular motions of the ascending node and the perigee and 
amplitudes of long-periodic terms from these orbital elements, I use data 
covering about one period of revolution of argument of perigee, that is, 
about 80 days for Vanguard satellites, for example. 

Secular accelerations in the mean anomaly or the mean longitude, and 
secular decreases in the semimajor axis due to air-drag, are then evaluated 
roughly; they can be used to compute theoretically secular variation in 
the longitude of the ascending node, the argument of perigee; and the 
eccentricity due to the air drag with sufficient accuracy, by assuming the 
rate of secular decrease of the perigee height. The computed secular 
variations in the three orbital elements are subtracted from the mean 
elements. 

After the corrections with long-periodic perturbations due to even 
zonal harmonic terms are made, the argument of perigee co, the longitude 
of the ascending node, fi, the inclination i, and the eccentricity e are 
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expressed by the following simple forms: 

O) cos O), 

n= flod” cos CO, 

f = fo+-4i sin CO, (1) 

e = eo-\-At sin co. 

By the method of least squares we can determine the constants appear- 
ing in the formulas (1) from a set of the corrected orbital elements. How- 
ever, when the eccentricity is very small, say less than 0.02, the corrected 
eccentricity and the ailment of perigee are more accurately expressed 
by the following formulas ; 

€ sin co = eo (1— a) sin (cood-cbO+^e, 

(2) 

ecosco = eo (1+a) cos (cood-cot), 

where a, w'hich is due to even-order harmonics, can be computed with 
approximate values of J„ as 

a = sin“ (14-15 sin^ 2)d-5 (6-7 sin^ i) 

— 10.9375 J6(16-48 sin“ i‘d-33 sin^ f)/a^}/{16a^ ^2(4-5 sin^ i)}. (3) 

By using the formulas (2) w’e can determine Co sin &)o, co cos oio, and 
a correction to an assumed value of cb from observations by the method of 
least squares. 

The relation between the anomalistic mean motion n and our semimajor 
axis a is given as 

nV = Gd/|ld-|^(l-e2)^(l-3 cos^ f)|, (4) 

where 


GM = 3.986032 X 10“ cmVsec=, 

p = a(l-e2). W 

Expressing the mean motion in revolutions per day and the semimajor axis 
in earth’s equatorial radii, we can use the following number for GM : 

V’^ = 17.043570, (6) 

where I adopt the following value of the equatorial radius: 

a, = 6378.165 km. (7) 
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• ^ The earth’s gravitational potential is expressed with Legendre poly- 
nomials as 



f >-=2 


•/n(a.A)" Pnisin /S) 


} 


( 8 ) 


The secular motions of the node and the perigee and the amphtudes 
of long-periodic terms with argument to derived from observations are 
compared with those computed from my previous value of (Kozai, 
1963), 


72 = 1082.48X10-', 

7, = -2.562X10-', 

74= -1.84X10-', 

7s =-0.064X10-', 

7s= 0.39X10-', 

7t= -0.470X10-', 

7g= -0.02X10-', 

7,= 0.117X10-'. 


(9) 


Of course we must include luni-solar secular terms and a term, 
which can be computed with an approximate value of Jt to compute 
secular motions. Therefore, each secular motion and amplitude pro- 
vides us with (0-C), which will make it possible to improve values of J,. 


3. DATA 

(a) 1959 Alpha 1— Table 1 lists fourteen sets of data for this satellite, 
and table 2 gives (O-C)’s referred to my previous values for J^. 

The standard deviations for the daily secular motions u and Q given 
in table 1 are determined from observations; those in table 2 are com- 
puted by adding uncertainties which come from those in Co and io. 
Weighted mean values for the fourteen sets are given at the bottom of the 
table. As can be seen, the scattering of (O-C)’s is much lai^er than that 
expected from the standard deviations assigned to the observed values. 
However, the standard deviations assigned to the mean values in table 2 
should be more reliable, and will be used in the determinations of J,. 

(b) 1959 Eta — Ten sets of data are given in tables 3 and 4 for this 
Vanguard satellite. However, its orbital elements are not essentially 
different from those of 1959 al and the mean values of (0-C) in table 4 
are almost identical with those in table 2, as expected. For the two Van- 
guard satellites (0-C) in tl and A„ are significantly large. 

(c) 1960 Iota 2 — Since the eccentricity is veiy small for this rocket of 
Echo I, the fonnulas (2) are used in the reduction. Since «-f- Q are very 
small for this sateUite, it is necessary to take special care to compute terms 
with arguments 2({o-l-Q— fl©) and 2(w-|-i2— flc) in the luni-solar per- 
turbations. 

Five sets of data are given in tables 5 and 6. For this satellite the 
scattering of (0-C) for secular motions is very large. The large scatter- 
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Table 2. — (0-C) Referred to Kozai’s Previous Constants for 1959 Alpha I . 



wXl0‘ 

QX10« 

A,X10« 

A.X105 


K10‘ 

AnXlO* 

1 

19° 

+ 17° 

-31° 

+ 18° 

12+7 

13° 

+ 17° 

90° 

+ 40° 

-9° 

+ 5° 

2 

0 

+ 

8 

49 

+ 23 

17 + 5 

-26 

+ 18 

2 

+ 23 

19 

+ 7 

3 

3 

+ 

8 

-23 

+ 23 

18 + 6 

-13 

+ 33 

33 

+ 36 

F- 

CO 

1 

+ 9 

4 

-15 

+ 11 

-22 

+ 21 

2+6 

41 

+ 16 

102 

+ 29 

-23 

+ 4 

5 --- 

3 

+ 

6 

-42 

+ 13 

7+4 

-3 

+ 7 

41 

+ 17 

-1 

+ 3 

6 

-3 

+ 

4 

-33 

+ 10 

0 + 4 

-31 

+ 14 

12 

+ 12 

10 

+ 3 

7 

6 

+ 

3 

3 

+ 14 

7+5 

44 

+ 13 

12 

+ 10 

11 

+ 6 

8 

-5 

+ 

7 

-27 

+ 17 

7 + 3 

24 

+ 21 

98 

+ 18 

-20 

+ 5 

9 

1 

+ 

5 

-24 

+ 16 

-4+5 

6 

+ 30 

37 

+ 13 

-14 

+ 4 

10 

-4 

+ 

3 

-8 

+ 15 

3 + 2 

18 

+ 11 

40 

+ 8 

1 

+ 5 

11 
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+ 

3 

-4 

+ 8 

-2 + 2 

-20 

+ 7 

31 

+ 10 

8 

+ 2 

12 

13 

+ 

3 

-33 

+ 7 

8 + 3 

-17 

+ 7 

51 

+ 8 

-2 

+ 4 

13 

3 

+ 

3 

-46 

+ 14 

5+3 

-13 

+ 14 

40 

+ 14 

-36 

+ 4 

14 

9 

+ 

2 

00 

1 

+ 8 

6 + 2 

38 

+ 12 

75 

+ 9 

-21 

+ 3 












Mean.. 

4 

+ 

2 

-26 

+ 6 

4±2 

-2 

+ 8 

42 

+ 8 

-5 

+ 5 


ing for w may be partly due to the fact that the radiation pressure effects 
in the argument of perigee are too large to handle accurately. Also, I 
suspect that the anomalistic mean motion cannot be determined with 
sufficient accuracy for a satellite of such small eccentricity. This might 
be one reason why we have large discrepancies in the secular motions of 
the node. 

However, (O-C)’s in cb, Q and A, are still significant. 

(d) 1961 Nu — For this satellite precisely reduced Baker-Nunn observa- 
tions are not available and observations must be used that are not pre- 
cisely reduced. However, since the satellite is close to the earth and the 
inclination is the smallest used in this paper, the node and the perigee 
move rapidly and the relative accuracies in the determination of the 
secular motions are fair. 

Four sets of data are given in tables 7 and 8, which show a wide scatter 
in the values of (0-C) in A, and Ai. The residuals in the two secular 
motions take large values. This satellite was not used in the earlier 
determination of at that time the smallest inclination was 32°. 9, for 
1959 al. 

(e) 1961 Omicron — There are two separate satellites for 1961 o. How- 
ever, since they have almost identical orbital elements, they are treated 
as one satellite here. The eccentricity is very small. Since the inclina- 
tion is rather close to the critical inclination, the argument of perigee 
moves very slowly. Therefore, one set of observations must cover more 
than 500 days. However, as the mean motion changes rather rapidly due 
to air drag, I have used one set of 400-day observations. 
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Table 4— (0-C ) for 1959 Eta 


SATELLITE GEODESY 
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Table 6. — (0-C) for I960 lota 2 


ZONAL HARMONICS COEmOENTS 



±22 ±41 ±8 ±24 ±61 ±12 




















SATELLITE GEODESY 


For this satellite, the mean height is rather low, about 900 km, and 
the inclination is high. Therefore, the object is rather difficult to observe 
from the Baker-Nunn stations due to visibility conditions, and there are 
many gaps in the observations, periods for which accurate orbital ele- 
ments are not available. As the Baker-Nunn stations are between -b35° 
and -35° in latitude, the inclination of this satellite is poorly determined 
although the longitude of the node can be well determined. This situa- 
tion is contrary to that of Vanguard satellites. 


Table 8.— (0-C) jor I960 Nu 



«X10" 

nxio« 

A,X10“ 

^iXlOs 

il„Xl0* 

O 

X 

5 

1 

-59° + 15° 

211° + 30° 

-17± 3 

-173° +44° 

-45° ±20° 

17° ± 5° 

2 _ 

-74 +10 

90 +20 

-14+ 3 

-140 +30 

-50 ±40 

16 ± 6 

3 

-51 +15 

71 +80 

-17+3 

— 105 +35 

-25 ±34 

21 ±11 

4 

-10 +22 

111 +45 

11+8 

-20 ±24 

64 ±61 

0 ±12 

Mean 

-48 +20 

131 ±40 

-9 ± 14 

-no ±70 

-14 ±50 

14 ±10 


The secular motion of the node is determined quite accurately, as we 
can see in table 9. However, we cannot compute theoretical values of 
the secular motions so accurately as the observed ones, because of uncer- 
tainties in the inclination. Therefore, the standard deviations in (0-C) 
of S2 in table 10 are large. But (O-C)’s in il themselves are quite large, 
as we can see in table 10. In the previous determination of J accurate 
orbital elements from Baker-Nunn observations were not available. 

The value of (0-C) in 11 for the epoch 4 is quite different from the 
others, and I suspect this scattering is due to some accidental errors in 
io for the epoch 4, and give small weight to this value in taking the mean. 

For this satellite the radiation pressure effect in the argument of perigee 
is too large for my program to compute it with enough accuracy. This is 
also true for other satellites of small eccentricity. 

(f) 1961 Alpha Delta 1— This satellite has a polar orbit. However, 
as the mean height is quite high, we can determine the orbit very accu- 
rately from Baker-Nunn observations. 

This satellite, and the three listed in tables 13-17, which were launched 
in 1962, were not used in my previous determination. 

The first set of data is determined from 300-day observations, and the 
second set is from 400-day observations, which cover one revolution of 
argument of perigee. 

To compute the solar-perturbations there arise three small divisors, 
namely, 2(n.Q — di), 2(w— Q-j-riQ), and 2(ii0 — 2d>— fi). 

Tables 11 and 12 show that the eccentricity is very small and that 
(C-C) in n is very significant. 
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Table 11 . — Orbital Data for 1961 Alpha Delta 1 


SATELLITE GEODESY 




Table 14. — (0-C) for 1962 Alpha Epsilon 


ZONAL HARMONICS COEFnCIENTS 
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Table 16. — (0-C) for 1962 Beta Mu 1 


SATELLITE GEODESY 



156 ° ± 7 ° 18 + 4 - 51 ° ± 27 ° 4 ° + 20 ° 88 ° + 6 ' 







ZONAL HARMONICS COEFnCIENTS 


(g) 1962 Alpha Epsilon — For this satellite three sets of data are given 
in table 13. However, observations in sets 1 and 2 are overlapped widely. 
Since to and — Q have nearly the same value, 2(to+ fi) and 2(w0— 2co— fi) 
take small values, as for 1962 Beta Mu 1. Therefore we must be careful 
to compute luni-solar perturbation terms with such arguments. 

All values of (0-C) in table 14 are significant. 

(h) 1962 Beta Mu 1 — This is a geodetic satelhte, and although the in- 
clination is not very much different from that of 1956 ae, the eccentricity 
and the mean motion take quite different values. 

The mean height of this satellite is not high enough for the Baker-Nunn 
cameras to track the object over a long arc. Therefore the accuracy of 
determination of the orbital elements is not high. 

(i) 1962 Beta Upsilon — Unfortunately, precisely reduced Baker-Nunn 
observations are available for this satellite only for 200 days, during which 
the argument of perigee moves by 240°. Therefore I will increase by a 
factor of five the standard deviations given in table 17 in the determina- 
tion of J„. 

4. DETERMINATION OF J, 

Table 18 gives for the nine satellites the semimajor axes in units of 
earth equatorial radii, the inclinations, the eccentricity, and the area- 
to-mass ratio in cgs units. The same table also gives J 2 ^ terms and luni- 
solar secular terms in a> and 12 (Kozai, 1962; Kozai, 1959). 

A previous paper (Kozai, 1962) gives the formulas used to compute 
secular perturbations and amplitudes of long-periodic tenns with argu- 
ment 0 } by including up to 8th-order harmonics. However, I include up 
to 14th-order harmonics in the present determination, and the additional 
formulas are given in the following: 


10 fn(63- 109202-^49140^- 79560' -f-419961®) 

X(128-l-2304e2-i-6048e<-f3360e'-t-315e') 

~ 67 ; i^|64p>2 ^"(231 - 577502-1-392700^-106,5900' 

-I- 124,3550'- 52,00301“) • (256-h7040e2-l-31,680e*-f 36,960e' (10) 

-fll,550e'-|-693ei“) 

- ... ■ „ w0(429- 14,58602-1-138,56701-554,2680' 

-fl,O62,3470'-965,77O0i“-l- 334,305012) • (1024-b39,936e2 
-H274,560e‘-h549,120e'-l-360,360e'-f72,072ei“-f3003ei2), 
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5co= —dSQ 

3465/io 


i n(63 - 3465«2+3O,O3O0< - 90,0909*+ 109, 395«* 


8,388,608pw 
-46,189fl»'') • (128+1152e2+2016c*+840e*+63e«) 

~ 268,W56p^^ ~ 18,0189^+225,2259*- 1,021,0209* 

+2,078,5059*- l,939,9389i«+676,0399i*) • (1024+39,936e* 
+274,560e*+549,120e*+360,360e*+72,072ei«+3003e‘*) (11) 

- 4,2^S f -- ’‘**^~**-°**'”+™'"**^'~*'***’***'’‘ 

+ 14,549,5359*- 22,309,287910+ 16,900,97591*- 5,014,57591*) 
X(1024+19,968c*+91,520e*+137,280e*+72,072e* 

+ 12,012ei“+429ei*), 

s 

Se= —sin i (1 — 59*)~i(l— e*)^ CjAjBj sin to, (12) 

i-* 

Si= —ed 5e/{sin i (1— e*)}, (13) 

6 

6J2 = e9 sin“i i (1 — 59*)~i ^ Cj{ —sin* i*D/ 

^ +(9-59*)(l-59*)-i^lBjCOSco, (14) 


5to= — 95Q— sin i-e~i*(l — 59*)~i CjAjEj cos to, 


(15) 


where 

9= cos i 

„ 105J, 

65,536/2P^ 

1155Ju 

' 4,194,304Jsp» 

^ 3003*/ 13 

‘“67,108,864/jpii 

444 = 7 - 3089*+20029*-400496+24319* 

^s=21-13659*+13,6509*-46,4109*+62,9859*-29,3939io 

^6 = 33-29709*+42,0759*-213,1809*+479,6559*-490,3149'o+185,7259i* 
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54 = 64 + 336e2 + 280e^ + 35e« 

56=128+1152e2+2016e^+840e«+63e8 

56 = 512+7040e2+21,120e<+18,480e6+4620e*+231e“ 

54 = 88(7-9102+2730^-2210'’) 

5s= 130(21 -42O02+214204-38760'’+22610*) 

5e = 6O(99-28O502+21,3180''-63,9540'’+81,7190'*-37,1450'“) 

54 = 64 + 1776g2 +4760e' + 2485e« + 21 Oe* 

£5 = 128+5504e2+26,208e^+30,072e«+8,967e®+504e"’ 

56 = 512+31,360e2+232,320e^+467,280e'’+300,300e8+57,982e’» + 23106^2 

(16) 

(a) Even harmonics — Table 18 gives equations of condition to deter- 
mine values of through J u. There are 18 equations with 7 unknowns. 

The equations can be solved by assigning to each a weight reciprocally 
proportional to the standard deviation. Actually, each equation is 
divided by its standard deviation, and then normal equations are con- 
structed. Before solving the equations, note that S(0-C)2 is 3882 
( = 18X14.72); that is (0-C) is bigger than the standard deviation by 
factor of 14.7. This value comes down to 23 = (18 — 6 ) X 1.42 after solving 
Ji 2 , and to 13.4= (18 — 7) XI. 12 after solving J m, whereas it is 93.5 = 
(18 — 5)X2.72 after Jw is solved. Therefore we can stop either at J 12 or 
at J 14 , although the solution including Ju is, of course, better. 

In table 19 residuals based on the solutions up to Jh and J 12 are given 
under headings I and II, respectively, in units of 10“® degrees. Under 
the heading KH, residuals based on King-Hele and Cook’s values (1964) 
are given; that is, 

J 2 =1802.70X10-®, J 4 =-1.40X10-®, 

Jt = 0.37 X10-®, Jg = 0.07 X10-®, 

Jio= -0.50X10-®, Ji 2 = 0.31X10-®. (17) 

In the node equations the residuals based on my new determinations 
for 1962 fiv are larger than the standard deviations. However, since this 
datum is not entirely reliable, being based on a single determination cov- 
ering an incomplete period of time, this may not be a weak point in this 
determination. 

In the perigee equations of 1961 v and 1962 at, the residuals are larger 
than their standard deviations. This may suggest that we must still 
include higher-order terms to express these data. 
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The two sets of solutions derived are the following: 


Solution I (in units of 10 


dJ2=1.65, dJi^ 1.81, 

dJo = 2.56, 

II 

00 

-2.50, 


±6 ±16 

+ 30 


+ 50 

(18) 

J 10— — 0.54, 
+ 50 

IsS 

II 

1 

1+ 00 

J u — 

1.79, 
+ 63 


Solution II (in units of 10 ’’) 

dJi = 1.50, 

dJ 4 — 2.03, 

II 

<o 

2.03, 


±5 

±18 


+ 31 

(19) 

dJi = -1.29, 
+ 34 

J 10 — 1.55, 

±45 

J \1 — 

-2.94. 
+ 49 



(b) Odd harmonics — As shown in table 20, we have 32 equations to 
determine 6 unknown coefficients of odd harmonics. At first 2(0-C)^ 
is 349( = 32X3.3^). This number comes down to 153( = 28X2.3^) after 
J, is solved, and to 42( = 27X1.252) and to 39( = 26X1.232) after Jn and 
Ji 3 , respectively, are solved. Therefore, the inclusion of Ju does not 
reduce the residuals too much. Two sets of solutions are derived, one 


up to Ju and one up to Ju; that is. 




Solution I (in units of 10“^ 




dJs = 0.31, dJ,=-\Al, 

II 

1— 

1.36, 


±20 ±25 


±39 


J9=-1.67, Jii= 3.02, 

J n — 

-1.14, 


±60 ±35 


±84 


Solution II (in units of 10“^) 




dJz= 0.07, dJ6=-1.22, 

II 

0.93, 


±11 ±17 


±22 


Jo =-0.75, 

J ii — 

2.96. 

(20) 

±17 


±35 



Table 20 gives the residuals based on solutions I and II for each datum. 
Residuals in the eccentricities of 1961 v and 1962 in the perigee of 
1961 V, and in the nodes of 1962 a€ and 1962 fiv have much larger values 
than the standard errors. This may show that still higher-order har- 
monics are significant. 


74 



. ZONAL HARMONICS COEITiaENTS 

In this analysis parallactic terms 

are neglected in computing lunar 

p>erturbations. However, in the parallactic disturbing function there is 

a term, 


45 . . 5 5 Q 

-g I'Sin €(1— -sin* i)(l— - sin* e)ee'(H-^*) sin wsin u', (21) 

where e is obliquity, e is lunar eccentricity, and w' is lunar argument of 
perigee. Since u' moves slowly, we must include this term if we treat 

observations of high-altitude satellites in the future. 

5. RESULTS 

The two sets of solutions derived in 

this paper are the following: 

Solution I (units of 10“*) 


J2 =1082.645, 

J, = -2.546, 

±6 

±20 

Ji = -1.649, 

Js =-0.210, 

±16 

±25 

Je = 0.646, 

Jy =-0.333, 

±30 

±39 

Js = -0.270, 

Js =-0.353, 

±50 

±60 

Jio= -0.054, 

Jji= 0.302, 

±50 

±35 

J12— — 0.3o7, 

Ji3= -0.114, 

±44 

±84 

Ji4= 0.179 


±63 

(22) 

Solution II 


Ji =1082.630, 

Jz =-2.559, 

±5 

±11 

Ji =- -1.627, 

Js =-0.185, 

±18 

±17 

Js = 0.593, 

Jy =-0.376, 

±31 

±22 
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Js = 

-0.149, 

±34 

J, = 

0.039, 

±17 


Jio = 

-0.155, 

±45 

Jn = 

0.296, 

±35 


J 12 = 

-0.294 

±49 



(23) 


A. H. Cook (1964) recently derived values of J 2 , ^4 and Je by using high 
satellites only, and his results show remarkable agreement with Solution I. 

The flattening of the reference earth ellipsoid based on this value of 

is 1/298.252. The theoretical value of Ji for the reference ellipsoid 
assumed to be in hydrostatic equilibrium is computed as —2.350X10"®. 
The deviation of the geoid computed on the geopotential based on solution 
I is expressed as a function of geometric latitude : 

fi=+0.8-18.3 sin /S-87.8 sin* ;8-119.1 sin® |3+1042.5 sin® 

+ 1191.7 sin® /3-5074.2 sin® ;3-3636.7 sin* /3+12, 668.0 sin® ^ 
+5230.8 sin® 16,676.3 sin'“ ;3-3556.4 sin“ /3+10,913.0 sin'* |3 
+926.8 sin'® |3— 2791.3 sin'® /3 (in meters). (24) 

Figure 1 shows the value of as a function of jS based on this equation. 
The value of geoid height h in the north pole is 13.5 meters, which is the 
maximum value, and is —24.1 meters in the south pole. 

In the solutions (22) and (23), the values of do not tend to converge 
to zero as n increases. However, if n is large enough, J „ should take a 
very small value. Otherwise the gravity expression, which is derived 
by differentiating the potential with respect to the radius, may give a 
very great difference of gravities between the equator and the poles and 
between the north and south poles. 

To determine how strong or weak the solutions (22) and (23) are, the 
correlation coefficients in my determinations are shown in tables 21 and 
22. The tables indicate that these solutions are derived from rather 
strongly correlated equations of condition. Therefore, in the future we 
must use both low and high satellites having the same inclination. 

However, to determine the orbital elements of low satellites with high 
inclinations we need observations from high latitudes. As I mentioned 
earlier, I could not assign a large weight to the node equation of 1961 o 
to determine even-order coefficients, because the inclination could not 
be determined with sufficient accuracy. Also, I must mention that I did 
not use satellites with inclinations below 28°, between 50° and 67°, or 
between 67° and 85° in this determination. 
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Figure 1 . — Geoid height (A) as a function of geometric latitude (J>). Solid line 
shows geoid height in Northern Hemisphere and broken line shows that in the 
Southern Hemisphere. 


However, I believe that the present determination is much more re- 
liable than the previous one, since the data themselves are more reliable, 
both because of the number of observations and because I included some 
satellites that were not used in the previous determination. 
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Part 3 

Derivation of the Earth's Gravity Field by 
Nonoptical Tracking 


Although precision optical tracking has yielded the 
most extensive geodetic information to date, electronic 
tracking methods are rapidly becoming competitive. 
Of the latter methods, the precise measurement of the 
frequencies received from stable oscillators in the 
satellite has been particularly fruitful. 

A section of a recent report by Dr. Robert R. Newton 
of the Applied Physics Laboratory, one of the group 
which has exploited this technique most successfully, 
provides an explanation of the basic technique. This 
is followed by a report of the analysis of some of the 
Doppler data and data from Syncom 2. 


* PXGI BLANK WOT RUBBD. 


The follomng is excerpted from a paper presented by Robert 
R. N ewton at the European Symposium on Satellite Geodesy, 
Paris, December 1964. 


Measurements oj the Doppler Ship in Satellite 
Transmissions and Their Use m Geometrical 
Geodesy 

Robert R. Newton 
Applied Physics Laboratory 
Johns Hopkins University 


1. CHARACTERISTICS OF EXISTING DOPPLER MEASURING EQUIPMENT 

T he applied physics laboratory exercises technical direction on 
behalf of the Department of the Navy of a network of doppler 
tracking stations that is called the TRANET System. This system has 
been in operation for approximately five years and currently contains 
twelve stations with another expected to be added in 1965. Many other 
sites have been occupied at various times. The Department of the Navy 
maintains this tracking system and supports geodetic analysis of the data 
obtained from it as a part of the national geodetic satelhte program.* 

Up to the present, the data from this system have been used almost 
entirely in studying dynamics of satellites and in deriving models of the 
earth’s gravity field from these dynamic studies. 

The method of operation of the existing doppler stations is described 
in Ref. 1 . Briefly summarized, a station makes a measurement of doppler 
frequency eveiy 2 or 4 seconds and thus typically obtains several hundred 
data points during each pass of a satellite. 

In these measurements the effects of radio noise and of instrumental 
contributions to error are virtually neghgible. The dominant errors in 
the raw data arise from refraction both in the ionosphere and in the 
troposphere. 

Of these two sources of refraction ionospheric refraction is by far the 

*The station in Antarctica, if it is established, wiU be supported by the National 
Science Foundation, and therefore is not strictly speaking a part of the TRANET 
System. However, data from it will be available for analysis so that it can be con- 
sidered together with the other stations for the purposes of this paper. 
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biggest problem. If no correction were made for refraction in the iono- • 
sphere at a transmitted frequency of 162 mh, for example, satellite posi- 
tions derived from the doppler measurements could easily be in error by 
as much as a kilometer. Fortunately, the refractivity in the ionosphere 
depends upon the frequency and to first approximation varies inversely 
with the square of the frequency. This fact is the reason for using two 
or more transmitted frequencies (see Refs. 1 and 2) which are both con- 
trolled from the same oscillator and which are hence always coherent. By 
measuring the apparent doppler shift at each of two coherent frequencies, 
it is possible to calculate the effects of ionospheric refraction and therefore 
to eliminate refraction, on the assumption that the refractivity does 
indeed vary with the square of the frequency. Since this law does not 
hold exactly, there is some residual refraction arising from neglected 
effects. At middle and high latitudes we believe that the residual effects 
amount to the equivalent of a few meters or less in satellite position. 
However, there is some indication that at low latitudes and under some 
circumstances the residual effect may amount to 100 meters or more. 
Fortunately, position errors from residual refraction have a tendency to 
average out when many satellite passes are used. 

Depending upon the level of accuracy desired, the uncorrected iono- 
spheric refraction may still be unacceptable. There are three ways of 
reducing the effect of refraction still further. One is by the use of a third 
frequency which will be provided on GEOS A (Ref. 3). A second method 
which shows promise involves a calculation of the higher order refraction 
effects based upon a model of the ionosphere, combined with the measured 
refraction effect obtained from two frequencies. There is considerable 
hope that the higher order refraction can be calculated in this way to an 
accuracy of about 10% of the high order effects. The third way is to 
use higher frequencies such as the 972 mh beacon provided on GEOS A. 

Tropospheric refraction effects can produce the equivalent of about 
50 meters in satellite position if no correction is made. Tropospheric 
refraction does not depend appreciably upon frequency and therefore 
cannot be removed by the techniques that are successful with the iono- 
sphere. Up to the present time the only way that we have found to 
correct for tropospheric refraction (Ref. 4) is to calculate this refraction 
using a model of the troposphere combined with suitable meteorological 
observations made at the time of each pass. This correction has been 
quite successful except when a weather front is near a station at the time 
of a pass. If high accuracy is desired it is necessary to eliminate any 
data obtained when a weather front is within about 30 km of a station. 

As we have already said, the stations that make up the TRANET Sys- 
tem supply several hundred data points during each pass. This volume 
of data is necessary for many research purposes such as research into the 
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effects of tropospheric refraction which were just discussed. However, 
if the presently used methods for ehminating refraction are considered • 
by an observer to be sufficiently accurate, and if his principal interest is 
in using the doppler data for geometric purposes, he can work with much 
simpler receiving equipment which provides a much smaller volume of 
data. 

This simplified equipment is now under development at the AppUed 
Physics Laboratory. 
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The following paper vxis originally published in Nature, 
vol SOO, No. 4902, October 12, 196S, pp. 124-125. 
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* 

Determination of the Non-Zonal Harmonics of 
the Geopotential From Satellite Doppler Data 

W. H. Guier 
Applied Physics Laboratory 
Johns Hopkins University 


A LEAST SQUARES DETERMINATION of valucs foT the Coefficients, 
C”, S", of the non-zonal harmonics of the geopotential : 

=~'(l+ 2 S — - [C^ cos mX+(Sr sin mX]l (1) 

' t n-2 m—0 r j 

has been made through n = m = i using Doppler data from the three 
Satellites 1961 ol, 1961 017!, and 1962 /3/xl. In equation (1): 

r = geocentric radius from the Earth’s centre of gravity; #5 = geocentric 
latitude; X= geocentric longitude from Greenwich meridian, 

and P’!Xz) are the associated Legendre polynomials: 

— d™ 

p:(z) = (l-,y (2) 

The values of C J and S' were constrained to be identically zero in the 
determination. 

Table 1 presents a brief summary of the Doppler data which were 
chosen for the determination. The data were chosen as being the best 
available from the point of view of: (1) quality with respect to noise- 
level, timing errors, presence of ionospheric refraction errors, and stability 
of the satellites’ transmitted frequencies ; (2) uniform distribution of sta- 
tion locations over the Earth’s surface with particular care taken to avoid 
an overweight of data from stations located within the continental 
United States; (S) uniform coverage of the satellite orbit by the data with 
respect to the angle between the ascending node and the Greenwich 
meridian. 

The principal ionospheric refraction contribution was removed by 
coherently combining two frequencies received simultaneously from the 
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satellite in the Doppler receiving equipment*’^. To ensure further mini- 
mum refraction errors, data were considered only if the maximum eleva-- 
tion of the pass was less than 83° and greater than 20°. Finally, data 
points which corresponded to instantaneous elevations below 13° were 
deleted, and the remaining data points corrected for tropospheric refrac- 
tion’. A typical noise-level for that data finally chosen was roughly 5 
parts in 10‘° of the transmitted frequency. Timing errors were judged 
to be below 0-003 sec (3-sigma value). 

The data which survived the foregoing criteria were then grouped by 
days, one group containing all passes for a given satellite which occurred 
between two successive midnights (g.m.t.). After grouping, only those 
were selected which contained ten or more high-quality sets of data from 
at least five separate stations with at least three being outside the conti- 
nental United States. From Table 1 it can be seen that the fewest groups 
of data for any one satellite occurred for Satellite 1961 aijl. This is due 
to the low inclination of the satellite precluding reception of data from 
high-latitude stations. While considerably more groups of data satis- 
fying the foregoing criteria were available for the other two satellites, the 
total number of groups chosen for each satellite was limited to being 
roughly equal to avoid overweighting the data corresponding to a particu- 
lar inclination. In the final choice of groups care was taken to ensure a 
reasonably even distribution of angles between the ascending node and 
the Greenwich meridian at the time of the first pass in each group. In 
addition, groups for two different satellites occurring on the same day 
were avoided to increase the a priori probability that the effects of unde- 
tected station instrumentation errors and refraction errors were random. 

Having made the final choice of groups, a least-squares determination 
of a satellite orbit for each group was made, resulting in no computed arc 
of the satellite motion lasting more than 24 h. Orbit parameters span- 
ning data over such short time-intervals were chosen because: (1) the 
satellite trajectory could be integrated numerically with negligible error 
and made unnecessary the very tedious task of solving the perturbed 
equations of motion to high orders ; (2) unknown fluctuations in atmos- 
pheric drag and errors in the values for the zonal harmonic coefficients 
introduce negligible errors. 

The computer programme used to determine the orbits^ includes the 
Sun and Moon forces, although they should introduce negligible errors 
for these short data spans. The resulting differences between the 
Doppler data points and the corresponding computed values for the 
Doppler shift (based on fixed initial values for the harmonic coefficients 
of the geopotential) were used as the final data for the least-squares 
determination of the differential changes in the values of the non-zonal 
harmonic coefficients. 
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NON-ZONAL HARMONICS OF THE GEOPOTENTIAL 


Average time 
span per 
group (h) 

16-9 

22- 7 

23- 3 

Average No. 
passes per 
group 

^ M O 
^ ^ 

No data groups 
(satellite 
orbits) 

8 

9 

10 

27 

No. distinct 
station loc. 

Ci 40 CO 


No. individual 
satellite 
passes 

c: 40 

« Cl 40 

441 

Eccentricity 

0-010 

0-007 

0-008 

Totals., 

Semi-major 
axis (km) 

7,410 

7,510 

7-320 

Inclination 

(deg) 

32-4 

50-1 

66-8 

Satellite 

designation 

6 0^ O 

*-H ^ 

cO CD CD 
Oi Oi O 
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Very briefly, the least-squares procedure used for determining the non- 
zonal harmonic coefficients was the following. The differences between 
the theoretical and experimental values for the Doppler shift were as- 
sumed to arise wholly from errors in the locations of the tracking stations 
and values for the non-zonal harmonic coefficients. It was further 
assumed that because of these errors there existed errors in the orbit 
parameters associated with each group of data. Perturbation equations 
were then derived to express these residuals as a function of differential 
changes in the three co-ordinates of each tracking station represented in 
the totality of data, differential changes in each of the six orbit param- 
eters for each group of data, and differential changes in the values of 
the C" and S" of equation (1). A special IBM 7094 computer pro- 
gramme was designed and coded in which the differential values of the 
station positions and orbit parameters are continually adjusted to yield 
an extremum of the mean-square data residuals with respect to these 
parameters (irrespective of whether the resulting values are reasonable) 
for any given values for the differential harmonic coefficients. Conse- 
quently, at any point in the computing process, changes in the values of 
the differential station and orbit parameters can produce negligible 
reduction in the mean-square data residuals. In this way the station 
co-ordinates and orbit parameters are effectively eliminated from the 
least-squares determination of the harmonic coefficients. The pro- 
gramme is then iterated until stable values for the differential non-zonal 
harmonic coefficients are obtained. 

Table 2 presents the results of the determination using the data sum- 
marized in Table 1. Without new data at a different inclination, it is 
very difficult to estimate the true error in these values ; and no attempt 
has been made to ascribe a realistic error to each coefficient. The fact 
that the station parameters were effectively eliminated from the deter- 
mination ensures that the principal error in the resulting values for the 
coefficients arises from correlated errors in the experimental data and not 
from erroneous values for the station locations. Consequently, the 
efficacy of this method depends principally on there being sufficient ex- 
perimental data of sufficient accuracy that the build-up of computer error 
combined with experimental errors does not destroy the sensitivity of 
the data residuals to the harmonic coefficients once the station and orbit 
differential parameters are eliminated. Several numerical experiments 
were performed with the amount of data intentionally reduced to the 
point where a realistic solution would not be expected. The objective of 
the experiments was to obtain criteria for the number of harmonic coeffi- 
cients which could be evaluated with the data available. As expected, 
these experiments indicated that if stable values for the coefficients were 
produced with successive iterations, and if their resulting values did not 
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change markedly with changes in the relative weighting of the data 
residuals, the final values should be trustworthy. 

Recognizing that the Tesseral harmonic, n = 2, wi = l, should have 
neghgibly small coefficients, a final check on the validity of the values in 
Table 2 was made by performing a new determination wfith C\ and S 2 
added to those coefficients fisted in Table 2. The resulting values for 
those coefficients in Table 2 were not materially altered and the values 
for the new coefficients were found to be : 

C 2 = 1-78X10-* 

S2=-3-48X10-* (3) 

The fact that they are negligibly small (compared with the values in 
Table 2 for m=l) lends considerable credence to the values given in 
Table 2. 

A comparison between the values given in Table 2 and the present 
literature*-* indicates agreement of roughly the same character as the 
agreement between any two previously published sets of values. Conse- 
quently, such comparisons add little additional information on the 
accuracy of the values presented in Table 2. 

Clearly the overall effort involved in obtaining these results was very 
large, and it is impossible to acknowledge all the people who made valu- 
able contributions. However, I wish to acknowledge especially the 
contributions of R. R. Newton and S. Yionoulis for their aid in develop- 
ing the satellite perturbation equations; G. Worsley, G. C. Weiffenbach, 
P. E. P. White, and their colleagues in connexion with obtaining the 
experimental Doppler data; and H. D. Black, C. Weisert, R. Henderson, 
and their colleagues for their help in coding and running the complex 
computer programmes involved. 

This work was supported by the Department of the Navy, Bureau of 
Weapons, under contract NOw 62-0604-c. 
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Most geodetic data have been derived from comparatively 
near-Earth satellites. However, some of the terms in the 
Earthl s gravity field are masked for these satellites by various 
diurnal effects. Conversely, a satellite in a Zl-hour orbit is 
particularly sensitive to these terms because it is constantly 
affected by the same perturbations. Because of its high 
altitudes, a 2^-hour satellite is difficult to track optically and 
the Minitrack interferometer does not have the angular accu- 
racy needed. To overcome this problem, an electronic trans- 
ponder has been developed which provides a measure of both 
the distance to the satellite and its radial velocity. This range 
and range^aie system was used on the Syncom satellite. 
The geodetic analysis of the tracking results is reported in the 
following paper which was originally published in the Journal 
of Geophysical Research, vol. 70, no. 6, 1965, pp. 1566-1568. 


A Determination of Earth Equatorial Ellip- 
ticity From Seven Months of Syncom 2 
Longitude Drift 

C. A. Wagner 

Goddard Space Flight Center., NASA 


T he 24-houk Syncom 2 satellite has been under periodic observation 
by range and range rate radar and Minitrack Radio Interferometer 
stations since mid-1963 (Wagner, 1964b). Seven months of longitude 
drift in the vicinity of two momentarily stationary configurations were 
analyzed for sensitivity to hypothetical longitude components of the 
earth’s gravity field which would be in “resonance” on such a satellite 
(Blitzer et al., 1962; Wagner, 1964a). This drift, in the region 54 to 
64° W (over Brazil), was derived from orbits calculated at the Goddard 
Space Flight Center. 

From mid-August 1963 to late November 1963, the figure-eight ground 
track of Syncom 2 drifted from 55° W to 59° W, with a mean acceleration 
of 

— (1.27 + .02) X 10~* degrees/day’ (1) 

The average growth of the semimajor axis for this period was estimated 
as (0.0993 ± .0042) km/day. The figure-eight configuration was mo- 
mentarily stationary at about 54.76° on September 6, 1963, at which 
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time the semi-major axis was (42166. 0 + . 2) km. On November 28, 
1963, the westward drift of Syncom 2 was stopped by ground comman,d 
firing of tangentially oriented cold gas jets on-board the satellite. 

From early December 1963 to mid-February 1964, the figure-eight 
ground track of Syncom 2 drifted from a momentarily stationary posi- 
tion at about 59.15° W (semimajor axis: 42165. 9 + . 4 km) to 63.5° W. 
The estimated mean geographic longitude acceleration of the ground 
track for this period was 

— (1.32 + .02) X 10“® degrees/day- (2) 

The average growth of the semimajor axis for this period was (0.0994 + 
.0080) km/day. Simulated Syncom 2 trajectories for these drift periods, 
starting with the initial orbital elements, show good agreement with the 
observed trajectories if an earth gravity field with longitude dependence 
is used in the particle program of the simulation. Sun and moon influ- 
ence, as well as earth zonal gravity influence, on drift acceleration (as 
assessed by this simulation) for the full seven months of Syncom 2 data 
analyzed appears to be negligible compared with hypothetical earth 
gravity with longitude dependence. Other possible causes of this ob- 
served long-term accelerated longitude drift, such as (1) selective out- 
gassing or leaking of on-board gas jets, (2) micrometeorite collisions, 
(3) solar wind or radiation interactions, or (4) geomagnetic field interac- 
tions in the environment of Syncom 2 all appear to be extremely unlikely. 

Wagner (1946b) showed that any ellipticity of the earth’s equator will 
cause the otherwise stationary figure-eight ground track of a 24-hour 
near-circular-orbit satellite to drift in longitude with an acceleration 
given by 

X= — d .22 sin 2y (3) 

where 

^ 22 = -727tV22 {Ra/af rad./sid. day- (4) 

Relations 3 and 4, without the inclination factor, can be derived from 
equation 14 in Allan (1963) and equation 57b in Wagner 1964a. 

7 is the nodal longitude of the 24-hour configuration east of the minor 
axis of the elliptical equator. J 22 is the amplitude of the first significant 
longitude-dependent term in the spherical harmonic expansion of the 
earth’s gravity potential. It is related to the difference between major 
and minor earth equatorial radii by (Izsak, 1961; Wagner, 1962; Kaula, 
1965) 

do — 5 o= — QRoJ 22 

Ro is the mean equatorial radius of the earth. Small a is the semimajor 
axis of the 24-hour satellite. Small i is its inclination. For Syncom 
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. 2, the inclination during both drift periods analyzed was close to 33.0°. 

Under the assumption that the above observed drift acceleration is 
sensing only the lowest order of longitude-dependent earth gravity (that 
associated with the ellipticity of the earth’s equator), the two drift 
accelerations of Syncom 2 (equations 1 and 2), satisfy (3) and (4) uniquely 
with the following values of the magnitude (J^) and phase angle (X 22 ) 
of equatorial ellipticity: 


j. ( unadjusted for higher order\ 

■'-( earth pavitykctB ) — a ™± 05)X10- (5) 

corresponding to a difference in major and minor equatorial radii of 


and, 


do— l>o = 65±2 meters 
X 22 = (19 + 6) degrees W’est of Greenwich 


(6) 


(locating the major axis of the earth’s elliptical equator). 

From Table 1 it is seen that these values for the second-order tesseral 
harmonic of earth gravity are in reasonable agreement with recent 
determinations of longitude gravity from lower-altitude satellite observa- 
tions and surface gravimeter data. From a consensus of recent and older 
geoids, which give tesseral field coefficients to higher order than the 
second, it appears that at the high altitude of the 24-hour satellite the 
second-order longitude effect accounts for about 85% of the full field 
effect at 54° to 64° W. The full longitude field is consistently depressed 
below the ^22 field at these longitudes for all the recent geoids of Table 1. 
On the basis of the recent measures of higher-order tesseral gravity, cal- 
culations show the Syncom 2 estimates 5 and 6 should be ; 


or, 


( adjusted for probable higher-orderX 
earth gravity effects ) 


(1.92±.2)X10-« 


and 


flo— 6o = 73 + 8 meters 
X 22 = (2 1+7)° west of Greenwich 


(7) 

( 8 ) 


That longitude-dependent earth gravity exists seems well established 
from a large number of recent gravity reductions on different bases 
(Table 1). This 24-hour satellite reduction (because the altitude is so 
high) appears to separate out the second-order effect almost entirely 
from the sum of all higher-order earth gravity effects. Unless the earth 
is far more inhomogeneous in longitude than is thought to date, (7) and 
(8) are to be considered absolute estimates of the bounds on the ellipticity 
of the equator. 
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Table 1 illustrates the present state of fluidity with respect to individual 
tesseral coeflScients. The adjusted bounds on the Syncom 2 determina- 
tion of equatorial ellipticity encompass “reality” (and not just the ran- 
dom errors of the experiment) as far as can be judged by a number of 
recent geoids based on different kinds of gravity data. A recent reduc- 
tion of Syncom 2 drift over the Central Pacific between 120° and 170° 
west longitude appears to confirm this conjecture. Preliminary results 
from the accelerated drift in this region yield: 

/m=-(1.75)X10-« 

\ 22 = (18)° west 

The experiment-standard deviation on the J 22 number, however, is a 
number of times larger than that on the J 22 value derived in the Syncom 
2-Brazil experiment (eq. 5). 
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Part 4 

Astronomical Constants 


Astronomers are very conservative about changing the 
basic constants which they use in their various ephem- 
erides and in many of their analyses. It is not that 
they do not recognize the improvements which take 
place from time to time, but they feel that the confusion 
resulting from the use of frequently changed values 
would more than offset the advantages of using improved 
values. Nevertheless, the rapid improvement in our 
knowledge of the size, shape, and mass of the Earth, 
and of the mass of and distance to the Moon and Venus, 
which has been derived from satellites and probes, as 
well as from modern radar observations, has been so 
significant that the following report was approved by 
the International Astronomical Union in 1964. 
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Report of the Working Group on the System oj 
Astronomical Omstants: Agenda and 
Draft Reports 

International Astronomical Union, Twelfth General Assembly 
25 August-3 September 1964 


INTRODUCTION 

T he principal recommendations of this report on the system of 
astronomical constants are in accordance with the resolutions 
passed at lAU S 3 unposium No. 21 (Paris, May 1963). We first of all 
give a reference list of the constants of the system and a set of explanatory 
notes. (See Appendix B.) We have used the term “primary constant,” 
rather than “fundamental constant,” since the latter has a connotation 
in astronomical usage that is inappropriate to the manner of selection of 
the primary constants. In choosing the values for the primary constants 
we have, perhaps, adopted a conservative view of the likely errors of 
their determinations, but even so the new system should be of adequate 
accuracy for astronomical studies for many years. Limits within which 
the true values are believed to lie are indicated in a later section, in which 
we also give expressions for differential corrections to the derived con- 
stants. Finally we discuss the manner in which this system should be 
introduced into the national and international ephemerides. 

We regret that owing to severe illness Professor A. Danjon has been 
unable to share in the preparation of this report ; Dr. J. Kovalevsky has, 
however, been co-opted to the Group in his stead. 

NOTES ON THE CONSTANTS 

1. The value given for the number of ephemeris seconds in the tropical 
year at 1900 is taken from the definition of the ephemeris second that 
was adopted by the Comit4 International des Poids et Mesures {Proc^ 
Verbaux des Seances, deuxieme scrie, 25, 77, 1957). It is, in fact, derived 
from the coefficient of T, measured in Julian centuries of 36525 days, in 
Newcomb’s expression for the geometric mean longitude of the Sun 
referred to the mean equinox of date. In the list “1900” refers to the 
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fundamental epoch of ephemeris time, namely 1900 January 0 at 12'' 
E.T., or to 1900.0, as appropriate; the values for constants 20-23 also . 
refer to the fundamental epoch. Throughout the list and this report the 
term “second” must be understood to mean the “ephemeris second.” 

2. The value of the Gaussian gravitational constant (k) is that adopted 
by the lAU in 1938, and serves to define the astronomical unit of length 
(a.u.) since the corresponding (astronomical) units of mass and time are 
already defined. (The unit of mass is that of the Sun and the unit of 
time is the ephemeris day of 86400 ephemeris seconds. The units of k 
are: (a.u.)* (ephemeris day)”' (Sun’s mass)”*). To simplify the later 
equations an auxiliary constant k', defined as /c/86400, is introduced and 
a rounded value is given in the list. 

3. The value for the measure of the a.u. in metres is a rounded value 
of recent radar determinations. 

4. The value for the velocity of light is that recommended by the Inter- 
national Union of Pure and Applied Physics in September 1963. 

5. The term “equatorial radius for Earth” refers to the equatorial 
radius of an ellipsoid of revolution that approximates to the geoid. (See 
also note 16.) 

6. The term “dynamical form-factor for Earth” refers to the coefficient 
of the second harmonic in the expression for the Earth’s gravitational 
potential as adopted by lAU Commission 7 in 1961. (See also note 16.) 

7. The geocentric gravitational constant (GE) is appropriate for use 
for geocentric orbits when the units of length and time are the metre and 
the second; E denotes the mass of the Earth including its atmosphere. 
Kepler’s third law for a body of mass M moving in an unperturbed 
elliptic orbit around the Earth may be written 

GE{\-\-M / E) =nV 

where n is the sidereal mean motion in radians per second and a is the 
mean distance in metres. The value of GE is based on gravity measure- 
ments and observations of satellites. 

8. Again the mass of Earth includes the mass of the atmosphere. The 
reciprocal of 81 ‘30 is 0*0123001. 

9. The value for the sidereal mean motion of the Moon is consistent 
with the value of the tropical mean motion used in the improved lunar 
ephemeris, less the general precession in longitude. 

10-12. The values of the principal constants defining the relative posi- 
tions and motions of the equator and ecliptic are those in current use. 
Secular terms and derived quantities are already tabulated elsewhere. 

13. The rounded value 8" *794 for the solar parallax should be used 
except where extra figures are required to ensure numerical consistency. 
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14. The value of the light-time for unit distance is numerically equal 
tg the numter of light-seconds in 1 a.u. Its reciprocal is equal to the 
velocity of light in a.u. per second. 

15. Apart from the factor Fi the constant of aberration is equal to the 
ratio of the speed of a hypothetical planet of neghgible mass moving in 
a circular orbit of unit radius to the velocity of light; it is conventionally 
expressed in seconds of arc by multiplying by the number of seconds of 
arc in one radian. The factor F\ is the ratio of the mean sjieed of the 
Earth to the speed of the hypothetical planet and is given by 

E, _”o 

k' (l-e2)i 

where Wg is the sidereal mean motion of the Sun in radians per second, Ug 
is the perturbed mean distance of the Sun in a.u., and e is the mean 
eccentricity of the Earth’s orbit. Newcomb’s values for Wg, Og and e are 
of ample accuracy for this purpose. The factor Fi and the constant of 
aberration take the following values 


F 1 k' 

1800 1.0001427 20.49583 

1900 1.0001420 20.49582 

2000 1.0001413 20.49581 


The rounded value 20"’ -496 should be used except where the extra figures 
are required to ensure numerical consistency. 

16. The condition that the reference ellipsoid of revolution for the 
Earth shall be an equipotential surface implies that three parameters are 
sufficient to define its geometrical form and external gravitational field, 
provided that the angular velocity (to) of the Earth and the relative mass 
of the atmosphere (/!„) are assumed to be known. The variability of the 
rate of rotation of the Earth can be ignored, and the mass of the atmos- 
phere is only just significant ; the required values are : 

co = 0'000072 921 radians per second; aio = 0-000001 

The expressions for the flattening (/) and the apparent gravity at the 
equator (gfe) in terms of the primary constants are, to second order: 

g, = {GE / (1— J^m+^ m^) 

where m = is obtained by successive approximations. The new 

values of these constants are not intended for geodetic use. 

17. The heliocentric gravitational constant corresponds to GE, but is 
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appropriate for heliocentric orbits when the units are the metre and the 
second. 

18-19. The derived values of the masses of the Earth and of the Earth 
+Moon differ from those currently in use, but will not supersede them 
completely until the system of planetary masses is revised as a whole. 

(See note 24.) ■ ■ ■ e 

20. The perturbed mean distance of the Moon is the semi-major axis of 
Hill’s variational orbit, and differs from that calculated from Kepler’s 
law by the factor Ej, which depends on the well-determined ratio of the 
mean motions of the Sun and Moon. (E. W. Brown, Mem. R.A.S., 

58, 89, 1897.) , 

21. The constant of sine parallax for the Moon is conventionally ex- 
pressed in seconds of arc by multiplying by the number of seconds of arc 
in one radian. The corresponding value of tt^ itself is 3422" '608. 

22. The constant of the lunar inequality is defined by the expression 
given and is conventionally expressed in seconds of arc. 

23. The constant of the parallactic inequality is defined by the expres- 
sion given; the coefficient is consistent with the corresponding quanti- 
ties in Brown’s Tables. 

24. The system of planetary masses is that adopted in the current 
ephemerides and the values given for the reciprocals of the masses include 
the contributions from atmospheres and satellites. The value for xsep- 
tune is that adopted in the numerical integration of the motions of the 
outer planets; the value used in Newcomb’s theories of the inner planets is 
l‘)700. In planetary theory the adopted ratio of the mass of the Earth 
to the mass of the Moon is 81-45 (compared with 81-53 in the lunai 
theory), and the ratio of the mass of the Sun to the mass of the Earth 
alone is 333432. This system of masses should be revised within the 
next few years when improved values for the inner planets are available 
from determinations based on space-probes. 

CORRECTION FACTORS AND LIMITS 

To first order, relative errors of the derived constants are given by: 
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The true values of the primary 
the following limits 

A : 149597 to 149601 X 10® m 

c ; 299792 to 299793 X 10® ms-> 

a,: 6378080 to 6378240 m 

J 2 : 0-0010824 to 0-0010829 

GE: 398600 to 398606X10* m®s-® 


constants are believed to lie between 
81-29 to 81-31 

n%: correct to number of places 
given 

p: 5026r40 to 5026^90 
€: 23°27'08fl6 to . . . 08."'36 
N: 9^200 to 9^210 


Correspondingly the limits for the derived constants are: 


yrQ : 8r79388 to 8T79434 

499?001 to499K)16 

k: 20f4954 to 20^4960 

GS: 132710 to 132721 
X 1015 m® S-® 

S/E: 332935 to 332968 

S/£'(1+m) : 328890 to 328922 


/-i: 298-33 to 298-20 
a*: 384399 to 384401 X 10® m 
sin 7T,: 3422f397 to 3422f502 
L : 6^4390 to 6T4408 
P,: 124^984 to 124T989 


IMPLEMENTATION OF THE REPORT 

We regard it as essential that the new system should be introduced into 
the national and international ephemerides as soon as possible after its 
adoption. Accordingly we are requesting the directors of the principal 
ephemeris offices to study the consequences of the introduction of the new 
system so that a firm timetable can be drawn up at the meetings of Com- 
mission 4 in Hamburg. We provisionally suggest that the new system 
be introduced into the almanacs for the year 1968; for the Sun, Moon and 
planets we suggest that differential corrections to the ephemerides based 
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on the current system be tabulated until such time as new or revised 
theories have been completed. 

We intend to meet again at the beginning of the 1964 Assembly so that 
we may consider any fresh observational evidence or theoretical argu- 
ments that may have been brought to our attention since our meeting in 
January 1964. We will then confirm or amend as necessary the list of 
constants given above, and this final list will be our recommendation for 
the “lAU System of Astronomical Constants.” We therefore request 
that the Executive Committee submit the following draft resolution for 
consideration by Commissions 4, 7, 8, 19, 20 and 31, with a view to its 
adoption by the General Assembly : 

The International Astronomical Union endorses the final list of constants prepared 
by the Working Group on the System of Astronomical Constants and recommends 
that it be used in the national and international astronomical ephemerides at the 
earliest practicable date. 
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Determination oj Relative Locations of Various Areas 

of the Larth 


Traditionally, geodesy has been concerned primarily 
with mapping. An understanding of the terrestrial 
gravity field has also been of interest, but largely because 
of the necessity of understanding the gravity field 
in order to interpret observations of astronomic or 
geodetic positioning. That satellite geodesy has re- 
versed this emphasis is due primarily to the relative 
ease of interpreting satellite orbits in terms of the 
terrestrial gravity field. 

Satellites also provide excellent targets for geodetic 
triangulation either through simultaneous observations 
or through the use of predicted orbits to connect ob- 
servations separated by short periods of time. Both 
methods require accurate tirning, good coordination 
between observing sites and, to avoid systematic errors, 
comparable or identical equipment at each site to be 
connected. Good calibration procedures are also im- 
portant and for orbit interpolation an accurate knowl- 
edge of the satellite motion between observation times 
is needed. For these reasons, work in this area is just 
starting. 

The use of Baker-Nunn observations to improve 
station locations is reported in the paper by Izsak 
included earlier. The Echo satellites have been very 
easy to observe because of their brightness. The U.S. 
•Coast and Geodetic Survey has used observations of 
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these satellites for triangulation on the North American 
continent and as a planning tool for a worldwide net; 
work. For the latter, a much higher satellite will be 
needed, however. The following papers include a 
report on the use of preliminary observations of Echo I 
by French geodesists to tie France and North Africa to 
the same reference system, a report on the observation 
of a flashing light satellite for geodetic triangulation, 
and the results of a test of the SECOR electronic 
ranging system for such interconnections. The SECOR 
system has the added advantage of providing scale as 
well. 



The following paper is excerpted from a translation of La 
. Jonction G6od6sique France Afrique du Xord par Photo- 

graphies Synchrones du Satellite Echo I. Symposium de 
GeodMe par Satelliies. Institut Geographique National, 2° 

Direction — Geodesic, Group d’Etudes Spaliales, December 
1964. 

W66 37353 

• • 

Geodetic Junction oj France and North Ajrica 
by Synchronized Photographs Taken From 
Echo I Satellite 

H. M. Dufour 

Institut Geographique National 

I. GENERAL PRINCIPLES OF TRI ANGULATION BY SATELLITE 

T he basic principle is quite simple: a luminous object is seen 
against a different background of stars from different points on the 
Earth. Photography of this luminous object, together with photography 
of the stars which surround it, gives information which permits determin- 
ing the position of points on the Earth’s surface, other points presumably 
being known. 

To refine this idea: the solid Earth, considered nondeformable, to which 
we attach a trihedron TXYZ, has a known motion on the stellar sphere as 
a function of time, a motion given by positional astronomy; knowing the 
time means fixing the trihedron in the stellar field; each star thus consti- 
tutes an absolute direction in the trihedron TXYZ; a direction given by 
the stellar ephemerides. 

In a similar fashion, we must fix the luminous object, which, practically 
speaking, requires making synchronous (or quasi-synchronous) observa- 
tions from the different points on the globe which we want to tie together. 

The luminous object can emit flashes of light, thus automatically as- 
suring simultaneity, or it may be continuously luminous, making it 
necessaiy to create “artificial flashes” by using synchronously revolving 
shutters at the various observation stations. 

Anna was such a flashing satellite : it has been useful for geodetic ex- 
periments. It was more costly than a continuously luminous saielhtc , 
it worked only on command, and it is no longer working. 

Echo I and Echo II are continuously luminous satellites and Echo I 
was chosen for the France-North Africa junction. 
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Practically speaking, photography involves: 

— A certain number of star photographs to calculate the orientation of 
the camera. 

— A certain number of light flashes. In the case of the French observa- 
tions, the shutters turned at the rate of 60 revolutions per minute, 
producing 60 luminous points per minute, which we call “flashes.” 

Essentially the problem is to interpolate the position of the flashes 
among the known star positions; let us mention immediately that the 
60 flashes are reduced to 1 central flash, and that the direction for this 
central flash is obtained in the trihedron TXYZ. 

Thus each photograph gives a known direction SF, S being the station 
and F the fixed point. We shall see how the geodetic problem is treated 
for the France-North Africa junction, starting from this fundamental 
principle. 

Several results are necessary for precision: 

If an angular precision of 1 second of arc (that is, 1/200 000 of a radian) 
is desired for SF, we must have 

— A 1 -millisecond precision tolerance for the revolving shutter; 

— Accuracy of sidereal time to one-twentieth of a second ; 

— Accuracy within 2 or 3 microns for the determination of the positions 
of the stars and of the flashes on the photographic plate (for a focal 
length of 30 centimeters). 

II. INSTRUMENTS 

The apparatus includes, essentially 

— A camera of 30-centimetcr focal length, furnished with a revolving 
shutter and with a leaf shutter. 

— A quartz chronometer, to regulate the shutter at the rate of one rps, 
for photography of the satellite; and to release the leaf shutter for 
the star photographs. 

The entire apparatus is extremely simple and highly portable. 

Notation of time is done 

— For the satellite, by direct reading on the revolving shutter scale: 
each second blip of the time signal illuminates the frame of the 
shutter, always at the same point, which permits the observer to take 
a correct reading (to almost 1 millisecond). 

— For the stars: according to the same principle — on a revolving gradu- 
ated disk incorporated into the quartz chronometer. 

III. THE FRANCE-NORTH AFRICA GEODETIC TIE (MAY 1964) 

Under the auspices of the National Center for Space Studies (CNES), 
the National Geographic Institute was charged with effecting the France- 
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Figure I.— Geodetic junction, France-North Africa, by Echo I, May 1964. Points 
of the sky photographed, with an indication of successful stations. ► suc- 
cessful photograph. + possibly usable photograph. The direction of 

the arrow identifies the station. P: Paris Observatory; B: Besancon Observatory. 


North Africa tie by means of synchronous photographic observations of 
the Echo I satellite. (See fig. 1.) 

Known points : 

Lacanau (near Bordeaux, on the Atlantic) 

Agde (near S^te, on the Mediterranean) 

Oletta (Corsica) 

Points to be determined: 

Hammaguir (West Sahara) 

Ouargla (East Sahara) 

Actually, the five points are already known in the European Compensa- 
tion system (called European Compensation 1950), and the work con- 
sists of comparing classical geodetic transfer with transfer by space 
geodesy. 
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Satellites Used 

The only observations retained are those made of Echo I, which, by 
reason of its culmination at around 47.5°, is well suited to observation 
from the Mediterranean area. 

The observations took place from the 4th to the 26th of May 1964. 

It was sometimes possible to obtain four useful passes in one evening. 

Echo II was similarly observed beginning May 20, but no observations 
have been retained. Echo II could have at most only one useful passage. 

Geodetic Configuration 

In the final calculations, all the observations of absolute directions will 
be replaced by observational relations, which will permit definition of the 
unknown points (with their table of errors) ; but for the observations 
themselves, it is necessary to define a general line of work which assures 
a priori a satisfactory overall configuration. 

One may reason as follows : 

First Point of View 

The central meridian of the observation zone has a longitude of about 
2° East; the central parallel is at 38° latitude. The most important 
points for the tie are the points of the satellite’s trajectory, on the mean 
parallel, to the West and to the East of the central meridian, respectively 
(the eastern point is in the zenith in the south of Italy, the western points 
are in the zenith in Spain). These points are determined by intersection 
from the three French bases. Their intersection determines the African 
points. Considering the diagram, it becomes apparent that this latter 
intersection is a good one, the first one not being as desirable, since the 
French points constitute a rather narrow base. 

Second Point of View 

The set of directions defined with reference to the stars determine a 
certain number of planes; these planes, combined, define the lines joining 
the points on the ground. 

Examining the diagram, one realizes that all the lines are thus well 
defined, except the line Ouargla-Hammaguir, as it was difficult to have 
African points outside the shadow cone. Here again, one notes that the 
African points are defined by an intersection starting from a French base 
which is rather narrow. 

Forecast of the Passages 

For each passage of the satellite, the aiming point on the satellite s 
trajectory is, in practice, imposed; the coordinator of the work chooses 
this point on the basis of various elements. 
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• ^ — First of all, is the necessity for a generally well formed geometric 
figure, as we have noted. 

— Observation time suitable for star photography: from this point of 
view, the stations at high latitudes are very unsuitable in summer. 
— It is also necessary that the satellite leave the shadow cone soon 
enough that the observers are not surprised (in practice, the time of 
the satellite’s arrival often cannot be predicted to better than about 
5 minutes). 

Theoretically, it should be possible to take into account previous success- 
ful evenings of observation, but practically, this is hardly possible ; in fact, 
the general coordinator establishes a plan and follows it more or less rigidly. 

For the predictions, the Echo I ephemerides furnished by the Smith- 
sonian Astrophysical Observatory in Massachusetts are used. 

Upon reception of an ephemeris, the passages are numbered, the 
ephemeris is extrapolated for about 15 days, and the points of passage 
are chosen (that is, the points in the sky at which the ballistic cameras 
will have to be aimed) . An electronic program then furnishes a statement 
indicating for each station all the elements of the observations. 

In practice, between the elements of the extrapolated ephemeris and 
the values from the following ephemeris, we currently encounter differ- 
ences on the order of 10 minutes in time and 2° to 3° in longitude. These 
differences, although disagreeable, are acceptable for the observ’ations. 
They underscore the fact that Echo I, being very light, is extremely 
sensitive to radiation pressure. 

Field Operations 

The task of taking the pictures requires only two operators, but in 
practice it is a good idea to provide a team of three, since, during the 
usable time p>eriod, photographs are taken every 2 hours during the night, 
every night. It is desirable under these conditions to provide one night’s 
rest every three nights for each operator. 

The picture-taking operation includes the following phases: 

— Placement of the camera in the station 2 hours before the observa- 
tion, and starting the chronometer 
— Loading a plate in the camera 
— Setting the chronometer, by time signals 
— Photography of the stars (first exposure) 

— Starting the shutter disk ; regulation by time signals 
— Photography of the satellite, at the rate of one exposiire per second, 
except for the fiOth-second blip of each minute, which is eliminated 
by a disk manipulated by an operator^ 

•The operator hears the time signals and places the disk in front of the objective 
when he hears the long signal of the full minute. 
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— Stopping the disk: photography of the stars (second exposure) 

— Photography of the synchronizing marks at the base of the camera 
— Replacing the plate 

— Setting up the camera for the next observation 

Operations also include: 

— Development of the plate 

— Recording the following information: direction of the satellite’s 
passage, notation of the minutes . . . direction of diurnal movement. 

The important elements of field operations are essentially 
— The stability of the apparatus during the interval of time which 
separates the two photographs of the stars ; this interval is presently 
15 minutes, and will be reduced to less than 10 minutes. 

— Recording the time on the revolving shutter, which must be done 
with an accuracy within 1 millisecond. 

Success of the Observations 

About 60 positions of the satellite were selected in advance to be photo- 
graphed (each comprising 60 points). Here is the tabulation of successful 
photographs : 

Hammaguir: 15 Ouargla: 40 

Oletta: 40 Agde: 35 

Lacanau: 25 

Naturally, every point photographed simultaneously by at least two 
stations is usable. Theoretically, the probability of one position in five 
successful negatives is: 

1 5 V/ 4 0 \y 4 0 sy 3 5 2 5 1__ 

(In fact, there were only two entirely successful points out of 60.) 

IV. PROCESSING THE PLATES— CALCULATIONS 

We cannot go into the complete details of this last phase of the work 
here. We only indicate the fundamental subdivisions, which involve the 
following operations : 

(A) Calculations of the theoretical coordinates of the stars on the plates 

(B) Measurement of the plates with a comparator 

(C) Calculations of the formulas giving the direction cosines of the 
flashes in the terrestrial Cartesian system 

(D) Intersection of the flashes from the approximate coordinates of the 
stations 

(E) Calculations of the unknown coordinates of the stations from the 
known coordinates. 
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The characteristics of these various operations are given in the follow- 
ing pages. All the calculations were made on the CAB 500 of the 
Geographic Institute. 

Other Calculations 

In addition, the coordinates of the tie points in the terrestrial Cartesian 
system were carefully established by means of classical geodesy. 

Planar coordinates, the geographic coordinates (X<^) for North Africa in 
the European Compensation system, are available; their accuracy is 
close to 10 meters. 

Altitudes, the altitudes H of the different points above the International 
Ellipsoid, are calculated by integration of the deviations of the vertical 
at the astronomic points. 

First calculation: French network 

Second calculation: Southern France, Italy, Spain, North Africa 

The French network is dense enough to make the results worthy of 
confidence. The second network is made up of points rather distant 
from each other. After cross-checking, it may be estimated that 1- 
decameter precision is preserved here as well. 

The calculation is carried out by astronomic leveling, establishing the 
relations by observations between each point and its nearest neighboring 
points. We assume that each difference of level has an error independent 
of the others and proportional to the distance ; these hypotheses are not 
absolutely rigorous. It is proposed to include calculations statistically 
more valuable, which are presently being undertaken. 

The approximate coordinates of the stations (calculations D and E) are 
furnished by the terrestrial Cartesian coordinates, derived from the 
coordinates calculated on the International Ellipsoid. 

(A) Calculation of theoretical coordinates of the stars on the plates 

Ecliptic coordinates 1950; XeYeZe 



Selection of stars: 

Calculation of XYZ ecliptic coordinates of the camera for 1950 
Rejection of stars by studying the scalar product : 


XXe+YYe+ZZe 
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Reduction of the selected stars to (p, q, r) direction cosines in the - 
terrestrial Cartesian system 
Transformation into the local horizontal system: 



Refraction 

Transformation into the plate system : 



(/= focal distance) 


Output 

Tape Printed form 

i i 

Operation C Operation B 

zenith distance 
azimuth 

universal time for the stars 
universal time for the flashes 
pressure 
temperature) 

(B) Measurement of the plates with the comparator 

The measurement of a plate requires, in practice, a day’s work for two 
operators. It must be done very carefully. 

The operators use the approximate coordinates (Xcijc), given by calcu- 
lation A, to position the plate under the comparator. 

Then they note the elements of the plate in the following order: 

— Synchronizing marks at the bottom 
of the camera 
— Stars 
— Flashes 

— Flashes 
— Stars 

— Synchronizing marks at the bottom 
of the camera 


Operator 1, operator 2 
acting as secretary 


Operator 2, operator 1 
acting as secretary 


DZ: 

AZ: 

UT: 

UT: 

P: 

T-. 
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The use of two operators has the following advantages : 

. — Systematic obser\’ational errors show up in the course of the work 
itself ; sometimes it is possible to correct these on the spot 
— The secretary figures the averages in the course of the work ; there- 
fore, any error is immediately corrected 
— The operators keep each other company (a rather boring task) 
and they can rest their eyes 50 percent of the time 

Following the comparator measurements, the flashes are smoothed; 
instead of 60 flashes in the minute chosen, only 5 flashes are kept, corre- 
sponding to seconds 10, 20, 30, 40, and 50. The smoothing is based on 
the hypothesis that the x and y coordinates on the plate are represented 
by polynomials of the fourth degree as a function of time. 

(C) Processing the plates 



Use of the stars : Xoyo — /: 

Correction of radial distortion 

Examination of the homographic formula, in least squares 



Application to the flashes ; 

Correction of distortion 

Application of the homographic formula 


>(Xc2/c-/) 

Return to terrestrial Cartesian system: 

Complementary corrections : 

Time of propagation of the radio signals 
Correction to reduce the position to the even second 
Correction of position on the plate 
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Printing of a tape : 

V 

q of the flashes and elements P, T 
r for refraction 


(D) Intersection. Several simultaneous plates give a series : 


Plate 1 


Plate 2 


Plate 3 



First least squares calculation, without weighting, and without correction for 
refraction 

>Approximate coordinates of the flashes. 

Second calculation, with weighting, and with correction for refraction 

>Coordinates XfYfZf, the most coherent as a function of the coordinates of 

the stations. 

Printing of a tape 

i 

Operation E 


(E) Geodetic calculation 

Regrouping of the series: each plate itself gives two observational rela- 
tions in a general system of least squares. 

For these relations, the smoothed values of the central flash is used. 


For each plate: 



*In a more elaborate program, 10 relations will be written corresponding to the 5 
smoothed directions (seconds 10, 20, 30, 40, 50) with their table of variances. 
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For each flash: 

Refraction correction 

Phase correction 

Estabhshment of the observational relation 

Normalization 

Elimination of the flash coordinates 
Summation of the equations : 

(1) resulting from a priori connections between the station 

points (known points, known distances, etc.) 

(2) resulting from the normal equations, after elimination of 

the flashes 

Table N of weights of the coordinates 

Table of the variance. Calculation of mean quad- 
ratic error of a direction 

Coordinates of the stations 

V. RESULTS OBTAINED AND PLANNED IMPROVEMENTS 

In a general way, the Geographic Institute has been occupied with the 
production of highly transportable equipment, and with effective expe- 
rience; we hope in the future to perfect the various op>erations progres- 
sively, according to the successive problems which increased precision 
poses. 

At present, we expect from the France-North Africa experience a 
theoretical angular precision of 1/100 000 of a degree (that is, two arc- 
seconds) which corresponds, in gross fashion, to a precision of 20 m be- 
tween the extreme geodetic points. It is not yet possible to say an 5 rthing 
definitive about the practical results, the first of which are just be ginning 
to appear (operation D).^ 

We hope in the future to increase the precision of the measurements 
as follows: 

— Regulation of the entire picture-taking operation by a quartz 
chronometer; currently, this clock assures regular rotation of the 
shutter disk, but not its time setting, which is done directly from the 
time signals. In the future, the chronometer will also set this time, 
being regulated itself, by different time signals, before and after 
photography of the satellite and of the stars. This latter operation 
may thus be cut down to less than 10 minutes. 

®The sightings of series 6 (longitude 13° E, latitude 40°) intersect with a radial 
uncertainty of 10 m, an excellent result if one considers that the mean distances are 
2000 km. 
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— To use thick photographic plates (6 mm), rigorously flattened. 

— To use fast emulsions, in order to increase the number of stars 
recorded in one session. 

— To take two exposures of the stars before and two exposures after- 
ward. 

— To place two cameras in parallel in each station systematically: 
this procedure improves the precision of the directions in the ratio 
\/2 and permits the easy disclosure of any discrepancy. 

— To air-condition, at least partially, the picture-taking apparatus 
and the quartz clock. 
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Th£, folhwing paper is excerpted from Measurements of t£e * ** 
Doppler Shift in Satellite Transmissions and Their Use in 
Geometrical Geodesy. Presented at the European Sympo- 
sium on Satellite Geodesy, Paris, December 1964. 


Geometric Geodesy hy Use of Doppler Data 

Robert R. Newton 
Applied Physics Laboratory 
Johns Hopkins University 

« 

I T SHOULD NOW EE POINTED OUT that Only two of the coordinates of a 
station can be obtained accurately from the data obtained during 
one pass. In order to obtain all three coordinates, it is necessary to 
use data from at least two passes. 

By far the dominant error made in a set of station coordinates derived 
in this way, with our present state of knowledge of the earth’s gravity 
field, is in the position of the satellite. This satellite position error 
shows up one-for-one in station coordinates. Typical errors in satellite 
position currently run about 75 meters, which is an unacceptable level 
for geometrical geodesy. 

However, suppose that two stations observe the satellite at about the 
same time and that each one derives a set of station coordinates. Since 
the satellite is nearly in the same position as used by both observers, the 
errors in satellite position almost cancel out and the differences in the 
station coordinates derived in this way should be quite accurate. The 
most accurate case of this use of station coordinates to obtain geometric 
information is that in which both stations are exactly in the same loca- 
tion. This was the case with the information that we showed in the 
preceding section where we found that relative position could be found 
from a single pass with an accuracy of a few meters. 

As the distance between the stations increases, the accuracy of the 
geometrical position decreases for at least two reasons. First, neglected 
refraction effects, which cancelled exactly when the stations were to- 
gether, will no longer exactly cancel. This error is probably small com- 
pared with the fact that the two observers do not use the satellite at the 
same times and hence the errors in satelhte position are not exactly the 
same for the two observers. Analysis indicates that the error in relative 
position obtained this way should increase about as the square root of the 
distance between the stations. 
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The method of simply extracting station coordinates is equivalent to 
using data from one station to correct the position of the satellite and 
using this corrected satellite position to give the position of the other 
station. The computed satellite orbit, of course, will contain velocity 
errors as well as position errors and it is these velocity errors which cause 
the relative station position error to increase with the distance between 
the stations. The accuracy of relative information can therefore be in- 
creased if data from three or more stations at about the same time are 
used. In effect, what we can do then is to use information from two of 
the stations to correct both position and velocity errors in the satellite 
orbit and thus improve the accuracy of the third station relative to the 
other two. In practice, of course, one would probably treat the data 
from all three stations simultaneously, varying the orbit and the station 
coordinates until the best fit to all the data is found. The differences in 
station coordinates would then represent accurately the lengths of the 
chords joining the three stations. 

Analysis of doppler data by this method has been carried out by 
Anderle and Oesterwinter (Ref. 5) using data from six TRANET stations 
within the continental United States. The comparison of the set of 
chord lengths derived in this way and that derived from ground survey 
is shown in figure 1, based upon data from about thirty passes. The 



Figure L.— Comparison of chord lengths derived from ground survey and from 

doppler observations. 
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mjmber accompanying each chord in this figure is the discrepancy in 
meters between the two values of chord length. Figure 1 is copied from 
reference 1 with the p>ermission of the authors. 

Results of at least equivalent accuracy should be obtained using the 
integral doppler receiver, since studies of the TRANET stations show 
that they have about the same instrumental errors. 

Except for the chords that terminate at Station 003 in New Mexico, 
no discrepancy exceeds 20 meters, and most are considerably less. The 
comparatively large errors involving Station 003 caused the location of 
that station to be investigated again. It has been found that the station 
antenna was moved several times, in order to improve signal reception, 
without notification being sent to the analysis center. At the present 
time, we are trying to trace the time history of the antenna location to 
see if Station 003 results can be made consistent with other results. 

We have considerable confidence in the accuracy of the satellite results 
because of an incident connected with Station 710. The original results 
for this station showed an east-west discrepancy of about 30 meters. 
This discrepancy was traced to a card-punch error of exactly 1' in the 
longitude of the station. The results shown in figure 1 were obtained 
after correcting this card-punch error. 

In summary, the relative horizontal position of two stations wuthin a 
few hundred kilometers of each other can be obtained, using the doppler 
data from a single pass, with a standard deviation of about 10 or 15 
meters. By the use of more passes, all three coordinates of the relative 
position can be determined, and with a higher accuracy. Across the 
dimensions of a continent, using four or more observing stations, relative 
positions can be obtained in all three dimensions with a confidence level 
of about 10 meters, using a modest number of passes. 

REFERENCE 

1. Anderle, R. J. ; and C. Oesterwinter: “A Preliminary Potential for the Earth 
From Doppler Observations on Satellites,” U.S. Naval Weapons Laboratory, 
May 1963. 
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On October 31, 1962, the Department of Defense launched the first satellite 
designed primarily for geodesy. This permitted the testing of several optical 
and electronic techniques for precise geodetic satellite tracking. The follow- 
ing papers describe the satellite design and describe the results of observations 
of the flashing light which this satellite carried. Observations of the doppler 
frequency measurements are included in the results of doppler tracking 
reported earlier. 



N66 37355 

The following paper is a chapter from The Use gf AftiSciaW # 
Satellites for Geodesy, G. Veis, ed., North-Holland Pub- 
lishing Co. {Amsterdam), 196S, pp. 255-260. 


Project ANNA 

Mark M. Macomber 

Bureau of Naval Weapons, Department of the Navy 
Department of Defense, Washington, D.C., U.S.A. 


P ROJECT ANNA IS A UNITED STATES geodetic satellite program. The 
name ANNA is an acronym for Army, Navy, NASA (the National 
Aeronautics and Space Administration), and Air Force, the agencies that 
originally collaborated on formulating the program or that were expected 
to participate actively in the observation program. When the geodetic 
satellite program was first proposed, the Department of Defense feared 
that the output of the program, in addition to being of a scientific nature, 
might be of critical military significance. The program was therefore 
temporarily classified, and as a result, NASA participation did not 
materialize. Recently a review of the program and its expected results 
has been made, and the classification has been removed, opening the pro- 
gram to participation by NASA and, through them, by the world-wide 
scientific community. 

The belief that the scientific community should participate in ANNA 
was so firm to the people who formulated the program that in mid-1960 a 
provision was made to allocate 50% of the available light-flashes from the 
optical beacon to a scientific program administered by NASA. 

The first ANNA satellite is scheduled to be launched into a near-circu- 
lar orbit inclined 50° to the equator, with an altitude of one megameter. 
At present, only two launches are authorized, the second being scheduled 
as a backup to the first in order to ensure a successful orbit and to provide 
a sufficiently extended period of observation to produce significant data. 

The satellite contains three basic types of instrumentation that will be 
used for obtaining positioning information. For range determination, a 
transponder in the satellite is coupled with ground instrumentation that 
makes a phase comparison between a modulating frequency as transmitted 
to and as returned from the satellite. Three frequencies in the vhf-uhf 
bands are used: one for transmission to the satellite, and two coherent 
frequencies for transmission from the satellite to ground. Analysis of the 
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difference in phase shift on the two returning frequencies permits a cor- 
rection to be made for refraction effects. To resolve any ambiguities m 
distance measurements, four different modulating frequencies are em- 
ployed. The satellite transponder receiver is energized at all times, but 
the transmitters are turned off except when the satellite is being inter- 
rogated. Power limitations within the satellite restrict interrogation to 
six or seven passes per 24 hours ; the satellite is thus adequate for use by 
only one complex of ground stations. In an attempt to keep the satellite 
instrumentation as simple as possible, the ground stations are tied together 
by a VLF timing net, and transmit in bursts such that the signals from 
the various ground stations are received at the satellite sequentially. 
The highly complex ground instrumentation necessary to accomplish the 
interrogation and measurement functions, coupled with the limitation 
on power available for the transponder, make this system unavailable to 
the scientific community. 

For optical determination of satellite topocentric direction, a high- 
intensity optical beacon is used; when activated, it produces a series of 
five light-flashes spaced 5.6 seconds apart. The main advantage of the 
optical beacon is that when it is flashed, an infinite number of observers 
within the circle of visibility can observe the light with relatively simple 
equipment, making this system of maximum benefit to the scientific 
community. One drawback of the optical system is that the power 
consumption of the beacon is very large, limiting the number of flash 
sequences per day to about 20 or less, depending on the exposure of the 
satellite to sunlight to recharge the battery. Because of this power 
limitation, light flashes cannot bo programmed for everyone who might 
desire to participate in the program, but must be programmed to provide 
the maximum return in geodetic knowledge. The result is that certain 
observers will not be able to use the light when most convenient for them, 
although they will be able to observe it at other times even though the 
geometry of position solutions will not be optimum. 

Range rate information is obtained by observing the doppler shift of 
ultra-stable transmissions of the satellite. Four frequencies will be 
broadcast continuously for this purpose. The pair of frequencies 
designed for geodetic measurements will be 162 and 324 Me, with another 
pair, 54 and 216 Me, reserved for refraction studies and for a possible 
backup in case of failure of the prime tracking frequencies. All four of 
the frequencies are coherent, so tracking could be accomplished using 
any two. Since transmitters are of low power drain, they can be left on 
continuously, and thus be available to observers throughout the world. 
The ground instrumentation for receiving the signals and counting the 
doppler shift imparted to the signal is complex and may be a deterrent 
to observations by the scientific community. 
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The temperature-controlled crystals that drive the ultra-stable trans- 
mitters are also used to run a satellite-borne clock. Approximately 
every 90 seconds, time signals are broadcast from the satellite as phase 
modulation on two of the doppler signals. By keeping track of the fre- 
quency drift of the crystals, as obtained by the doppler tracking stations, 
it is possible to correct the time of transmission of these signals to better 
than 0.5 millisecond. Since the same clock that provides these timing 
signals also initiates the flash sequences for the optical beacon, the 
doppler and optical observations are tied together in time. 

The small satellite memory consists of 22 sixteen-bit words. Injected 
into it are the two’s complement of the identification of the 5.6-second 
timing pulse that should initiate each light-flash sequence. Every 5.6- 
second timing pulse generated by the satellite clock is added to each of 
the first 21 words, and when each word overflows a light-flash sequence 
is initiated. 

The satellite is oriented with the earth’s magnetic field so that the pole 
of the satellite visible from earth is a function of the geomagnetic latitude 
of the satellite. Two beacons face the north pole of the satellite, and two 
face the south pole. Only 15 bits of the 16 in each word in the memory 
are used for light-flash timing. The 16th bit has injected into it an 
indication of whether the north or south lights should be flashed. The 
22nd word of the memory will not initiate flash sequences, but is used 
for telemetry purposes to indicate whether or not the lights did flash 
as planned. 

In addition to the geodetic systems instrumentation, the satellite con- 
tains various minor experiments that test the environment and the 
attitude of the satellite. 

During the first three months after launch, an intensive calibration 
program will be undertaken, wherein the three types of measurements 
will be compared among themselves, and with terrestrial survey results. 
This calibration program is considered of prime importance to prove that 
no biases exist in any of the instrumentation used or in the methods of 
data handling that are utilized. If any biases do exist, they must be 
eliminated or corrected before any geodetic program can be undertaken 
on a world-wide scale. Upon successful completion of the calibration 
phase, a world-wide geodetic observation prograpi will be undertaken with 
the intent of refining our knowledge of the earth’s gravitational field and 
of providing the location of tracking stations relative to the earth’s 
center of mass. 

Now that the entire program has been declassified, NASA can partici- 
pate in the observation program, and can coordinate the participation of 
foreign and domestic scientific interests. Neither details of observation 
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equipment to be used, nor a modus operand! has yet been crystallized, 
but the following plans for participation are under active consideration.: 

1. use of the Smithsonian Astrophysical Observatory optical network; 

2. use of the MOTS (Minitrack Optical Tracking System) cameras 

located at most of the NASA minitrack stations ; 

3. use of minitrack interferometry information; 

4. use of the NASA-Jet Propulsion Laboratory radar systems for 

space probes; 

5. voluntary participation of interested observatories. 

Because of the power limitation for the optical beacon, only selected 
observatories that would add strength to the present planned deployment 
of tracking equipment can be considered for use in this program as far as 
special programming of the light is concerned. This does not in any way 
limit participation by all observatories within the area of visibility of each 
flash sequence, but it does mean that certain observatories will not derive 
the benefit of optimum geometry. All contact with civilian observers will 
be through the NASA. Details of this observation program will be 
promulgated to all interested persons as soon as they are made firm. 

At the present time, two satellites have been fabricated. They have 
completed comprehensive tests, performing perfectly all the while. These 
two satellites are currently undergoing extensive checkout after the 
shipment and are being fitted out to the launch vehicle. 

The activity collecting data will reduce and preprocess the data to 
eliminate any spurious data points, il the facilities are a\aihible. Since 
many independent observers will be unable to carry out this function, 
provisions must be made to handle reduction of data and elimination of 
spurious points at some central activity. Data will then be processed 
to produce a unique reference frame within the earth, to determine the 
earth’s gravitational field, and to provide geocentric locations ol tracking 
stations. 

The two planned launches, which together assure essentially one orbit 
of sufficient lifetime, are of course insufficient to meet the needs of any 
geodetic program. Some planning has taken place on the variety of 
orbits that would be of value to any program of this sort, and that are 
feasible from existing launch sites with existing boosters in the moderate 
price field. 
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Results From Satellite (ANNA) Geodesy 
Experiments 

Owen W. Williams, Paul H. Dishong, 

AND George Hadgigeorge 

Terrestrial Sciences Laboratory 

Air Force Cambridge Research Laboratories 


INTRODUCTION 

O N 31 October 1962 with the launch of the Satellite ANNA 1-B into 
a nearly circular orbit of approximately 600 nautical miles, a new 
era in the sciences of photogrammetry and geodesy was started (ref. 1). 
This geodetic satellite (designed for six months’ operation) continues to 
be reproductive some 24 months later. The name ANNA reflects its 
four United States sponsors, the Air Force, Navy, NASA and the Army. 
The Army was responsible for the electronic ranging system known as 
SECOR ; the Navy was responsible for the electronic doppler system and 
integration of the satellite payload; the Air Force was to develop the 
flashing light system; and NASA was to assist in the optical observm- 
tional program. This paper will cover only those aspects contributing 
to the optical phases of the active ANNA program. 

In January 1963, after more than two months of excellent p>erformance, 
a defective capacitor bank caused a malfunction in the optical system. 
The output of the light system was reduced to abotit 25% to 30% of its 
original value and this condition remained in effect into July. Therefore, 
between January and July only a limited program of flash transmissions 
transpired and while usable camera data was still received it was for the 
most part not up to original expectations. 

Prior to the capacitor problem in the satellite, geodetic stellar cameras 
obtained an excellent series of photographs. The image diameter of the 
flash recorded on 103-F emulsion averaged 70 microns on the PC-lOOO’s 
and approximately 50 microns on the BC-4-300 plates. In mid-July 
1963, the light output returned to normal and good data were again 
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obtained. Because of normal solar cell deterioration, only 7 flash 
sequences a day can currently be programmed as opposed to the original 
30 sequences. 

Due to the condition of the main battery controlling the command 
system and the lack of doppler tracking data, the Applied Physics Labora- 
tory, Johns Hopkins University, suggested to Air Force Cambridge 
Research Laboratories in October 1963 that we operate the flashing light 
by employing our alternate optical logic. The emergency over-ride 
system (EMOS) was developed to provide the necessary redundancy for 
the optical operations and this bit of foresight continues to pay dividends. 
The EMOS system consists of a World timing system accurate to one 
millisecond, a transmitter, a linear amplifier, and an antenna system. 

DESCRIPTION OF BEACON 

The ANNA optical beacon, developed for AFCRL by Edgerton, 
Germeshausen and Grier (EG&G), consists of two pairs of xenon filled 
stroboscopic lamps with reflectors, one pair on the north face of the solar 
cell panel and one pair on the south face (Figure 1). When cither set 



Figure 1 . — A close view of one of the four 
strobe lights located on the ANNA geo- 
detic satellite. 


of lights is triggered by the satellite memory (no longer possible) or l>y 
EMOS, a scries of 5 flashes is produced. An enlargement of a part ot 
one of the first photographs of ANN.V clearly demonstrates the above 
mentioned spacing (Fig. 2). 

LIGHT INTENSITY 

The light intensity of the ANNA beacons is not constant over the 
light angle but varies as shown by a dotted line in Figure 3 while the 
solid line shows the intensity used to predict the image sizes for ANNA. 
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Figure 2 . — ANNA I-B strobe light images. Copy of photograph taken during 
test to check optical system developed by Air Force Cambridge Research Lab- 
oratories. Satellite photographed at 0350 hr e.s.t. as it crossed Boston area 
with BC-4, 300-mm FL camera. Images on plate. 

VALIDITY OF EQUATIONS 

During the development phase of AXXA, the image sizes on Air Force 
PC-1000 geodetic stellar camera plates were predicted for conditions of 
moderate haze. Currently, two groups, EG&G and Duane Brown 
Associates, Inc., have e.xamined some of the AXXA PC-1000 plates to 
determine the validitj- of the equations. 

Duane Brown Associates analyzed 88 images on 20 different plates 
taken from six different PC-1000 cameras (ref. 2). The plates w^ere 
selected to be representative of different cameras, different locations, and 
a wide range of image diameters. The key results of the study are 
summarized in Table 1. Both groups obtained good agreement between 
computed and measured image sizes. There were four cases in which 
cameras 104 and 121 were employed in side-by-side operations. The mean 
image diameters from 104 turn out to be consistently and significantly 
larger than those from 121. This convincingly demonstrates that signifi- 
cant differences may exist in the capabilities of cameras of the same type. 

On the whole, agreement between theory and observation is considered 
to be sufficiently good (13.4 micron mean error) for the theory to be 

ncoH for r»nrr*OQ/a< nf apnprnl nlfirmiTiD* 

Q ^ ^ 

INTERVISIBLE OBSERVATIONS 

Simultaneous photographic observations were made by ten Air Force 
PC-1000 cameras during the period September 1963 to January 1964. 
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Figure 3. — Light beam pattern. 


Geodetic position determinations using these ANNA 1-B observations 
indicate that the PC-1000 camera system used with the intervisible 
technique is capable of extending geodetic control to a proportional 
accuracy of better than 1/100,000. 

Sixty observation nets were obtained in the Gulf Test. A net is 
defined as a flash sequence which was successfully photographed from 
three or more camera stations. Six nets were used to test the capabilities 
of the PC-1000 cameras. (See Table 2.) 

The 1381st Geodetic Survey Squadron of APCS provided the coordi- 
nates of the camera sites which were tied to first order control of the North 
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Table 1 . — Summary of Key Results of Image Diameter 

[RMS difference = 13. 4] 


RESULTS FROM SATELUTE GEODESY EXPERIMENTS 


Difference between 
observed and com- 
puted diameters 
(microns) 

— oot^coo 

I r-v-tlT-l 1 'T-H 1 

1 1 -1-^ 1 ‘ ^ 1 -b 1 ^ + 1 ' 

Standard error 
of individual 
image diameter 
(microns) 

5D'^G0WC^C^00*^’^C»?0b-0005C0C^C0OOCC 

cOro^^CDW^OOOWOciot^cDl>!:Do6cDTp 
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image 

diameter 

(microns) 
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Computed 
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Lf5if50COQ’-^COT}<tf540t^C5CJ:OC^-^i^5CD;D»— 

iCiCOiOOCOcOcDcOcDcOcDCOt^t>r«t^t>toOO 

Slant range 
normalized 
to zenith 
(miles) 

C0c0c0?0i0CD«D'^C05O'^4r5'^t»C000rt>OO«0 
QQaOOOC«GOt^'>^CCCOOOOCOt^CO’— ©OSOi'^ 
oo©©©©QOoooooocor«t«t«r«t«t>-©©cD 

Av. zenith 
distance 
(degrees) 

«CO^’-'Or-b-COXaOOOii^COC^©b-iOOOOOCO 

©0©©©u:3C»C^XCOOCC’-'Tt<iO»C«^iCir50 

ioiOTt<««t'^««i^Tj<'^cceocccococONo»c^i-H»-'i-' 
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(miles) 
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©©©©XXXOOXOOl>t^t^b>t*t«©©©© 

Camera 

No. 

2222^®2^22222 22222^ 

Plate 

No./(No. 

flashes) 
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American Datum (NAD 1927). ACIC supplied the geoid-spheroid 
^parations necessary to convert heights above mean sea level to heights 
above the reference spheroid. 



Figure 4.— The coordinates of station 640 are determined from those of 
stations 648 and 649 by nets 30 and 34. 

Table 3 . — Geodetic Position Determination of Station 640 From Stations 

648 and 649 

Reference stations 648 and 649 <rKs=<rB-w=<rH=0 

Input standard deviations station 640 trjv_s = 100 m, <rj_ «, = 100 m, = 5 m 

Input error station 640 0 = -|-3'; x = — 4' 

Observations: All stations observed 5 flashes each on nets 30 and 34 
Average observation <r used: 648, ."85; 649, i:08; 640, lf56 


North American Datum 27 
ANNA Data Reduction 

A 

a 


X 

ff (meters) 

R(meters) 

29°33'44.'80 

44.'78 

90°40'44.'19 

43r93 

7.0 

7.6 

7.1 

8.7 

.'02 

rii 

."26 

f25 

-.6 

4.8 

K = [(Ax)»+(Ai/)* + (A 2)*1 i/* 
Rr — [af 

Distance: 987 km 

Proportional Accuracy (NAD Standard) 1/140,000 
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GULF TEST REDUCTION 

The coordinates of station 640 (Table 3 and Figure 4) were determine^ 
from intervisible observations of 10 flashes in nets 30 and 34, using 
648-649 as a baseline. The mean observation standard deviations were 
0f85 for station 648, lf08 for station 649, and 1!'56 for station 640. 

Stations 648 and 649 were regarded as being perfectly known in 
location, whereas station 640 was given an input position error of +3" 
in latitude and —4" in longitude and standard deviations of 100 m in the 
North-South and East-West directions and 5 m in height above the 
spheroid. The sizeable station 640 position error was removed in the 
reduction with a resulting R; NAD-ANNA, of 7.1 m which corresponds 
to a proportional accuracy of 1/140,000. 

The overall reduction considered Station 647 as unknown and six 
surrounding stations (640, 641, 643, 648, 649, and 650M) as references 
(Table 4 and Figure 5). Again, the reference stations were held fixed 
and the unknown station had horizontal standard deviations of 100 m 

Table 4 . — Geodetic Position Determination oj Station 647 From Stations 
640, 641, 643, 648, 649 and 650M 


Reference stations 640, 641, 643, 

648, 649 and 650M ; aN-s=aE-w=<rH=9 

Input standard deviations stations 647; ass = 100 rn, oe-w = 100 m, o-h =5 m 

Observations; (1) All stations observed 5 flashes on Net 34, (2) Stations 641, 643, 
648 and 647 observed 3 flashes on Net 8 and 5 flashes on Net 22, (3) Stations 640, 
643, 648 and 647 observed 4 flashes on Net 19, (4) Stations 640, 643, 648, 649 and 
647 observed 2 flashes on Net 28, (5) Stations 640, 648, 649, 650M and 647 ob- 
served 5 flashes on Net 30. 

Average observation (T used; 640, 2IT1; 641, ITOl; 

643, r78; 648, f81; 

649, f99; 650M, IIW; 

647, r90 




X 

H (meters) 

R (meters) 

North American Datum 27 

30°14'48r28 

88°04'42"51 

5.2 


ANNA Data Reduction 

48H6 

42"38 

4.8 



A ri2 

.'13 

.4 

5.2 


<7 ro3 

ro3 

2.8 

3. 


R=[(^xY+{^yy+{^^yv'^ 
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and a vertical standard deviation of 5 m. The geometiy is favorable 
for a strong determination of the coordinates of the imknown station. 
Mean observation sigmas were as follows: 2'11 for 640; l.'Ol for 641- 
Of 78 for 643; 0f90 for 647; 0f81 for 648 ; Of 99 for 649 ; and 1 f07 for 650M. 
The resulting R of the reduction was 5.2 m and the R, was 3.2 m. Be- 
cause of the geometiy, proportional accuracy has no meaning in this case. 

LONG LINE AZIMUTH TEST 

Geodetic azimuths of 4 lines ranging from 264 km to 1365 km have 
l^n determined. These hnes were from the Gulf Test Nets. The test 
lines involve first order off-set stations 640, 641, 643 and 648. The 
azimuths between the stations were computed by AGIO ii.<;ing Rudoe’s 
method for normal section azimuths. 

RESULTS 

Table 5 shows the good agreements, ± 1 second of arc, that can be 
obtained. 
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Back azimuths between 648-641, 648-643 and 648-640 compute to 
essentially the same degree of accuracy. 

These results are encouraging considering that : 

a. Distribution of observations relative to the baselines is mediocre 
at best, being nonuniformly distributed about the line between the 
stations. 

b. The number of observations is very small. 

PLANS 

The remaining nets are now being analyzed. It is anticipated that 
all possible nets will eventually be combined in one overall reduction 
and published. 


Table 5.—ANNA-Determined Azimuths Compared to the Normal 
Section Computed Inverses 



Forward azimuth 

Line 

Distance (m) 

ANNA 

T.ni'Sh.?" 



1,365,165 

Normal Section 

TSMl'SV.S" 

641- 

-648 

Ad, 

0.8" 

1 



ANNA 

71°04'33.6" 

I 


950,628 

Normal Section 

71°04'34.5" 

640- 

-648 

Ad. 

0.9" 

1 



ANNA 

83°01'09.1" 

1 


1,080,640 

Normal Section 

83°01'10.1" 

643 

-648 

Ad. 

1.0" 





CONCLUSION 

The ANNA 1-B satellite is part of a research program, and it has 
proven the feasibility of using multi -angulation space surveying tech- 
niques to obtain geodetic data previously not obtainable with a sufficient 
degree of accuracy. The results of the Long Line Azimuth and the Gulf 
Test reductions demonstrate that the PC-1000 geodetic stellar camera 
system is operationally capable of extending geodetic control to a pro- 
portional accuracy of better than 1/100,000 when cameras in a network 
simultaneously observe a flashing satellite beacon. 
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SECOR for Satellite Geodesy 

T. J.. Hayes 

Topography and Military Engineering Office 
Chief of Engineers 

THE SECOR CONCEPT 

S ECOR IS AN ACRONYM formed from the method hy which the system 
operates: Sequential Collation of Range. The system is essentially 
an electronic distance measuring system designed to make use of the 
intermediary space positions of artificial satellites for deter mining 
geodetic positions on earth. The application of the geodetic data that 
is collected is based on the trilateration principle whereby three or more 
ground stations are placed at surveyed geographical locations of known 
geodetic coordinates. These individual locations can be as much as 
2300 miles apart, depending on the altitude of the satellite. A fourth 
station is placed at some geographic location for which the coordinates 
are desired. When the satellite appears above the radio horizon of the 
ground stations, each ground station sequentially ranges to the space- 
borne satellite on a time-shared basis in exact synchrqnism as estabhshed 
by the station selected to be the master timekeeper. 

The precise slant range between each of the four stations and the 
transponder in the satelhte is determined from phase comparisons of 
modulations transmitted to and returned from the transpKinder on radio- 
frequency electromagnetic waves. Each station records its range data 
on magnetic tape simultaneously with signals from an electronic clock 
(accurately set by means of Bureau of Standards timing signals) which 
permits fixing the location of the satellite in space. Since several posi- 
tions (at least three) of the satellite relative to the ground stations are 
then known, the position of the fourth station relative to the first three 
stations can be calculated. The recorded data are translated into a 
format suitable for use by an electronic computer located at the Army 
Map Service. Using basically a space resection computation from satel- 
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lite position triangles, the computer calculates the desired geodetic 
coordinates independent of the satellite’s orbital parameters. The 
accuracy of the calculations for the position location is enhanced by the 
vast data redundancy provided by the system. For a typical satellite 
pass of six minutes duration, for example, the total number of ranges 
measured by four ground stations is approximately 29,000. Because of 
considerations of geometry it is necessary to collect data during at least 
two satellite passes to determine the geodetic position of the unknown 
station. In practice, a larger number of passes are used to afford the 
best geometry. 

The computational procedure is repeated until the required precision 
criteria are met. When the unknown station’s geodetic coordinates are 
considered acceptable, this fourth station can then be designated as a 
station of known geographical location. This permits one of the other 
stations to leapfrog to another location whose position is desired. Thus, 
control will be extended by this progressive method. 

EQUIPMENT COMPRISING THE SECOR SYSTEM 

The new SECOR system has been developed by the Cubic Corporation 
under the guidance of the U.S. Army Engineers. The basic equipment 
for operation of the SECOR system consists of at least four identical 
ground stations and an earth-orbiting satellite carrying a SECOR 
transponder. 


SECOR GROUND STATION COMPONENTS 

The most complex part of the SECOR system is the ground station 
element. Each SECOR ground station consists of three air-transport- 
able shelters; The Radio Frequency (RF) Shelter, the Data Handling 
(DH) Shelter, and a Storage Shelter. Each ground station also carries 
its own support equipment which includes test equipment, generators, 
air conditioners, ranging antenna, communications antenna, single side 
band communication equipment, and wheel adapters for each shelter. 
All shelters and equipment are ruggedized for extensive and varied 
field operations. 

Major units in the RF Shelter include receivers, transmitters, power 
supplies, antenna control assembly, and intercommunications equipment. 
Although the antenna is normally mounted on the RF Shelter, it may be 
operated from a ground mount if required. After the satellite signal is 
acquired (usually from prediction data provided to each station ahead 
of time) the RF operator tracks the satellite by keeping the antenna 
pointed in the direction of maximum signal as shown by signal strength 
meters mounted on the antenna control panel. 
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The digitizing of the ranging and timing data is accomplished in the 
Data Handling (DH) Shelter. The major units in this shelter include 
the magnetic tape recorder, quick-look recorders, control console, ti min g 
assembly and master intercommunications. The data handling equip- 
ment controls the modulation and timing of the pulses transmitted by 
the RF equipment and converts the transponder replies into a form suit- 
able for recording on magnetic tape. The quick-look recorders, by 
presenting a continuous visual plot on moving paper tape, provide the 
operator wdth a constant check of system operation during a pass. 

The Storage Shelter contains radio equipment for communication 
with the other ground stations and with Field Headquarters. It also 
accommodates test equipment and working space for servicing electronic 
equipment. 

Since the SECOR system is planned for use in a worldwide, leap- 
frogging operation to tie continents, geodetic datum, and islands into 
one geodetic control net, the ground-based components have been de- 
signed to be compatible with several modes of transport. SECOR 
ground stations may be air transported by aircraft and helicopter. The 
stations also can be carried over water by ship and landed by landing 
craft. Overland travel is easily effected by towing. 

Each ground station and the satellite transponder comprise an elec- 
tronic distance measuring unit w'hich is largely automatic in operation. 
Distance measurements are made to the satellite-bome transponder by 
determining the round-trip phase shift of a series of modulated signals 
transmitted by the ground station and retransmitted by the transponder. 
The RF signal carrying the modulation frequencies is transmitted by the 
ground station to the satellite whose transponder demodulates these 
frequencies and retransmits the modulations on two offset carrier fre- 
quencies to the ground station. The shift in phase between outgoing 
and returned modulation signals is measured and then recorded on 
magnetic tape. By measuring the phase delay on several related 
modulation frequencies and augmenting these data with the time delay 
of a radar-like pulse, unambiguous measurements are obtained. The 
two offset carrier frequencies provide correction for ionospheric refraction, 
a problem common to all techniques using electronic magnetic waves. 

SECOR SPACEBORNE COMPONENTS 

The spacebome portion of the SECOR system consists of a SECOR 
transponder which may be mounted in its own separable satellite or 
installed as an integral part of a lai^er satellite. 

The SECOR transponder is composed of a receiver, transmitters, and 
power supply w’hich converts the satellite battery voltage to the voltages 
required. 
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In operation, the transponder remains in a standby condition until . 
activated by a select call signal generated by a SECOR ground station. 
Upon activation, the transponder receives and retransmits the ranging 
frequencies. This sequence is repeated for each ground station in turn, 
one complete cycle from the four ground stations occurring every 50 
milliseconds. The same ranging and carrier frequencies are used by all 
four stations. 

SECOR SATELLITE 

To attain a flexible launch capability required to carry out an opera- 
tional program, the U.S. Army Engineers dc\eloped a self-contained, 
separable SECOR satellite that could operate in its own orbit and would 
be of a size, shape and weiglit that would permit the satellite to ride as a 
hitchhiker on one of many available boosters. This resulted in the 
currently orbiting Type II SECOR satellite. It is 10 in. wide, 13 in. 
long, 9 in. high, and weighs about 40 pounds. In addition to the SECOR 
transponder, the satellite carries solar cells, batteries, a telemetry system, 
an antenna system, and a magnetic orientation device. 

OPERATING MODES 

Eigure 1 illustrates the operation of the SECOR system in its two 
principal modes of operation: simultaneous and orbital. The simulta- 
neous mode will be used except when geography forces the use of the 
orbital mode. 



Figure I. — Simultaneous and orbital tracking modes. 
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. In the simultaneous mode, range measurements are made to the 
satellite in rapid sequence (essentially simultaneously) from each of four 
ground stations during the period when the satellite is above the radio 
horizon of the ground stations. Three of the ground stations are located 
at known geodetic positions, and the fourth ground station is positioned 
at a site where geodetic coordinates are desired. Range measurements 
obtained in the simultaneous mode permit the determination of the 
desired geodetic coordinates independent of the satellite’s orbital param- 
eters. 

When the geographic locations are such that the satellite can be viewed 
simultaneously by the three known stations only, and the unknown 
station is able to view the satellite shortly before or after on the same 
pass, the orbital mode will be used. In this case an arc of two or more 
satellite orbits precisely determined by the three known stations wall be 
extrapolated over the vicinity of the imknown station. Observed ranges 
from the unknown station to satellite positions on the extrapolated 
sections permit the desired geodetic coordinates to be estabhshed. 

The line crossing mode is a special case of the simultaneous mode. 
Geography may force the spread of the stations in a less-than-optimum 
ground geometry with associated low-look angles and long distances. 
The simultaneous solution for the position of the unknowm station in the 
above case may be strengthened by the direct determination of the 
geodetic length of the lines from the unknowm station to the three knowm 
stations by incorporating line crossing techniques. The known stations 
accurately determine the height of the satellite at the time when the sum 
of the ranges from the unknown station to the satellite and from the 
satellite to a known station is a minimum. This provides an additional 
mathematical constraint on the geodetic solution. 

TESTING THE SECOR SYSTEM 

To insure that the SECOR system as developed would fulfill the 
requirements of this research and development task, tests were performed 
on the various units and subsystems which comprise the system, then 
upon the composite system as a whole. Essentially, these tests may be 
described as tests of the ground station equipment prior to satellite 
launch, tests on the satellite prior to launch, and system tests conducted 
both prior to and after launch. 

To insure that the satellite would carry out the mission for which it 
was designed and developed, it was subjected to a full program of tests 
covering all the parameters and conditions considered importaut to 
system operation. The results indicated that the satellite would perform 
quite favorably within the anticipated launch and space environments. 

The satellite tracking phase of this development was planned to eval- 
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uate the test data obtained in the different operating modes to enable a . 
determination of SECOR system capabilities and thus to minimize the 
training, operating, maintenance, and data processing problems inevita- 
bly encountered after a new system goes operational. At the time that 
the satellite tracking plan was prepared it was envisioned that most of 
the simultaneous mode test data would be obtained on a small (500-mile) 
quad, with ground stations based at Stillwater, Oklahoma; Las Cruces, 
New Mexico; Austin, Texas; and (the unknown site) Fort Carson, 
Colorado, and that the orbital mode test data would be acquired utilizing 
the small quad plus the station located at East Grand Forks, Minnesota. 

Orbital mode operations were conducted for a period of about three 
weeks on the small quad with the East Grand Forks station. The track- 
ing data were obtained in such a manner that the data could be evaluated 
in both the simultaneous and orbital modes of operation. The results of 
simultaneous data reduction were used to check the results of orbital 
data reduction, clearly indicating where the orbital data reduction was 
erroneous. This comparison was necessary due to the unprecedented 
accuracy required in the projection of the satellite trajectory. 

It is believed at this time that reduced accuracy in the orbital mode 
accrues from three primary sources; namely, 

(1) Relatively small orbit fitting spans. 

(2) System errors. 

(3) Base site survey errors. 

(4) Internal timing errors. 

Small fitting spans allow any error in the data to upset the vector 
fitting and give less accuracy in the injection vector determination; the 
forward prediction then deteriorates rapidly. System and survey biases 
give a slight misorientation of the injection vectors and, therefore, affect 
the forward predictions also. Timing error arises from lack of exact 
synchronization of clocks at the ground stations, and any time offset will 
mean that the predicted satellite positions and the measured ranges will 
not be coincident. Time synchronization problems are not encountered 
when the Geodetic SECOR equipment is used in the simultaneous mode. 

The station at Las Cruces, New Mexico, was then moved to Larson 
Air Force Base, Washington, to complete the large quad shown in Figure 
2. The longer ranges and lower elevation angles to the satellite obtained 
on this quad provide the difficult operating conditions more typical of 
those likely to be met in operational usage. Data obtained in the simul- 
taneous mode on this quad were also reduced and evaluated in the orbital 
and line crossing modes to compare the capabilities and limitations of 
each mode for given satellite/ground station geometries, etc. 


144 



SECOR FOR SATEI.T.ITE GEODESY 



In the large quad data reduction, Larson Air J^ orce Base, a-suiiigton, 
was the unknown station. 


[Table 1 lists the majority of the unknown station solutions 
processed for the large quad. SECOR survey differences result 
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from averaging a sequence of actual Geodetic SECOR solutions 
and subtracting the U.S. Coast and Geodetic Survey site coordinates 
from this average solution. The standard deviations shown per 
solution are computed from the residuals which are the differences 
between each SECOR-survey offset and the average of the offsets 
for one particular solution. RSS refers to “root sum square” and 
indicates composite bias and noise error. The RSS is computed by 
squaring the mean offset and the standard deviations, adding and 
taking the square root of the sum.] 

Utilizing generally the same line crossing technique currently employed 
by systems such as HIRAN (but in this case employing the satellite as a 
super-elevated spaceborne station in place of the HIRAN airborne 
station), line crossings were initiated on a comparatively short line (1150 
miles) between San Diego, California, and Stillwater, Oklahoma, for a 
period of about one week to resolve operational problems, determine the 
accuracy obtainable with the line crossing technique, and finalize the 
computer program. Geometry was selected to allow good satellite 
altitude determination by data obtained from three of the four stations 
on the small quad. Data from the line crossings then were used to 
determine satellite position and geodetic distances between San Diego 
and Stillwater. 

Line crossing operations were then switched to a long line (2300 miles) 
between San Diego and Herndon, Virginia, for a period of about three 
weeks to determine the ultimate operational and accuracy limitations of 
the SECOR system for line crossing operations. Measurement errors 
from tropospheric refraction, ionospheric refraction, and multipath effects 
increase as elevation look angles are decreased. The data gathered 
should determine the accuracy of SECOR measurements at the long 
ranges obtainable when using satellites for line crossings. 

Table 2 contains a summary of the line crossing solutions and com- 
parisons computed from the satellite data. The results of the line cross- 
ing mode are commensurate with the theoretical results with the exception 
that all but 2 of 17 lines measured to the Herndon, Virginia, site are 
longer than the inverse lengths. 

The line crossing solution in this experiment is made possible by the 
use of ground stations to track the satellite during the crossing. The 
line crossings processed here represent the longest lines ever measured 
and clearly demonstrate the potentiality of the technique. 

A simulation study was performed by the contractor, Cubic Corpora- 
tion, to aid in the evaluation of the different operational modes (simul- 
taneous, orbital and line crossing). Effects of base site geometry, errors 
in site survey, tropospheric and ionospheric refraction, scale factor, and 
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(1) Error model used: ffTropo=6%; <rs<^.i„= 1 ppm; <Tsurvey=4 ppm; <n„„„=5%. 

(2) Error model used; <rsy.= 15'; <ri,opo= 5%; ascaie= 1 ppm; <rBurv.y=4 ppm; <nono=5% 
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(1) Error model used; <r8y.= 9'; <rTtopo=5%; (r8cai«= 1 ppm; <rBorv.y=4 ppm; <riono=5%. 
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equipment limitations have been considered with respect to their' effect 
on solution accuracy. In the computer simulation, geometries were 
selected to correspond to the actual site locations used in satellite tracking 
tests. 

LINE CROSSING ERROR ANALYSIS 

The method employed to obtain geodetic line crossings from Geodetic 
SECOR ranging data is not directly analogous in computational tech- 
nique or operational procedure to that performed with HIRAN or 
SHIRAN observational data. A distinctly different approach is neces- 
sary because the path of the satellite is not essentially concentric with 
the earth’s surface. The occurrence of a minimum range sum derived 
simultaneously from two tracking sites will, therefore, not be coincident 
in time with the minimum geodetic line crossing. However, if the range 



Figure 3 . — Geodetic distance error vs. elevation angle. 
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obsenrations are projected onto the earth’s surface at each point along 
the satellite s orbit, then the sum of the geodetic distances from two sites 
•to the subtrace of the satellite will have a niinimiiTn coincidently and 
approximately equal to the minimum distance between the two tracking 
sites. 

In the modified version of the classic HIRAN-SHIRAN line crossing 
schemes, it is assumed that the heights above mean sea level of the end 
points of a line are known. It has also been assumed that the geocentric 
coordinates of the vehicle will be known too, either by direct solution 
from three simultaneously-tracking known base sites or by orbital pre- 
dictions. The distance of the satellite from earth center is computed 
from the geocentric coordinates. With the station heights, satelhte to 
earth center distance, and approximate knowledge of the geodetic coordi- 
nates of the end points of a line, the projections onto the reference sphe- 
roid of the range measurements to each site can be computed. The 
minimum geodetic distance is now computed from the composite set of 
adjacent projections. 

In the error analysis for the geodetic satellite fine crossing, the effects 
of three major factors on the geodetic distance computation are evalu- 
ated; namely, 

(1) Ranging errors due to sources such as refraction, scaling, system, 

etc. 

(2) Station height errors, which are errors in independent height- 

above-mean-sea-level measurements. 

(3) Vehicle height errors, or errors in earth center to satelhte distance. 


Table 3 . — Relative Contribution oj Error Sources to Geodetic Line 

Crossing Error 


Elevation 

angle 

Total ESS 

Vehicle 

height 

Hanging 

Station 

height 

10° 

18.3' 

7.7' 

16.3' 

2.6' 

20° 

17.7 

9.1 

14.2 

5.5 

30° 

19.8 

11.3 

13.8 

8.6 

40° 

24.0 

14.5 

14.7 

12.5 

50° 

30.8 


16.6 

17.6 

60° 

44.0 


21.4 



64.0 

40.5 

29.3 



Error Model 



^Tropo 

<^IoXK> 

<^Seale 

<^T*pfet 

trsite 

9' 

0.05 

0.05 

1 ppm 

15' 

15' 
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A set of geodetic distance versus elevation angle error curves' were _ 
processed using a composite typical error model and different station and 
vehicle height errors. Figure 3 illustrates these error propagations,. 
Table 3 shows the amount each major error source contributed to the 
total geodetic distance error for one case shown on Figure 3. A satellite 
height of approximately 500 NM was used in this propagation. 

The line crossing analysis demonstrates that ranging accuracy controls 
the geodetic distance error propagation at low elevation angles. When 
elevation angles exceed 30°, station height and vehicle height errors pre- 
dominate. In all cases, solutions should be constrained to lines with 
elevation angles less than 45°. 

All errors propagated in this evaluation are carried as random and un- 
correlated errors. The results are therefore conservative. In actual 
solutions, and in particular under the scheme used to compute the satel- 
lite geodetic line crossing, errors will be correlated and will give either 
equal or more accurate results. 

CONCLUSION 

The advent of artificial satellites has opened doors to technological 
advances that were only dreams less than a decade ago. In January 
1964 a United States Army Corps of Engineers artificial satellite carrying 
a SECOR transponder successfully achieved orbit. This event made 
possible a series of tests of this new geodetic electronic ranging system. 
The tests, which were completed on 1 May 1964, clearly established that 
the SECOR system has the full capability of a first order geodetic tool. 
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Summary and Conclusions 


S ATELLITE GEODESY is Conveniently considered as divided 
into gravitational or dynamical geodesy and geometrical 
geodesy; in practice, of course, there is some overlap. 

Dynamical geodesy comprises the study of the Earth’s 
gravitational field through the use of satellites as in situ probes. 
The techniques employed are the traditional techniques of 
celestial mechanics modified and updated to apply to close- 
Earth satellites in the presence of a perturbing nonspherical 
Earth, the Moon, and the Sun, and particularly adapted to 
rapid calculation with modern high-sp>eed electronic com- 
puters. 

Starting with the determination of the improved value of 
the flattening of the Earth in 1958, this field has matured 
until, in the past 18 months, consistent values have been de- 
rived for a number of zonal and tesseral harmonics by several 
analysis techniques and from several types of data. The 
lower terms, which are best determined, are summarized in 
the following table for the zonal harmonics, and in the table of 
tesseral harmonics (Pt. 3, paper by Wagner). By contrast, 
only the flattening value was well established before the 


Zonal Harmonics 



Kaula 

Kezai 

12 

1082.20X10-* 

1082.63 X10-* 

J3 

-2.57 

-2.56 


-2.01 

-1.63 

J5 

-0.066 

-0.185 

J6 

+0.324 

+0.593 

J7 

-0.465 

-0.376 


Anderle 

1082.65X10-* 

-2.59 

-1.53 

-0.165 

-t-0.793 

-0.426 
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satellite era and, as was mentioned above, this proved to be - 
significantly in error. Additional satellites in a variety of orbits 
will be needed before the dynamical analysis of the Earth’s 
gravity field reaches the point of diminishing returns, and 
particularly before we understand major differences between 
terms clerived by diverse methods. However, the quality of 
the data is already good enough to stimulate geophysical 
interpretations and attempts to correlate the satellite results 
with surface measurements. 

Geometrical satellite geodesy, that is, the use of satellites 
as reference points for determining precisely the relative loca- 
tions of various points on the surface of the Earth and partic- 
ularly for interconnecting datums and placing all observing 
stations and all individual datums on a common Earth- 
centered reference system, is in its infancy. Analyses of the 
Baker-Nunn observations of various satellites by reflected 
light, simultaneous photographs of the Echo satellite from 
widely separated portions of the Earth, analyses of Doppler 
data by various TRANET stations and analyses of the ranging 
data obtained using the SECOR system have all indicated 
that we are on the verge of significant progress in the area of 
geometric geodesy, but it will probably be at least 5 years 
before such observations contribute in a major way to map- 
making. Not only will placing the entire world on the same 
coordinate system be useful to geographers and navigators 
but it will also aid in space-probe tracking and in geophysics. 
At present, observations from various stations in NASA’s deep 
space net cannot be combined without introducing an addi- 
tional error from the uncertainty in the relative coordinates 
of the various stations. World gravity observations also lose 
some of their value because of the difficulty in interconnecting 
observing sites accurately. The geodetic satellite program 
should signficantly improve this situation. 

The launch of Geos in 1965 into an orbit specifically selected 
for dynamical geodesy should contribute appreciably to both 
geometrical and dynamical geodesy. Geos will carry a 
flashing light, a Doppler transmitter, two ranging systems. 
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and comer reflectors for use with laser transmitters. The 
•ranging systems and the laser reflections, if the latter prove to 
be a satisfactory op>erational system, should greatly improve 
the determination of scale a major problem in geometrical 
geodesy. The launch of the Pageos satellite, an Echo-type 
balloon, in a 3000-km orbit will provide an excellent target 
for triangulation across intercontinental distances and can 
serve as a basis for a worldwide geodetic network. Finally, 
for a detailed understanding of the mass distribution in the 
outer regions of the Earth, the gravitational field derived 
from satellite observations must be combined with local 
measurements. An appreciable extension of ground-based 
gravity surveys will be needed for this purpose, an extension 
which undoubtedly lies outside of the field of satellite geodesy. 

On a broader scale, both satellites and probes have given 
us better values of the various constants needed to compute 
ephemerides for astronomical bodies. The newly adopted 
values are given in appendix B. The selection and use of 
these values is discussed in Part 4. As the tracking of inter- 
planetary probes improves and as Ixmar and then planetary 
orbiters become operational, our knowledge of the shapes, 
masses, and gravity fields of the Moon, Venus, and Mars 
should be well established. This will, in addition to providing 
information on the present structure of these bodies, provide 
important boundary conditions for theories of the origin of 
the solar system. 
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SUMMARY 

r iR RECOMMENDED that the parametric values which are cur- 
rently most used in orbital computation be adopted as 
provisional standards, rather than those which may be the best 
available, because the “most used” values differ only slightly from 
the “best” values and further improvements in the values are 
expected within the next 4 years. Some of these values are: 

GMq =3.986032X10*0 cm* sec”*, 

Jt =1082.30X10-0, 

J, =-2.3x10-0. 

Ji =-1.8X10-0, 
a, =6,378,165.0 m. 

With parameters such as the foregoing the most serious geodetic 
errors affecting astronomy are tracking station positions. 

Standard methods of describing and transforming positions 
are suggested. 

INTRODUCTION 

This review recommends which geodetic parameters should be adopted 
as standard, the manner in which the parameters should be expressed, 
and the values which should be adopted. In making these recommenda- 
tions, current practice, available determinations, and anticipated im- 
provements will be considered. 

GRAVITATIONAL PARAMETERS 

For the notation of the earth potential, recommendations have already 
been made by Commission 7 on Celestial Mechanics, of the International 
Astronomical Union (Reference 1) : 

= ^ P," (sin cos sin (1) 
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where n = GM@, r is the distance from the center of the earth, R is the 
mean equatorial radius of the earth, is the associated Legendre poly- 
nomial, /3 is the latitude, and X is the longitude. Alternative notations . 
recommended for the gravitational coefficients are 

Jn=-C„.,. ( 2 ) 

and 



These two additions are suggested: 


1. Define 


(n+m ) ! 

(2n+l) {2-8J) 


where the Dirac delta 8m° is 1 for m = 0 and 0 for m^^O. The Cn.m, S„,n 
are coefficients of harmonics which have a mean square amplitude of 1 
for all values of n and m. 

2. Define the mean equatorial radius more precisely as the equatorial 
radius of the mean earth ellipsoid, i.e., the ellipsoid of revolution which 
best fits the geoid. This definition is consistent with geodetic practice 
and involves the equatorial radius with only two of the set of orthogonal 
parameters defining the radius vector of the geoid — the zeroth and second 
degree zonal harmonics. (The more literal definition of the mean equa- 
torial radius as the radius of the circle which best fits an equatorial section 
through the geoid would connect the radius to the infinite set of even 
degree zonal harmonics.) An alternative possibility for the equatorial 
radius in Equation 1 is the mean radius of the entire earth which, since it 
differs by a factor of 10”®, would affect the value of J 2 . The mean radius 
seems slightly preferable aesthetically, but current practice overwhelm- 
ingly favors the equatorial radius ; a perusal of some papers on close satel- 
lite dynamics and orbit analysis found ten workers using the equatorial 
radius but none using the mean radius (in addition, five theoreticians did 
not define their radius). 

To be consistent with the connection of equatorial radius to the mean 
earth ellipsoid, it is recommended that the following be the relationships 
between the astronomical parameters ti = GM® and —C 2.0 and the 
geodetic parameters R = Oe, the equatorial radius; 7e, the equatorial 
gravity; /, the flattening; and to, the rate of the earth’s rotation with 
respect to inertial space (References 2, 3, and 4) : 

GMe = o.®7«j^l+|m-/-^w/-~^m/®-0(f)j, (5) 
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‘where 

w*fl« 
m = — 

7« 

The v^ues of GM® and J 2 which are probably the most extensively used 
at orbit computation centers in the United States are (References 5 6 
and 7) : ’ ’ 

GM® =3.986032 ± 0.000030 X 10» cm» sec"®, 

J 2 =1082.30X10-*. 

In the alternative notation of Herrick, Baker, and Hilton (Reference 8) 
and Makemson, Baker, and Westrom (Reference 9) : 

k,= (GAfe)*^^= 0.019965049 megameter*/* sec~^. (9) 

The values of GM® and J 2 in Equation 8 are consistent with these 
values for the geodetic parameters: 



( 6 ) 

(7) 


Oe = 6,378,165.0 + 25.0 meters, 

7« = 978.0300 + 0.012 cm sec”*, 
f= 1/298.30, ’ 

£0 = 0.729211585X10-* sec-b 

The value for a, is a compromise between the solutions of Fischer (Refer- 
ence 10), and Kaula (Reference 11), and other values which are unpub- 
lished. The y, value differs from that of the International Formula and 
the Potsdam System (978.0490 cm sec-^) in three ways : 

1. Correction to Potsdam System absolute g (Reference 12) = 

-0.0128 + 0.0003; 

2. Change of flattening from 1/297 to 1/298.3 = —0.0051; 

3. Change of mean gravity over the earth’s surface (Reference 11) = 

-0.0005 + 0.0012. 

The correction to absolute g is a provisional value and has not been 
adopted by the International Union of Geodesy and Geophysics; an 
improved value should be forthcoming within the next few years from 
several determinations in progress (Reference 13). The correction to 
mean gravity is negative, mainly because correlation between gravity 
and topography was used to estimate anomalies for the areas without 
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observations, which are predominantly oceans. Solutions by Untila 
which fit observed gravimetry and do not use correlation with topography 
give positive corrections ranging from +0.0004 to +0.0019 cm sec . 
(Reference 14). Rather slow improvement is expected; problems m 
observing gravity at sea are not entirely solved (References 15 and 16). 
Somb improvement may also come from using the better statistical 

techniques which larger capacity computers permit. ,.c , i 

The value of GM® may also be obtained through the modified Kepler 
equation by using the radar mean distance of the moon A and the moon s 
mean motion n: 


GM®=- 


^ (i+^y 

1 +^ 

' I^E 


( 11 ) 


where (3 is the solar perturbation of the mean semimajor axis and 
is the ratio of the moon’s mass to the earth’s mass, equal to the lunar 
inequality (Reference 17). The most recently published value for ^ is 
384,402.0+1.2 km (Reference 18). As pointed out by Fischer, this 
value should perhaps be corrected because it is dependent on an exces- 
sively rounded-off lunar radius of 1740 km (Reference 19). The mean 
radius of the lunar limb is 1737.85 + 0.07 km. Geometrical determina- 
tions of the radius toward the earth vary considerably; Baldwins 
conclusion (Reference 20) leads to 1740.05 km, whereas Schrutka- 
Rechtenstamm (Reference 21) concludes that the bulge is too smaU to 
be determined. However, we are not interested in just the long axis ot 
a best-fitting triaxial ellipsoid, but rather in the mean radius of the area 
contributing to the leading edge of the radar return pulse, which would 
fall within the +7 degree area of libration. Contour maps of the moon 
(Reference 22, for example) indicate that the average radius of this ±7 
degree area could differ by as much as 2 km from the best-fitting ellipsoid. 
If the lunar surface is assumed to be an equipotential surface, then using 
the moments of inertia obtained from the physical libration yields 
1738 57 km as the radius toward the earth. Letting A =384,400.5+ . 
km, ^ = 0.0090678, n = 2.6616997 XlO"® sec~^ (Reference 23), and 
Mm/m£= 1/(81.375 + 0.026) (Reference 24) gives 

GM© = 3.986094 + 0.00004X10“ cm^ sec“^ (12) 

Using the 1/81.219 of Delano (Reference 25) reduces GM© to 

3.986001X10“ cm® sec"®, so the difference from solutions based on 
terrestrial data seems largely explicable as an error in the lunar inequality. 
The larger computers of today permit the application of more elaborate 
statistical techniques than it was possible to apply in 1950 (the year 
Delano and Rabe published their work). However, since the stellar 
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posifions are a major suspect for systematic error, it seems premature to 
reanalyze the Eros observations before the revised reference star systems 
.are available (Reference 26). Meanwhile, improved determination of 
the lunar inequality may be obtained from radio tracking of space 
probes such as Mariner II (1962 apl). Also, since spacecraft have been 
launched into high, nearly circular orbits such as those of Midas III 
(1961 c) and Midas IV (1961 a5), it may be worthwhOe to tiy to 
determine GMe from close satellite orbits. 

In addition to GMm and Jj, standard orbit computation programs 
usually incorporate and J 4 . The values which are probably most 
common at United States computation centers are (Reference 6 ) : 

J,= -2.3X10-«| 

J4=-1.8X10-«.j 

At present the best values of the zonal harmonics are undoubtedly those 
of Kozai (Reference 7 ) : 

J 2 = 1082.48 ± 0.06 X 10“', J* = — 2.562 ± 0.012 X 10“', 

^4= -1.84 + 0.08X10“', Js= -0.064 + 0.019X10“', 

Ji = 0.39 + 0. 12 X 10“', Jj = - 0.470 + 0.021 X 10“', 

Jg = - 0.02 + 0.02 X 10“', J9 = 0.117+0.025X 10“'. 

Note that the and J 4 now used, given in Equations 8 and 13, each 

differ from Kozai ’s improved values by less than 0.3X10“'; and that the 
coefficients J 5 and higher are all very small in absolute magnitude. 
Therefore, it does not seem worthwhile to adopt values, other than 
those already in general use, before 1966 or 1967, when analysis of 
geodetic satelhte orbits observed during the International Year of the 
Quiet Sun will be completed. 

Most of the current close satellite orbit analyses for geodetic purposes 
seek tesseral harmonic perturbations. In view of the smallness of these 
perturbations, it does not seem appropriate to adopt standardized values 
for the tesseral harmonics C,,„, The one exception might be 

C 2 . 2 , 82 , 2 , for which an upper limit would be useful because of its effect 
on supplemental energy requirements for 24 hour orbits. The most 
recent, unpublished determinations of Izsak, Kaula, Kozai, and Newton 
0.9X 10“' to 1.8X 10“' in amphtude (VC2.2^+S2.20 and from 
S' to 25" W in the direction of the principal axis [(1/2) tan“' (52.2/C2.2)]. 

GEOMETRICAL PARAMETERS 

As shown by analyses involving large systems of observations (Refer- 
ences 10, 11, and 19), the equatorial radius is a derived, rather than a 
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fundamental, quantity : accurate knowledge of the radius is not necessary 
to obtain other parameters, such as the lunar distance, geoid undulations, 
or datum positions by fitting of the astro-geodetic to the gravimetric, 
geoid However, for astronomical purposes, it is desirable to have a 
reference ellipsoid correct within +50 meters in order to obtain reason- 
ably correct positions of isolated tracking stations from astronomic 
latitude and longitude. Also it is convenient to have a unit of leng 
approximating the earth’s radius for use in the potential formula (Hqu - 
tion 1) and for use as a base line to compare or combine parallax observa- 
tions. For these astronomical purposes, the value of 6,378,165.0 meters 
given in Equation 10 should be entirely adequate. Marked iinprove- 
ment is not expected for about 5 years, by which time satellite observa- 
tions should contribute significantly to the strengthening of tnangulation 
systems and to the interconnection of geodetic datums. 

By far the most annoying problems in the astronomical application ol 
geodetic data pertain to tracking station positions. Errors in the 
adopted values of station positions, in conjunction with drag and non- 
uniform distribution of observations, prevent accurate determination o 
tesseral harmonics and are even believed to be a major cause of discrep- 
ancies in space probe trajectories (Reference 27) . These station position 
errors are due to both inadequate data and mistaken treatment ot data, 
in descending order of reprehensibility they include . 

1. Weak, erroneous, or nonexistent connection of tracking stations to 
local geodetic control (this includes the moving of antennas by stations 

without informing the computing center) ; , . 

2. Failure to state the datum or ellipsoid to which tracking station 

positions refer; 

3. Use of obsolete or erroneous standard datum and ellipsoid; 

4. An incomplete or ambiguous statement about how datum or 

ellipsoid transformations were made ; . u ■ . 

5. Failure to provide for geoid-ellipsoid difference in calculating heigh s, 

6. Neglecting systematic error due to incorrect observation ( or 
example, no Laplace stations) or incorrect adjustment (for example, 
arbitrary scale changes or rotations) of geodetic control connecting 
tracking stations more than, say, 1000 km apart; 

7. Actual observational error of position. 

In view of the number of geodetic datums and corrections thereto, they 
do not seem to be appropriate parameters to be adopted as standard y 
an international organization, except possibly for the large continental 
triangulation systems. The corrections to coordinates u, . w with posi- 
tive axes directed respectively toward latitude and longitude (0 , U ), 
(0°, 90° E), (90° N) obtained in the world geodetic system solution oi 
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Table 1 . — Corrections /o u, v, w From Reference 11 (Meters) 


Datum shift 

Au 

Av 

Aw 

WGS-NAD .. 

-23±26 
-57 + 23 
-89+40 

4-142 ±22 
-37+29 
4-551 ±53 

4-196±22 

-96±23 

4-710±40 

WGS-ED 

WGS-TD 



Kaula (Reference 11) are listed in Table 1, where NAD, ED,< and TD 
refer to the North American, European, and Tokyo datums, respectively. 
The uncertainties in this table are based on estimates of the errors due 
to interpolation and representation in the astro-geodetic and gravimetric 
geoids, and are probably a fair measure of item 7 on the above list, but 
may neglect significant effects falling under item 6. The relationships of 
the rectangular coordinates u, v, w to the geodetic latitude 4>, longitude X, 
and elevation h, referred to an ellipsoid of parameters a, and/, are: 

u=(v-\-h) cos<l>cos\, \ 

v=(p+h) cos<l>sm\, > (15) 

u?=[(l— e^) v-{-h\ sin <f),/ 

where »' = a,/(l— sin* <^)i« and e* = 2/— /*. 

To help minimize the number of unnecessary errors in categories 1 
through 5 on the above list, it is suggested that organizations be urged 
to pubhsh the following information pertaining to each tracking station 
for which they publish any precise observations of artificial satellites or 
probes, or orbital data based thereon : 

1. The names and coordinates of local geodetic control points, both 
horizontal and vertical, to which the tracking station is connected; 

2. The geodetic datum and ellipsoid to which the horizontal coordinates 
refer; 

3. The organization which estabhshed the local geodetic control points; 

4. The manner in which the horizontal and vertical survey connections 
were made from the local control points to the tracking station ; 

5. The date of the survey connection and a description of the termina- 
tion point of the survey; 

6. The geodetic (<^, X, h) and rectangular (u, v, w) coordinates of the 
station referred to the local geodetic datum: 

7. A statement of the geoid height, if any, estimated for the station 
and the basis for the estimate; 

8. If the tracking station position has been shifted for the purpose of 
referring observations (direction cosines or altitude and azimuth) or 


167 



SATELLITE GEODESY 


calculating orbits, the geodetic and rectangular coordinates after the 
shift and the ellipsoid to which the new coordinates refer. 

Every item on this list is an action which must be accomplished for 
any tracking station, but thus far the Smithsonian Astrophysical Institute 
Baker-Nunn camera network is the only one for which even part of the 
list has been published (Reference 28). It is symptomatic of the diffi- 
culties which occur that, since this publication, the coordinates for at 
least four of the twelve Baker-Nunn cameras have been found to be m 
error by 20 meters or more. These geometrical details of tracking station 
position are rather uninteresting, but they must be examined carefully 
and determined correctly if the full potentialities of modem tracking 
techniques are to be realized. 
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Appendix B 

Reference List of Recommended Constants 


ff'ECEOtNQ PAGE BLANK ROT RLttCD. 


Rejerence hist oj Recommended Constants 


Defining Constants 

1. Number of ephemeris seconds in 1 

tropical year (1900) s =31 555 925-9747 

2. Gaussian gra\dtationaI constant, defining 

the a.u k =0-017 202 09895 


Primary Constants 


3. Measure of 1 a.u. in meters A =149 600X10® 

4. Velocity of light in meters per second. . . c =299 792-5X10® 

5. Equatorial radius for Earth in meters .. . a, =6 378160 

6. Dynamical form-factor for Earth J 2 =0-001082 7 

7. Geocentric gravitational constant 

(units ; m® s“®) GE = 398 603 X 10® 

8. Ratio of the masses of the Moon and 

Earth ^ =1/81-30 

9. Sidereal mean motion of Moon in 

radians per second (1900) n^* =2-661 699 489X10-® 

10. General precession in longitude per 

tropical century (1900) p =5025*64 

11. Obliquity of the ecliptic (1900) t =23°27'08!26 

12. Constant of nutation (1900) N =9!210 


Auxiliary Constants and Factors 


fc/86400, for use when the unit of time is 1 

second A;' =1-990 983 675X10“^ 

Number of seconds of arc in 1 radian =206 264-806 

Factor for constant of aberration (note 15) . . . F, = 1 - 000 142 
Factor for mean distance of Moon (note 20) . . ^2 = 0-999 093 142 
Factor for parallactic inequality (note 23) . . F, = 49853!2 
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Derived Constants 

13. Solar parallax arcsin (aJA) = itq = 8" 79405 (S'! 794) ^ 

14 Light-time for unit distance il/c' = T.4=499.012 

• ^ =170-002 003 96 

15. Constant of aberration F ifc'TA = « = 20M958(20.196) 

16 Flattening factor for Earth /= 0-003 352 9 

= 1/298-25 

17. Heliocentric gravitational 

constant (units: m’s”^) AW^ = GS= 1^2 718X10 

18. Ratio of masses of Sun and 

Earth (GS)/(CE)=S/E = 332 958 

19. Ratio of masses of Sun and Earth 

+ Moon S/E(1 +m) = 328 912 

20. Perturbed mean distance of Moon, in meters 

FAGE{1 +ti)/n^ *2)4 = a, = 384 400 X 10^ 

21. Constant of sine parallax for 

• Moon Qe/a, = sin 7 t, = 3422" 451 

22. Constant of lunar inequality I^~A = (6.440) 

23. Constant of parallactic Fj ^~^ =P,= 124-986 

inequality 1+y.A 


System of Planetary Masses 
Reciprocal mass 


24. Mercury 6 000 000 

Venus 408 000 

Earth+Moon 329 390 

Mars 3 093 500 


Reciprocal mass 

Jupiter 1 047-355 

Saturn 3 501-6 

Uranus 22 869 

Neptune 19 314 

Pluto 360 000 
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